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ABSTRACT 


This  thesis  explores  the  possibility  of  designing  an 
electronically  steerable  phased  array  antenna  system  using  a 
LC  impedance  matching  network  and  sections  of  transmissiom 
line  as  an  interconnection  medium  between  pairs  of  radiating 
elements.  This  interconnection  network  would  be  used  to 
control  the  direction  of  the  array’s  radiation  pattern. 
Adjustments  of  the  parameters  of  the  LC  network  would  attempt 
to  enforce  a  desired  current  distribution  on  the  elements  of 
the  array,  resulting  in  the  desired  radiation  pattern. 

Two  design  procedures  have  been  investigated  and  are 
discussed  with  presentation  of  results  and  sample  radiation 
patterns  obtained.  The  first  design  method  uses  a  Taylor's 
series  expansion  to  linearize  the  network  equations 
describing  the  array.  The  other  design  method  utilizes  an 
optimization  routine  to  systematically  adjust  the  parameters 
of  the  impedance  matching  networks  until  the  desired  current 
distributions  are  realized  as  closely  as  possible. 
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I.  INTRODUCTION 


An  antenna  is  a  device  used  in  the  transmission 


and  reception  of  electromagnetic  energy.  Sometimes,  as  is 


the  case  with  most  broadcast  radio  stations,  transmission 


is  done  most  effectively  and  economically  with  a  single 


radiating  element.  However,  when  an  application  requires 


higher  gain,  directivity,  steerability,  or  other 


characteristics  that  a  single  element  cannot  provide,  an 


antenna  array  consisting  of  several  radiating  elements  may 


be  necessary  to  satisfy  one's  requirements.  In  the  field 


of  amateur  radio,  phased  array  antennas  are  used  quite 


frequently  to  provide  the  user  with  the  type  of  radiation 


pattern  desired.  Also,  the  military  has  used  large  antenna 


arrays  for  defense  against  ballistic  missiles  and  in  the 


surveillance  and  tracking  of  objects  in  space. 


Another  group  that  has  recently  become  interested 


in  the  use  of  phased  arrays  is  the  air  traffic  control 


community.  Since  World  War  II,  simple  mechanically  steered 


antennas  have  been  used  for  air  traffic  control,  but  with 


the  enormous  increase  in  the  number  of  takeoffs  and 


landings,  types  of  aircraft,  and  increase  in  their  speeds. 


the  potential  use  for  computer  controlled  phased  arrays  is 


quite  evident. 


In  an  antenna  array  consisting  of  a  single  type  of 


radiating  element,  there  are  generally  four  parameters 
that  can  be  adjusted:  the  number  of  radiating  elements, 
the  spatial  distribution  of  the  elements,  the  amplitude, 
and  the  phase  of  the  currents  used  to  excite  the  elements. 

When  considering  the  synthesis  of  a  phased  array 
antenna  system,  there  are  two  basic  requirements  that  must 
be  met.  First,  determination  of  the  current  distribution 
on  each  element  of  the  array  that  will  produce  the  desired 
radiation  pattern  must  be  done.  The  second  task,  which  is 
extensively  explored  in  this  thesis,  is  the  determination 
of  the  proper  impedance  transforming  and  phase  shifting 
networks  needed  to  produce  the  desired  current 
distribution  on  the  array.  Traditionally,  in  the  simple 
arrays  found  in  broadcast  and  amateur  radio,  the  impedance 
transforming  and  phase  shifting  medium  is  a  transmission 
line  connected  between  the  individual  array  element  and 
the  transmitter. 

In  this  thesis,  a  method  for  calculating  the 

current  distribution  for  a  desired  radiation  pattern  was 

used  as  a  starting  point  for  the  development  of  the 
.  .  2 

interconnection  networks.  This  procedure  allows  a  user 
to  specify  a  direction  of  maximum  radiation  (main  lobe) 
and  one  or  more  directions  with  a  minimum  amount  of 
radiation  (nulls).  This  method  of  beam  steering  was 
adapted  to  the  triangular  geometry  of  an  array  consisting 
of  three  equally  spaced  radiating  elements  as  shown  in 


»  ^  J 


Figure  1.  The  number  of  steerable  nulls  is  dependent  upon 
the  number  of  elements  in  the  array:  the  number  of  nulls 
is  equal  to  half  the  number  of  elements.  Therefore,  the 
radiation  patterns  in  this  thesis  will  be  specified  by  one 
null  direction  and  a  main  lobe  direction. 

As  an  example  of  the  beam  steering  procedure,  a 
radiation  pattern  and  the  corresponding  current  excitation 
on  the  elements  is  presented  in  Figure  2.  This  pattern 
was  generated  from  an  input  of  0m  =  degrees  for  the 
direction  of  maximum  radiation  and  0n  =  320  degrees  for 
the  null  direction.  The  radiation  pattern  shown  in  Figure 
3  was  generated  by  a  5  element  array  with  the  direction  of 
maximim  radiation  at  50  degrees  and  the  two  nulls 
specified  at  +180  degrees  and  270  degrees.3  Note  that 
additional  nulls  may  be  produced  in  unspecified 
directions . 

This  thesis  will  explore  an  unconventional  method 
of  enforcing  the  desired  current  distribution  on  the 
array.  The  elements  are  interconnected  in  pairs,  using  LC 
impedance  1  .tching  networks  and  short  sections  of 
transmission  line  as  an  interconnection  medium  between 
pairs  of  radiating  elements.  The  LC  impedance  matching 
network  is  in  <-!.e  form  of  a  PI  network  which  is  assumed 
lossless.  This  PI  necwork  has  capacitors  in  the  two 
parallel  branches  and  an  inductor  in  the  series  branch,  as 
shown  in  Figure  4.  Figure  5  shows  the  3  radiating 


Element  #3 


FIGURE  1  -  Three  Element  Triangular  Array 


ee  Element  Array  Pattern 


elements  and  the  proposed  interconnection  medium  which 
together  define  the  phased  array  antenna  system  to  be 
analyzed  and  designed  in  this  thesis. 


Two  independent  design  methods  for  determining  the 
parameters  of  the  matching  networks  have  been  developed. 
The  first  method  developed  is  based  on  solving  a  set  of 
nonlinear  design  equations  by  an  extension  of  Newton’s 
root-finding  method.  The  word  "design"  as  used  in  this 
context  means  that  the  current  distribution  on  the 
elements  are  given,  while  the  parameters  of  the  PI 
networks  are  the  unknown  values  that  need  to  be 
calculated.  Since  the  design  equations  are  nonlinear,  a 
Taylor’s  series  expansion  procedure  is  used  to  linearize 
the  set  of  complex  equations.  The  Taylor's  series 
expansion  is  performed  about  an  initial  guess  for  the 
nonlinear  variables  and  the  linearized  equations  are  then 
solved  using  matrix  inversion  techniques.  A  new 
approximation  to  the  desired  solution  results  from  the 
first  iteration  of  this  process;  this  new  approximation  is 
then  used  as  the  starting  point  for  the  next  iteration. 
This  process  will  be  referred  to  from  now  on  as  the 
iterative  matrix  solution  method. 

The  procedure  described  above  is  similar  to 
Newton's  Method  for  finding  a  root  of  f(x)=0,  where  an 
initial  guess  is  used  as  a  starting  point  for  what  one 
hopes  will  be  a  convergent  sequence  approaching  reasonable 


A  graphical 


values  for  the  root  in  question.4 
representation  of  Newton's  iterative  root-finding 
procedure  is  shown  in  Figure  6.  If  the  iterative  matrix 
solution  method  converges,  it  does  so  very  quickly,  as  is 
the  case  with  Newton’s  Method.  Newton's  Method  can 
diverge  if  the  initial  guess  is  not  fairly  close  to  an 
acceptable  solution  and  this  is  also  true  with  the 
iterative  matrix  solver  routine.  Considering  the  system 
of  equations  to  be  evaluated  Figure  6,  is  quite 
misleading;  it  misrepresents  the  complexity  of  the 
equations  actually  being  dealt  with  in  this  thesis.  The 
24  complex  equations  that  describe  the  phased  array  system 
of  this  thesis  are  presented  in  Chapter  Four. 

Consequently,  extensive  use  must  be  made  of  large-scale 
computational  techniques  and  resources. 

The  second  design  method  presented  in  this  thesis 
is  based  upon  the  quantitative  "analysis"  of  the  phased 
array  system  of  Figure  5  with  the  reactive  parameters  of 
the  PI  networks  being  given,  while  the  current 
distributions  on  the  array  become  the  unknowns  to  be 
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However,  since  a  particular  radiation  pattern  corresponds 
to  a  specific  set  of  relative  current  distributions,  a 
method  for  finding  the  parameters  of  the  PI  networks  that 
will  in  fact  produce  the  desired  current  distributions 
must  be  developed. 

A  direct  search  optimization  procedure  which 
adjusts  the  value  of  the  parameters  of  the  PI  networks 
until  the  element  currents  are  as  close  as  possible  to 
their  desired  values  was  developed  as  a  solution  to  this 
problem.  Initial  guesses  for  the  PI  network  parameters 
are  used  to  solve  the  system  of  complex  linear  equations 
describing  the  equilateral  triangular  array.  Three  of  the 
unknown  complex  variables  in  this  set  of  equations  are  the 
current  distributions  for  each  element.  The  difference 
between  these  current  distributions  and  the  given 
(desired)  current  distributions  becomes  the  criterion  upon 
which  the  optimization  procedure  bases  its  decisions. 

This  difference  is  then  squared  and  is  designated  the  hill 
height  magnitude.  The  procedure  is  a  type  of  feedback 
enforced  optimization  routine  where  the  present  values  of 
the  current  distribution  on  the  antenna  elements  are 
continually  fed  back  to  the  decision  module  of  the 
optimization  routine  in  order  for  the  ~I  network 
parameters  to  be  adjusted  accordingly. 

The  unique  direct  search  optimization  method  used 
in  the  constrained  "analysis"  design  method  of  Chapter 
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Seven  was  developed  in  this  thesis  as  a  possible  tool  for 
solving  the  nonlinear  and  linear  equations  describing  the 
equilateral  triangular  array.  It  was  developed  out  of 
discussions  with  the  author's  thesis  advisor  concerning 
multivariable  optimization  techniques.^  This  routine  will 
adjust  the  value  of  the  variables  of  any  function 
dependent  upon  the  magnitude  of  the  function.  Since  the 
function's  present  value  is  based  on  the  current  values 
(within  the  procedure)  of  the  variables,  these  variables 
will  progress  in  whatever  direction  will  lower  the 
magnitude  of  the  function.  The  routine  is  adaptive,*  it 
has  the  capability  of  increasing  and  decreasing  the  step 
size  of  the  variables.  This  gives  the  routine  the  power 
to  move  the  variables  in  small  or  large  increments.  Also, 
this  routine  is  able  to  detect  when  a  variable  is  close  to 
a  possible  solution.  A  variable  is  close  to  a  possible 
solution  when  the  routines  decision  module  begins  to 
periodically  adjust  the  variable's  magnitude  above  and 
below  the  particular  magnitude  that  has  (up  to  this  point 
in  the  process)  resulted  in  the  best  solution.  When  this 
oscillatory  behavior  is  detected,  the  magnitude  of  the 
changes  in  the  variable's  step  size  is  reduced. 

Early  development  of  the  first  design  method  was 
performed  in  APL  on  a  IBM  PC  computer,  but  once  the 
computation  burden  of  the  iterative  matrix  solver  problem 
went  beyond  the  capabilities  of  the  IBM  PC,  developement 
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of  the  procedure  was  moved  to  FORTRAN  code  on  a  VAX  11/780 
computer  system.  The  second  design  method  was  completely 
developed  in  FORTRAN  on  a  VAX  8600  cluster.  The  linear 
equation  solver  LEQ2C  in  the  standard  library  subroutines 
from  the  IMSL  scientific  subroutine  package  was  used  to 
solve  the  matrix  equations  of  the  analysis. 

In  Chapter  Two,  the  basic  principles  of  radiation 
patterns  are  discussed.  The  particular  characteristics 
used  to  determine  the  radiation  pattern  in  a  phased  array 
are  listed  and  reviewed.  The  interaction  of  the  radiating 
fields  from  each  element  in  the  array  is  demonstrated  and 
a  procedure  for  its  calculation  is  presented.  Chapter  Two 
is  completed  with  an  in  depth  evaluation  of  a  method  of 
beam  steering  which  allows  a  user  to  specify  directions  of 
maximum  and  minimum  radiation. 

The  methods  using  reactive  networks  for  impedance 
matching  an  antenna  driving  point  impedance  to  the 
characteristic  impedance  of  a  transmission  line  are 
discussed  in  Chapter  Three.  The  numerous  ways  of 
impedance  matching  are  listed,  as  well  the  factors  used  in 
deciding  which  method  is  best  suited  for  the  particular 
application.  Specific  examples  using  PI  and  TEE  impedance 
matching  networks  are  presented  with  calculations.  A  two 
element  vertical  phased  array  is  matched  using  PI  networks 
in  an  extensive  presentation  of  the  calculations  involved 
in  the  analysis  of  this  simple  array. 
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In  Chapter  Four,  a  detailed  description  of  the 
physical  and  electrical  characteristics  of  the  proposed 
phased  array  system  of  Figure  5  will  be  presented.  The 
equations  that  characterize  the  transmission  lines,  the  PI 
networks,  the  antenna  mutual  impedance  relations,  and 
Kirchoff's  Current  Law  are  reviewed.  The  quantitative 
analysis  of  these  24  equations,  which  result  in  a  solution 
for  the  currents  on  the  radiating  elements,  is  performed 
at  the  end  of  Chapter  Four . 

Presentation  of  a  interconnection  network  design 
procedure,  using  an  extension  of  Newton's  Method  for 
finding  the  roots  of  a  quadratic  equation,  will  be 
demonstrated  in  Chapter  Five.  Developement  of  the 
procedure  entails  linearizing  a  set  of  14  nonlinear  design 
equations  using  a  Taylor's  series  expansion  about  the 
initial  guesses  for  some  of  the  complex  unknown  variables. 
This  linearization  process  is  shown  in  detail  and  the 
resulting  linear  equations  are  listed.  This  chapter 
concludes  with  example  designs  and  the  patterns  produced 
by  those  designs. 

Chapter  Six  is  an  in-depth  discussion  of  the 
unique  direct  search  optimization  method  developed  by  the 
author  as  part  of  this  thesis.  The  fundamental  decision 
control  procedures  are  explained  using  a  quadratic 
function  of  two  variables  as  an  example.  This  method  has 
been  developed  into  a  n-variable  optimization  routine  in 
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the  form  of  a  FORTRAN  subroutine.  Example  functions  are 
shown  along  with  the  results  generated  in  attempting  to 
minimize  the  magnitude  of  their  hill  heights. 

The  direct  search  optimization  routine  is 
utilized  in  Chapter  Seven  where  it  is  applied  to  a 
constrained  "analysis"  design  technique.  The  design 
technique  of  Chapter  Seven  uses  the  analysis  equations 
that  describe  the  equilateral  phased  array  to  try  and 
enforce  a  desired  current  distribution  on  the  array 
elements.  Since  solving  the  analysis  equations  results  in 
a  current  distributions  on  the  elements  of  the  array,  this 
design  technique  will  use  the  optimization  routine  of 
Chapter  Six  to  adjust  the  choosen  PI  network  parameters 
until  the  desired  current  distribution  is  realized. 

Sample  runs  of  this  design  technique  and  the  resulting 
radiation  patterns  are  presented  and  discussed. 
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II.  BEAM  STEERING  IN  PHASED  ARRAY  ANTENNA  SYSTEMS 


The  electrical  field  produced  by  an  array  of 

radiating  elements  can  be  determined  by  vector  addition  of 

the  fields  produced  by  each  individual  element.  Phased 

array  antenna  systems  produce  a  pattern  with  directivity 

when  the  fields  from  the  elements  interfere  constructively 

in  the  direction  of  interest  and  cancel  each  other  out  in 

the  remaining  areas  (null  directions).  From  a  theoretical 

and  computational  standpoint,  this  can  be  done  very 

easily.  However,  from  a  practical  point  of  view,  the 

potential  performance  of  the  phased  array  can  only  be 

approached.  There  are  four  items  in  a  multi-element 

phased  array  that  determine  the  radiation  pattern  of  the 

6 

array  as  a  whole.  These  items  are: 


*  1. 

* 

* 

The  relative  displacements  between  each  of  the 
individual  elements. 

1  2* 

The  amplitude  of  the  current  flowing  into  each 
element . 

c 

\ 

The  phase  of  the  current  flowing  into  each 
element . 

jj  4. 

The  radiation  patterns  of  each  individual 
elements . 

In  this  thesis,  items  1  and  4  are  assumed  fixed, 
and  attention  is  focused  entirely  on  the  production  of 
feed  currents  with  appropriate  amplitudes  and  phases. 
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The  method  for  radiation  pattern  control  used  in 
this  thesis  allows  a  user  to  manually  input  the  pattern 
desired  by  specifying  directions  of  maximum  radiation  (0m) 
and  minimum  radiation  (@n),  and  receive  as  output  the 
complex  currents  required  on  each  element  to  produce  the 
desired  pattern.  Of  course,  even  if  the  proper  current 
distribution  is  on  the  array,  the  prospects  for  production 
of  exactly  the  desired  pattern  are  marginal  at  best.  The 
inaccuracy  of  the  self  and  mutual  impedance  values  are 
caused  by  a  non-infinite  ground  plane.  The  inaccuracies 
associated  with  the  measurement  of  the  elements'  input 
current  also  contributes  to  the  degradation  of  the 
radiation  pattern. 

When  evaluating  radiation  patterns  of  vertical 
phased  arrays,  one  considers  the  radiating  elements  to  be 
isotropic  point  sources.  Looking  down  upon  a  single 
element  radiating  in  free  space  one  sees  an 
omnidirectional  pattern  in  the  horizontal  (azimuthal) 
plane.  The  field  intensity  in  any  direction  from  the 
element  can  be  represented  by  a  vector  R/Q  ,  where  R  is 
the  magnitude  of  the  field  in  the  direction  Q  (see  Figure 
7).  The  angles  is  defined  as  the  angle  the  vector  R 
makes  with  the  positive  X  axis  in  a  counterclockwise 
direction. 

When  two  or  more  vertical  elements  are  configured 
in  a  geometrical  pattern  or  array  and  powered  from  a 
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common  source,  the  combination  of  each  element's  radiation 


results  in  more  radiation  in  some  directions  (constructive 
interference)  and  less  in  others  (destructive 
interference).  A  radiation  lobe  results  from  constructive 
interference  while  a  null  is  caused  by  destructive 
interference.  Also  the  radiation  pattern  is  directly 
dependent  on  the  time  delay  resulting  from  the  physical 
separation  between  the  elements  of  an  array.  It  is  this 
time  delay  which  causes  the  phase  shifts  and  leads  to 
constructive/destructive  interference. 

This  time  delay  is  also  referred  to  as  a  geometric 
phase  shift  and  an  equation  defining  this  effect  can  be 
derived  using  some  geometric  principles.  Figure  8  is  a 
representative  drawing  of  two  radiating  elements  in  the  XY 
plane  and  from  it  one  can  observe  quite  clearly  the  phase 
shift  from  a  infinite  observation  point.  From  Figure  8  it 
can  V •  hown  that  the  angle  equals 
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Tan 


-1 


Yi 


Xi 


(2.1) 


where 
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(2.2) 


From  both  elements  in  Figure  8  there  extends 
dashed  lines  which  represent  the  line  of  observation  from 
infinity  to  each  of  the  elements.  Through  the  element  #2 
located  at  (X^y^)  and  the  line  of  observation  from 
element  #1  there  is  an  imaginary  plane  of  reference  which 
is  perpendicular  to  both  lines  of  observation.  This  plane 
of  reference  is  used  to  determine  the  electrical  distance 
between  it  and  element  #1.  This  distance  represents  a 
time  delay  and  indicates  that  the  EM  radiation  from 
element  #1  is  behind  or  lagging  the  EM  radiation  from 
element  #2.  The  distance  can  be  defined  as 


X  =  RCos 


where 


R  =  V  xi  +  Yi 


(2.3) 
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FIGURE  B  -  Geometric  Phase  Shift 


2.  Provided  the  necessary  phase  shift. 

3.  Provide  the  desired  current. 


4.  Provide  the  desired  power. 


: 

i 


R-^  and  R2  are  the  input  impedance  to  the  PI  network  and 
they  are  purely  resistive.  The  desired  phase  shifts 

across  the  two  networks  are  B-^  and  B2?  1^  and  I2  are  the 
driving  point  currents  at  the  base  of  each  antenna,  and 
ZL1  and  ZL2  are  the  driving  point  impedances.  The  driving 
point  impedances  for  this  two  element  array  can  be 
calculated  using  the  mutual  impedance  equations 
characteristic  of  this  array. ^  The  mutual  impedances 
equations  for  this  array  are 


V1_I1Z11  +  I2Z12 
V2=Iiz22  +  I2Z21 

With  quarter-wave  spacing  of  quarter-wavelength  verticals 

over  an  infinite  infinitely  conducting  ground  plane,  the 

self-impedances  (Z^»Z22)  and  the  mutual  impedances 
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(Z12,Z2i)  frave  the  following  complex  values, 


(3.7) 

(3.8) 


Z11  =  Z22  =  <36-5  +  321) 


FIGURE 


located  at  points  in  a  horizontal  plane  defined  in  XY 
coordinates.  The  field  intensity  F  at  a  distant  point  in 

the  direction  0  is  the  sum  of  fields  E^  generated  by  each 
element, 


F  -  I,Ei 

i=l 


(2.5) 


where  N  is  defined  to  be  the  total  number  of  radiating 
elements  in  the  array,  and  is  9iven  by 


Ei  =V(Exi)  +  (Eyi) 


(2.6) 


where 


with 


Ex  =  X^  Ai  Cos(  ©e  +  & g  ) 


Ey  =  X  A^  Sin(  0e  +  0g  ) 


0e  =  electrical  phase 
0^  =  geometric  phase 


The  amplitude  of  the  ifc  element's  current  is  and  i 
the  phase  of  the  ith  element's  input  current.  ®g  is  the 


geometric  phase  shift  caused  by  the  elements  location  m 
the  xy  plane  and  is  defined  by  equation  (2.4). 

The  magnitude  of  the  field  intensity  is  based  in 
part  on  the  location  of  the  elements  but  more  importantly 
on  the  angle  of  observation.  Therefore  the  radiation 
pattern  at  any  point  around  the  array  can  be  computed  by 
calculating  the  total  field  magnitude  F  as  a  function  of 
(the  observation  angle).  The  calculation  implied  by 
equation  2.5  is  easily  implemented  in  FORTRAN;  a  computer 
program  named  HILPAT  will  calculate  and  plot  the  radiation 
pattern  of  N  vertical  antennas  spaced  at  any  desired 
locations  in  the  XY  plane  with  arbitrary  excitation 
currents  is  given  in  Appendix  A.  This  program  assumes  the 
antenna  placed  at  the  origin  (0,0)  in  the  XY  plane  is  the 
phase  reference  (9el  equals  0  degrees). 

By  definition,  a  null  in  a  radiation  pattern  is  a 
direction  in  which  the  field  intensity  equals  zero.11 
Thus  a  null  in  the  direction  0n  occurs  when  the  net  field 
intensity  F  satisfies 


J  Ej  =  0 
1=1 

0=0. 


(2.8) 


Since  each  is  a  complex  quantity,  the  above 


equation  implies  that 


Xx  Ai  Cos  (  9e  +  0g  )  =  0 


(2.9) 


Xx  Ai  Sin(  0e  +  0g  )  =  0 


(2.10) 


Therefore  a  single  null  in  direction  0n  can  be  produced  if 
equations  2.9  and  2.10  are  satisfied. 12  The  2N  unknowns 

(Ai»  i=l  to  N  and  ©ei#  i  =  1  to  N)  in  this  nonlinear 
transcendental  set  of  equations  can  make  solving  them 
quite  difficult. 

Solution  of  equations  (2.9,  2.10)  can  be  greatly 
simplified  if  the  phases  ©ei  are  known,  for  then  (2.9)  and 
(2.10)  become  two  linear  equations  in  the  variables  A^ . 

For  an  N  antenna  system  with  a  given  set  of  phases  0 


N/2  null  directions  can  be  produced  in  the  radiation 


attern  by  solving  N  linear  equations  for  the  current 


amplitudes  A^ . x 0 

Since  one  is  also  interested  in  the  specific 
direction  of  the  "main  lobe"  of  the  pattern  (particularly 
when  the  antenna  is  being  used  for  transmission),  it  also 
will  be  necessary  to  find  a  set  of  phases  that  will 
produce  as  much  radiation  as  possible  in  the  desired 
direction.  The  current  amplitudes  Ai  are  constrained  by 
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the  linear  equations  describing  the  null  direction  and  can 
not  be  directly  adjusted  to  provide  the  desired  "main 
lobe"  direction.14  However,  the  A^ 1 s  depend  on  the  0ei's, 
which  allows  some  measure  of  indirect  control  over  the 
main  lobe. 

The  method  for  maximizing  the  radiation  pattern  in 
a  specific  direction  0m  uses  a  direct  search  optimization 
routine.  This  optimization  routine  was  used  to  adjust 
phases  0e^  to  whatever  values  would  provide  maximum 
radiation  in  0m  under  the  constraints  imposed  by  the 

specified  0n.  Starting  values  for  the  phases  0ei  are 
input  by  the  user  and  after  a  specified  number  of 
iteration  the  phases  that  produce  maximum  radiation,  in 
conjunction  with  the  current  amplitudes  defined  by  0 
will  be  output  and  may  be  used  in  the  plotting  of  the 
resulting  radiation  pattern. 

The  described  beam  steering  procedure  was  adapted 
to  the  three  element  equilaterally  spaced  triangular  array 
in  Figure  1.  A  computer  program  (NULL)  was  developed  to 
solve  equations  (2.9)  and  (2.10)  for  the  A^  (Appendix  A) 
and  a  direct  search  optimization  routine  (MAX_HC)  was 
utilized  in  the  constrained  amplitude,  phase  maximization 
procedure  described  in  the  previous  paragraph.  Since 
there  are  only  three  (N  =  3)  elements  i*  the  equilateral 
triangular  array,  only  one  null  direction  can  be  specified 
along  with  a  direction  of  maximum  radiation.  This 
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direction  of  maximum  radiation  is  of  interest  to  a  user 
for  both  transmitting  and  receiving  purposes.  Knowing  the 
direction  of  maximum  signal  strength  allows  one  to 
determine  the  direction  of  the  signal  being  sent  out  or 
received. 

The  program  developed  for  determining  radiation 
patterns  for  the  3  element  array  of  Figure  1  was  run  with 

varying  Qm  and  0n  inputs  and  some  resulting  patterns  are 
shown  in  Figures  12  and  13.  The  main  program  MAX_RAD  uses 
the  subroutines  NULL  and  MAX_HC  to  generate  those 
resulting  radiation  patterns  an  are  listed  in  Appendix  A. 
The  specified  null  directions  must  be  kept  a  reasonable 
number  of  degrees  away  from  the  desired  main  lobe  in  order 
for  acceptable  patterns  to  result.  Null  and  main  lobe 
seperation  angles  of  greater  than  25  degrees  provide 
patterns  having  the  best  main  lobe  maximization  results, 
with  the  specified  null  always  being  present  (because  the 
null  locations  are  treated  as  constraints  that  the  direct 
search  optimization  procedure  has  no  choice  in  adjusting). 

The  current  amplitudes  and  phases  resulting  from 
this  beam  steering  method  are  used  as  input  for  both  of 
the  design  procedures  developed  in  this  thesis.  In  the 
design  by  linearization  procedure,  the  current  amplitudes 
and  phases  for  a  particular  pattern  are  used  as  the  given 
current  distributions  for  the  array  elements  with  the 
solution  of  the  nonlinear  equations  resulting  in  values 
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for  the  unknown  PI  network  parameters.  The  current 
distributions  are  used  in  the  constrained  “analysis" 
design  procedure  of  Chapter  Seven  as  the  given  (desired) 
solution  to  the  linear  version  of  the  equations  describing 
the  test  array.  Since  the  PI  network  parameters  are  the 
controlling  components  in  this  design  method,  they  will  be 
adjusted  by  the  direct  search  optimization  routine  of 
Chapter  Six  until  the  desired  current  distributions  are 
realized. 

The  beam  steering  procedure  or  the  means  of 
determining  the  current  distribution  for  a  desired 
radiation  pattern  is  of  primary  concern  during  the  early 
phases  of  designing  a  phased  array  antenna  system. 

Networks  that  will  in  fact  produce  the  desired  pattern  are 
in  themselves  a  separate  topic  the  designer  must  consider 
once  the  current  distributions  for  the  array  elements  are 
known.  An  approach  to  designing  the  phasing  networks  is 
the  primary  concern  of  this  thesis.  The  beam  steering 
procedure  described  above  was  simply  adapted  to  the 
author's  three  element  triangular  array  in  Figure  1. 

This  chapter  has  provided  the  reader  with  a 
brief  introduction  into  the  theory  behind  the  radiation 
patterns  generated  in  phased  array  antenna  systems.  The 
following  topics  were  discussed:  the  four  principle 
determinants  of  a  phased  array's  radiation  pattern,  the 
derivation  of  the  equation  describing  the  geometric  phase 


shift,  and  the  use  of  the  direct  search  optimization 
routine  to  maximize  the  radiation  in  a  specified  direction 
(  0m)  while  maintaining  the  constrained  direction  of 
minimum  radiation  (  0n).  The  programs  used  to  implement 
the  calculations  and  plotting  of  radiation  patterns  for 
the  equilateral  array  were  also  discussed. 
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III.  PEED  NETWORKS  FOR  FIXED  PATTERN 
PHASED  ARRAYS:  PRIOR  ART 

When  designing  an  antenna  system,  in  addition  to 
calculating  the  current  distributions  corresponding  to 
the  desired  radiation  pattern,  one  must  construct  a 
network  or  system  of  networks  that  supplies  these  desired 
current  amplitudes  and  phases  to  each  element. 
Traditionally  this  is  done  using  impedance  matching 
networks  connected  to  the  base  of  each  antenna.  These 
networks  produce  a  desired  phase  lag  or  lead  to  each 
particular  antenna  while  also  matching  the  driving  point 
impedance  of  the  antenna  to  the  characteristic  impedance 
of  the  transmission  line. 

There  are  many  different  methods  one  could  use  to 
match  the  antenna  driving  point  impedance  to  the 
transmission  line  supplying  the  current.  The  decision  to 
use  a  particular  method  depends  on  many  factors.  Some  of 
the  factors  that  need  consideration  are:  the  element 
operating  currents  and  voltages,  the  frequency  of 
operation,  the  degree  of  initial  impedance  mismatch,  the 
equipment  available  for  measurement  and  construction  of 
the  system,  the  required  bandwidth,  the  available  area  for 
operation  and  construction,  and  economics.  The  most 
commonly  used  methods  are:1^ 


36 


1.  Design  of  load  and  feeder  to  have  equal 


impedance. 

2.  Use  of  coupling  reactance  networks. 

3.  Use  of  tapped  transmission  lines. 

4.  Using  a  series  section  of  transmission  line  as  an 
impedance  transformer. 

5.  The  use  of  a  stub  section  of  transmission  line  as 
a  reactance  in  parallel  with  the  power  source 

at  a  point  where  the  impedance  at  that  particular 
point  will  be  equal  to  the  lines  charactersitic 
impedance. 

6.  The  use  of  a  reactance  component  in  place  of  a 
stub  line  and  electrically  equivalent  to  it. 

7.  The  use  of  a  coupled  section  of  line  in  parallel 
with  the  power  source  with  the  length  needed  to 
reflect  the  correct  amount  of  reactance  into  the 
main  power  source  at  the  point  where  an  impedance 
match  would  occur. 

8.  The  use  of  a  tapered  transmission  line  as  a 
broadband  impedance  matching  transformer. 

A  vertical  antenna  of  fixed  frequency  is  used  most 
often  at  broadcast  radio  stations  and  is  usually  matched 
using  a  coupling  network  of  reactances  in  series  with  the 
coaxial  transmission  line  from  the  transmitter.  Also,  an 
amateur  radio  operator  interested  in  only  one  frequency  of 
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operation  would  probably  find  this  method  of  impedance 
matching  most  practical. 


A.  Impedance  Matching  Networks 


With  this  in  mind,  an  example  of  an  impedance 
matching  network  designed  to  match  the  50-ohm  coaxial 
transmission  line  from  a  transmitter  to  an  antenna  having 
a  purely  resistive  impedance  of  100  +  jO  ohms  is  shown 
below.  A  phase  difference  of  +30  degrees,  between  the 
transmitter  current  and  the  current  flowing  through  the 
element  is  incorporated  into  this  design.  Therefore,  the 
phase  at  the  output  lags  the  phase  of  the  input  by  30 
degrees.  A  PI  network,  as  shown  in  Figure  14,  will  be 
used  as  an  example. 

Since  we  are  assuming  the  network  is  ideally 
reactive  with  zero  losses  the  components  of  the  network 

ZA«  ZB,  and  Zc  can  be  directly  computed  from  the  formulas 
below:  * 


Za  =  j 


Zr  =  j 


Rx  R2  Sin  B 
R2  Cos  B  —  R-^  R2 
rx  r2  Sin  B 


R1  Cos  B  -  R1  R2 


=  j  R^  R2  Sin  B 


(3.1) 


(3.2) 


(3.3) 


ll 
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B  is  the  phase  difference  between  the  input  and  output 
terminals  of  the  PI  networks,  while  and  R2  are  the 

image  impedances  with  R^  =  R^n  and  R2  =  Rj/  The  driving 
point  impedance  is  a  more  specific  term  for  the  image 
impedance  and  it  is  defined  as  the  impedance  at  the  base 
of  a  radiating  element  (vertical  antenna).  For  our 
example, 

100*50  Sin( 30 ) 

z  -  j - -  =  jl57 . 31 

50*Cos ( 30 )  -  100*50 

100*50  Sin( 30.0) 

Z  *  j -  =  -  j9i  .21 

50*Cos( 30.0)  -  100*50 

Zc  *  j  100*50  Sin(30.0)  =  j35.36 

and  the  resulting  network  is  shown  in  Figure  15. 

If  a  TEE  network  as  shown  in  Figure  16  is  used  in 
place  of  the  PI  network,  the  following  equations  would 
apply: 17 


CosB  R^  R2 
SinB 

R^  CoSi3  —  R]_  R2 
SinB 


(3.4) 


(3.5) 


Z3  =  -J 


R1  R2 
SinB 


(3.6) 


H 


wvwvwu av.v.w/.' ~wwv  ■u».wiA"J‘_v‘v  .V-W, wwlwvi 'J'.'ucwaw  >  >"‘>'y,’,>‘ v  v rj- 


and  the  resulting  reactances  being 

Zl  =  j54. 82 

Z2  =  - j31 . 78 

Z3  =  - jl41 . 42 

with  the  resulting  network  shown  in  Figure  17. 

Usually  in  the  design  of  PI  or  TEE  matching 
networks  the  desired  input  and  output  impedances  as  well 
as  the  phase  shift  are  known  and  the  network  components 
are  the  unknowns. 

Another  design  example  will  be  presented  with  the 
desired  input  and  driving  impedance  being  given  as  ZL  =  75 
“  j30  and  Z^n  =  600  +  jl50.  The  phase  shift  of  the 
network  will  be  the  difference  between  the  phase  of  the 
input  current  of  the  network  and  the  phase  of  the  base 
current  of  the  antenna  of  a  positive  60  degrees  (phase 
lead).  Since  the  B  in  equations  (3. 1-3. 6)  represents 
phase  lag,  a  negative  B  will  be  used  in  the  calculations 
in  order  to  keep  the  signs  of  the  network  components 
correct.  For  the  PI  network  of  figure  14  the  following 
values  were  calculated  using  equations  3.1,  3.2,  and  3.3 
after  transforming  the  complex  series  input  and  load 
impedances  to  their  parallel  equivalent  circiuts. 


ZA  =  +jl55 


(a  155-ohm  inductor) 
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TCrk/'X^V.  V K/’A.JVL  V1. ^yT.^-vy^yv ^ 


K.V.V.V.V.T 


ZB  =  -jH3 
Zc  =  - j!08 


(a  133-ohm  capacitor) 
(a  108-ohm  capacitor) 


and  the  resulting  network  is  shown  in-  Figure  18 
(remebering  that  the  required  network  is  lossless). 

Application  of  the  PI  network  to  a  specific  phased 
array  antenna  system  will  show  how  these  networks  are  used 
by  an  array  designer  to  provide  efficient  impedance 
matching  with  the  desired  phase  shift  in  all  transmission 
lines  leading  from  the  transmitter  to  the  array  elements. 


B.  Conventional  Phased  Array  Design 


A  two  element  vertical  array  with  impedance 
matching  networks  is  shown  in  Figure  19.  The  following 
design  approach  is  most  often  used  by  radio  broadcast 
stations  and  amateur  radio  operators  in  the  LF  to  HF 
frequency  ranges  using  vertical  antennas  positioned  above 
the  surface  of  the  earth  (non-infinite  ground  plane).  The 
transmitter  will  supply  power  to  both  elements  and 
therefore  a  PI  network  will  be  placed  in  series  with  the 
transmission  lines  leading  from  the  transmitter  to  the 


antennas . 


The  PI  networks  must  accomplish  four  things: 


1.  Provide  an  impedance  match  between  transmitter 
and  the  antenna. 


MRS 


I 


The  distance  X  is  expressed  in  wavelengths  and  can  be 
converted  to  radians  by  multiplying  by  27T/^. 

Classical  examples  of  basic  radiation  patterns 
producible  by  such  simple  arrays  are  the  endfire  and 
broadside  configurations.  The  broadside  configuration 
consists  of  two  vertical  antennas  fed  in  phase  yielding  a 
radiation  pattern  as  shown  in  Figure  9. 7  The  endfire 
array  also  has  two  verticals  (with  the  same  spacing),  but 
they  are  fed  180.0  degrees  out  of  phase.  The  radiation 
pattern  produced  by  this  arrangement  is  shown  in  Figure 


The  endfire  array  can  be  varied  in  spacing  and 
phase  in  such  a  way  as  to  produce  unidirectional  patterns. 
When  fed  90  degrees  out  of  phase  the  direction  of  maximum 
radiation  is  always  in  line  with  the  two  verticals  towards 
the  direction  of  the  antenna  receiving  the  lagging 
excitation.  A  classical  example  of  this  phenomena  is  the 
cardiod  pattern  shown  in  Fiqure  11.  This  pattern  was 
produced  with  two  quarter  wave  verticals  spaced  a  quarter 

g 

wave  apart  and  excited  90  degrees  out  of  phase. 

We  now  proceed  to  analyze  a  vertical  phased  array 
consisting  of  N  antennas.  The  vertical  elements  are 
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z12  =  Z21  =  *20-4  "  j14*2) 


The  driving  point  impedances  of  each  antenna  are 


JL1  ■ 


Vl_ 

II 

V. 


JL2 


(3.9) 


(3.10) 


and  the  power  delivered  to  each  of  the  elements  can  be 
defined  in  terms  of  the  driving  point  currents  or  voltages 
using 


pi  •  Real(Vi,Ii  ) 


(3.11) 


or 


Pi  -  Ji  RLi 


(3.12) 


where  (i)  defines  the  particular  element  in  question  and 

RLi  is  the  resistive  component  of  the  ifc^  element  liiving 

21 

point  impedance. 

Assuming  each  PI  network  will  consist  of 

tVi 

purely  reactive  components,  the  input  power  to  the  i 
network  equals  the  output  power  (P^)  of  each  network.  The 


*5 


total  input  power  to  the  array  can  then  be  expressed  as 


PT  -  X,  »i 
1  =  1 


(3.13) 


S^nce  one  would  want  the  input  impedance  to  the  array  to 
match  the  characteristic  equation  of  the  transmission 
line,  we  have 


R.  —7, 
in  ^o 


(3.14) 


and  by  using  50-ohm  coaxial  line  the  the  expression  for 
the  input  power  becomes 


PT  “ 


(3.15) 


with  Vin  being  the  voltage  from  the  transmitter.^ 

The  input  power  of  the  ith  network  can  be  used  to 
find  the  input  impedances  for  each  network  by  defining  the 
phase  of  the  input  voltage  (V^n)  as  0.0  degrees.  Equation 
3.15  then  takes  the  following  form 


/  N 

i=iPi 


(3.16) 


I 


Pi  =  Real  vin 


(3.17) 


With  the  phase  of  Vin  equal  to  0.0  degrees,  the  input 
resistance  to  each  networks  must  be 


4  =  lviri 


(3.18) 


The  PI  networks  must  be  able  to  match  any  complex 
driving  point  impedance  at  the  base  of  each  antenna 

with  any  input  resistance  .  By  using  the  PI  network  in 

this  example,  one  can  utilize  equations  (3. 1-3. 3)  to 
calculate  the  components  of  the  networks  once  the  driving 
point  and  input  impedance  that  result  from  a  given  current 
excitation  have  been  calculated  using  equations  (3.12)  and 
(3.18).  The  network  must  also  produce  the  desired  phase 
shift . 


C.  Feed  Network  Design  Example 


The  design  method  discussed  will  be  used  in  the 


$1 

I 


52 


following  specific  example.  In  Figure  20  there  are  two 
quarter-wave  vertical  antennas  spaced  a  quarter  wave  apart 
above  an  assumed  infinite  ground  plane.  The  base  currents 
I  and  I2  are  given  as  desired  values  and  in  phasor  form 
they  equal 


:i=  iA  =  i 


i2  =  izko  =  ji 


With  these  currents  the  driving  point  impedances,  using 
the  mutual  impedance  equations  3.7  and  3.8  and  equations 
3.9  and  3.10,  are  found  in  the  following  manner; 


ZL1  =  Z11  +  z12 


=  36+j21+(20.4-jl4.18) (jl) 

=  (36  +  14.2)  +  j ( 21 . 0  +  20.4) 


ZL1  =  50.2  + j41 . 4 


ZL2  =  z22 


(3.19) 


(3.20) 


=(36  +j21)(-jl)+(20.4-jl4.18) 

=  (21.0  +  20.4)  -  j ( 36  +14.18) 

ZL2  =  41,4  "  j50 . 2 

The  power  for  the  each  antenna  is  then  solved  for  using 


f 

$ 

r§ 


s 

S'' 

g: 

& 


Element  #2 


equations  3.13  and  3.14  in  the  following  manner; 

Px  =  |ljj2  RL1=l*Real[ZL1]=l*50.2=50.2  watts  (3.21) 

2 

P2  =  l1^  RL2=1*Real£ZL2^=1*41*4=41'4  watts  (3.22) 

Once  the  power  to  each  network  is  known  the  input 
resistance  to  each  antenna  can  be  found  using  equation 
3.18  as  shown  below. 

2 

RX  =  I  Vin| 
pl 

where 

Vin  =  y  50*  ( PL  +  p2)  =^50*91.6  =  67.5  (3.23) 

So 

=  91.235  ohms  R2  =  110.628  ohms 

The  phase  shift  between  the  input  and  output  port 
of  the  PI  networks  can  be  calculated  by  converting  the 
complex  values  for  ZL1  and  into  phasor  form.  The 

driving  point  impedances  are  now 


ZL1  =V(51-2>  + 


51.2 


2 

41.4)  +  (51.2) 


51.2 


41.4 


with 


ZL1  =65.8  /39.5  ZL2  =  65.8  /50.4 

From  the  previous  calculations  the  values  for  Rj/ 

R2'  zL1'  zL2  are  now  known.  The  load  impedances  are 
converted  into  their  equivalent  parallel  circuits  and  the 
parallel  resistance  components  are  be  plugged  into 
equations  (3.1  -  3.3)  for  R2  in  order  to  solve  for  the 
reactive  components  of  the  PI  networks.  The  parallel 
reactive  components  are  coupled  with  the  parallel  branch 
of  the  PI  networks  and  the  resulting  reactance  becomes  the 
value  of  the  particular  parallel  branch  on  the  network. 
Substitution  of  R-^  and  ZL1  for  R]^  and  R2  of  equations  3.1, 
3.2,  and  3.3  results  in 

(91.2) (124.01)Sin(39.5) 

ZA1  =  3 


(124.01) (Cos(39. 5) )  -  ( 124 . 01 ) ( 91 . 2 ) 


(91.2) (124.01)Sin(39.5) 


(91.2) (Cos(39.5) )  -  (124.01)(91.2) 


ZC1  =  j  (124.01)01.2)  Sin(  39 . 5  ) 


upon  completion  of  the  computations  above 


ZA1  =  - j674. 98 


ZB1  *  -j302.46 
ZC1  =  367.65 


The  substitution  of  R2  and  Z^  of  the  other  PI  network  for 
Rl  and  R2  in  equations  3.1,  3.2,  and  3.3  take  the 
following  form 


ZA2  "  3 


JB2  =  3 


(110.6) (69.56)Sin(50.4) 

( 69 . 56 ) ( Cos (50.4))  -  ( 69 . 56 ) ( 110 . 6 ) 


(110.6)(69.56)Sin(50.4) 

( 110. 6 ) ( Cos (50.4))  -  (69.56) (110.6) 


ZC2  =  3  (69.56) (110.6)Sin(50.4) 


The  above  computations  result  in 


Z.,  =  ~ jl36 . 67 


ZB2  "  -3259.56 

ZC2  =  j67 . 58 


The  inductance  and  capacitance  values 
corresponding  to  calculated  reactance  values  for  the 
paramaters  of  the  two  PI  networks  can  be  found  assuming  a 
frequency  of  7  MHz  with  the  following  equations: 


C  = 


wXc 


where 


XL  =  a  positive  j(  ZA  or  ZB  or  ZC) 
Xc  =  a  negative  j(  ZA  or  ZB  or  ZC) 
w  =  27Tf  (  f  =  frequency) 


Therefore,  the  PI  network  parameters  for  the  two  networks 
are 


c 

ZA1 

=  33 . 68pF 

CZA2 

=  166 

.  7pF 

CZB1 

=  75 . 64pF 

CZB2 

=  87 

.  6pF 

LZC1 

=  1 . 53uH 

LZC2 

=  1. 

53uH 

and  the  networks  are  shown  in  Figure  21. 


FIGURE  21  -  PI  Networks  of  2  Element  Antenna  System 


D.  Conclusion 


This  chapter  listed  many  of  the  currently 
acceptable  methods  used  for  impedance  matching  in  phased 
array  antenna  systems.  Discussion  of  PI  and  TEE  matching 
networks  was  presented  in  some  detail  for  two  reasons:  PI 
and  TEE  networks  are  an  often  used  form  of  impedance 
matching,  and  the  PI  network  will  be  used  in  the 
equilateral  array  studied  in  this  thesis  as  the  circuit 
for  controlling  the  radiation  pattern  of  the  array.  A 
simple  two  element  array  with  series  PI  networks  was  used 
to  correctly  match  the  antennas  feed  point  impedances  with 
the  impedances  being  seen  from  the  transmitter  side  of  the 
network  looking  towards  the  individual  element. 


IV.  ANALYSIS  OF  A  THREE  ELEMENT  PHASED  ARRAY 


A  three  element  equilateral  triangular  phased 
array  antenna  system  was  chosen  to  test  the  two  design 
methods  developed  in  Chapters  Five  and  Seven  because  of 
the  comparatively  small  number  of  radiating  elements 
involved.  In  this  chapter  the  analysis  of  the  three 
element  triangular  array  will  be  presented. 

The  "analysis"  procedure  to  be  demonstrated  in 
this  chapter,  as  well  as  the  "design"  procedures  of 
Chapters  Five  and  Seven,  can  be  extended  to  larger  arrays. 
Figure  22  shows  an  NxN  element  array  typical  of  larger 
designs. 

The  three  elements  of  the  equilateral  array  used 
throughout  this  thesis  are  quarter-wavelength  vertical 
antennas  spaced  a  quarter-wavelength  apart.  The  coaxial 
transmission  line  connected  between  the  pairs  of  elements 
and  in  series  with  the  PI  networks  has  a  characteristic 
impedance  of  50  ohms.  The  array  is  excited  by  a  50-ohm 
transmission  line  from  the  transmitter  to  element  #1. 


Parameters  of  the  PI  networks  are  assumed  to  be  purely 
reactive  in  all  calculations. 


Transmission  Line  Equations 
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Coaxial  transmission  lines  are  inherently  lossy. 
The  physical  characteristics  of  the  lossy  transmission 
lines  can  be  described  as  follows :2^ 


VrCosh y 

1 ZYL 

+ 

IrZQSinhy 

'ZYL 

(4.1) 

XrCosh  y 

/  ZYL 

+ 

V 

— —  Sinh- 
Zo 

VZYL 

(4.2) 

Z  =  R  +  jwL  (series  impedance  per  unit  length) 
Y  =  G  +  jwC  (shunt  admittance  per  unit  length) 


Vr  =  receiving-end  voltage  (load) 

Vg  =  sending-end  voltage  (generator) 

Ir  =  receiving-end  current  (load) 

Ig  =  sending-end  current  (generator) 

R  =  reactance  L  =  inductance 

C  =  cpacitance  G  =  leakage 


However,  for  some  communications  applications  in  the 
amateur  amd  broadcast  radio  fields,  the  transmission  lines 
can  be  assumed  to  be  lossless  because  the  actual  losses  in 
the  lines  are  quite  small.  For  lossless  transmission 
lines,  equations  (4.1  and  4.2)  take  on  the  following 
form: 25 


Z 

o 


V 


s 


Characteristic 

27TL 


VrCos 


X 


Impedance 

jIrZoSin 


=  50+j0  ohms 
2TTL 

X 


(4.3) 


(4.4) 


t-V  -  ■ r,  c,  /, 


S" 


■.  •.  .  *.  ■  \  V  ", 

•  v-V-  ----- 


>  .V  .V  . 


■V" 


2TTl 

T 
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Ig  =  IrCOS 


27TL 


+  j - Sin 

Z  o 


(4.5) 


Equations  of  the  above  form  will  describe  each  segment  of 
transmission  line  used  in  the  interconnection  networks. 


B.  Mutual  Impedance  Equations 


An  antenna  with  current  flowing  in  it  will  induce 
a  voltage  in  any  other  antennas  in  the  nearby  area.  The 
antennas  will  behave  as  though  they  were  coupled  to  each 
other  and  the  induced  voltage  in  a  second  antenna  divided 
by  the  current  flowing  in  the  excited  antenna  represents  a 
mutual  impedance.  Denoting  the  mutual  impedance  between 
elements  1  and  2  by  Z±2  an  alternative  form  of  equation 
3.7  can  be  written  as 


Ilzll 


‘12  = 


(4.6) 


where  E ^  is  the  induced  voltage  in  the  second  antenna  and 
is  the  excitation  current  in  the  first  antenna  (see 
Figure  21).  The  mutual  impedance  is  expressed  in  terms  of 
the  current  at  the  base  of  the  first  antenna  and  the 
voltage  (induced  by  the  current  in  the  first  antenna)  at 
the  base  of  the  second  antenna.  This  mutual  impedance 
between  two  elements  also  varies  with  the  geometry  of  the 


r»t, 

§ 

4 


TOW 


antenna  and  geometry  of  the  array.  The  voltage  and 
current  relationships  in  arrays  of  elements  can  be 
represented  by  summing  the  voltage  generated  in  the  nth 
antenna  caused  by  currents  in  all  n  antennas  (extension  of 
equations  3.7  and  3.8),  as  follows: 

V1  =  Ilzll  +  I2Z12 

v2  =  Ilz12  +  I2Z22 
•  •  * 

•  •  • 

Vn  *  ^ln  +  I2z2n 

where 

Vn  =  voltage  applied  to  the  base  of  element  n. 

In  =  base  current  flowing  in  antenna  n. 

Znn  =  self  impedance  of  element  n. 

Z—  =  mutual  impedance  between  antennas  i  and  j. 

The  mutual  impedance  (Z12)  between  two  antennas  is 
defined  as  the  ratio  of  the  voltage  at  the  base  of  the 
second  antenna  to  the  current  flowing  at  the  base  of  the 
first  antenna.  The  self  and  mutual  impedance  values 
representing  the  elements  of  the  equilateral  triangular 
array  are  calculated  from  theory,  using  an  assumed  current 
excitation  on  elements  of  the  same  geometry  and  spacing  as 


Inzln 

Inz2n 


Inznn 


(4.7) 

(4.8) 


(4.9) 


Mutual  impedance  measurements  are  difficult, 


the  array.26 

as  they  are  easily  affected  by  the  surroundings  in  the 
area  of  the  array  and  the  conductivity  of  the  earth's 
surface.  Therefore,  these  values  should  not  be  assumed 
accurate  and  measurements  should  be  made  at  the  actual 
constructed  array  with  the  measured  values  being  used  in 
the  re-evaluation  of  any  calculations. 


C.  PI  Network  Equations 

The  PI  network  shown  in  Figure  4  is  the  form  to  be 
used  in  the  interconnection  medium  of  the  phased  array 
antenna  system  of  Figure  5.  The  PI  network  is  connected 
in  series  with  two  short  sections  of  transmission  line 
between  each  pair  of  elements  in  the  array.  Kirchhoff's 
Voltage  Law  is  used  in  the  analysis  of  the  PI  network 
where 


1  _  1  "j 

' '  r  _  — — —  f 

sC  jwC  wC 


1 

with  Xr  =  - 

wC 


sL  =  jwL  =  jXL  ,  with  XL  =  wL 

L  =  inductance  C  =  capacitance 

XL  =  reactive  inductane  in  ohms 
Xc  =  reactive  capacitance  in  ohms 


The  KCL  equations  at  nodes  #1  and  #2  then  can  be  written 


(1/ jxc31) 


<V3A  -  V1A} 


3XL13 


After  transferring  X^^  and  Xq^i  to  the  numerator  of  the 
second  component,  equations  (4.10)  and  (4.11)  become 


X13A  -  V1A(3XC13>  + 


r31A  “  V3A^XC31}  + 


(V1A  "  V3A} 


<V3A  -  V1A> 


3xL13 


(4.12) 


(4.13) 


The  currents  (I^3A  and  I31A)  and  voltages  (V1A  and 
V3A^  (^.12)  and  (4.13)  are  complex  quantities,  but  the 
impedance  representations  for  the  inductor  and  capacitors 

are  only  shown  with  the  imaginary  components  (“jxLi3' 
^XC13'  jxC3l^  because  the  real  components  are  assumed  to 
be  zero;  that  is,  only  lossless  networks  are  being 
considered  here. 

Combining  the  transmission  line  equations,  the 


a 


effects  of  mutual  impedance  between  the  elements,  the 
Kirchhoff's  Current  Law  equation  at  the  base  of  each 
antenna,  and  the  PI  network  equations  yields  a  set  of  24 
complex  equations  describing  the  entire  array  as  follows: 


Mutual  Impedance  equations: 


V1  =  Ilzll  +  I2Z12  +  I3Z13 
v2  =  Ilz21  +  I2Z22  +  I3Z23 
V3  =  Ilz31  +  I2Z32  +  I3Z33 


Transmission  Line  equations: 


p=( 27TL) *3/2,  L=l/8,  Zo=50.0 


V1  =  VlACos(p)  +  jZ0I13ASin(p) 


r13  =  Ii3ACos(p)  +  —  V1ASin(p) 

z 

o 


v3  =  V3ACos(p)  +  jZ0I31ASin(p) 


I31  =  I3lACos^P^  +  ~T~  V3ASin(p) 

Zo 


V1  =  V1BCos(p)  +  jZ0I12BSin(p) 


X12  =  I12BCos(P)  +  —  V1BSin(p) 

Zo 


(4.14) 

(4.15) 


(4.16) 


(4.17) 


(4.18) 


(4.19) 


(4.20) 


(4.21) 


(4-22) 
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IO 

II 

V2BCos(p)  + 

^ZoI21BSin(P) 

(4.23) 

j 

*21 

=  ^21B('os 

+  -  V2BSin(p) 

zo 

(4.24) 

II 

CM 

> 

V2cCos(p)  + 

^ZoI23CSin(P) 

(4.25) 

j 

I23 

=  I23cCos(p) 

+  -  V2cSin(p) 

zo 

(4.26) 

V3  = 

V3cCos(p)  + 

^ZoI32CSin(P) 

(4.27) 

j 

I32 

=  l32CCos(P) 

+  -  V3cSin(p) 

zo 

(4.28) 

Kirchhoff's  Current  Law  (at  base  of  each  antenna): 


X1 

+ 

II2 

+ 

Z13  =  ZT 

(4.29) 

X2 

+ 

*21 

+ 

Z23  =  0 

(4.30) 

Z3 

+ 

*31 

+ 

z32  =  0 

(4.31) 

PI  Network  equations: 


I13A  "  VIA(^XC13)  + 

I31A  =  V3A(jXC31)  + 


^V1A-V3A) 

^XL13 

(V3A"V1A> 


jX. 


LI  3 


(4.32) 


(4.33) 


V1B<JXC12> 


+ 


*  v1B~v2B/ 
3XL12 


(4.34) 


I12B  " 


I21B  =  V2B<  3XC21 ^  + 


^V2B~V1B) 

^XL12 


I23C  =  v2C^XC23^  + 


( v2C“V3C^ 
3xL23 


I32C  ~  v3c( 3Xc32^  + 


(V3C“V2C) 


jX. 


L23 


(4.35) 


(4.36) 


(4.37) 


In  order  to  simplify  the  implementation  of  these 
equations  into  FORTRAN  code,  the  transcendental  components 
of  the  transmission  line  equations  are  (from  this  point 
on)  represented  in  the  following  form: 


2  7TL 

K,  =  COS  — 


27TL  27TL 

K2  =  ZQSin  — —  K3  =  Sin  - - 

A  z0A 


The  variables  in  this  se*-  of  24  complex  equations 
are  labeled  on  the  equilateral  phased  array  diagram  in 
Figure  5.  The  variables  V1 ,  v2 ,  V3  are  the  voltage  at  the 
base  of  each  of  the  three  elements  of  the  array.  1^,  I2' 
I3  are  defined  as  the  current  distributions  (current 


c 


A 


flowing  through  each  element).  The  variables  V^A,  V3 A» 

V1B'  V2B'  V2C'  V3C  are  the  v°lta9es  at  the  ends  of  the 
transmission  line  sections  that  are  not  connected  to  the 

base  of  the  antenna,  but  to  the  PI  networks.  The  I12* 

I21'  I13.  I31*  I23'  I32  variables  are  the  currents  flowing 
from  the  array  elements  through  each  section  of 
transmission  line  towards  the  PI  networks.  The  currents 
flowing  at  the  end  of  the  transmission  line  sections 
connected  to  the  PI  networks  are  the  variables  I13A,  I31A' 

I12B»  I21B'  I23C'  I32C*  The  iast  variable  of  this  set  of 
equations  is  IT*  which  is  defined  as  the  input  current 
from  the  source  of  power  for  the  array  and  is  input  at 
element  #1  of  the  array. 

In  the  "analysis"  of  the  equilateral  array  the 
reactive  components  of  the  PI  networks  are  given  as  known 
values  in  the  set  of  equations  describing  the  system.  The 
currents  flowing  in  the  elements  of  the  array  become  the 
unknown  variables  to  be  calculated.  Upon  close 
examination  of  this  linear  set  of  equations  the 
possibility  of  substituting  one  of  the  groups  of  equations 
into  one  of  the  other  groups  of  equations  is  realized. 
Substitution  of  the  mutual  impedance  and  PI  network 
equations  into  the  transmission  line  equations  was 
performed  to  bring  about  a  reduction  in  the  number  of 
equations  as  well  as  in  the  complexity  of  the  set  of 
equations  needing  to  be  solved.  Upon  completion  of  the 


above  mentioned  simplification  measures  for  the  "analysis" 
situation,  the  24  complex  equations  (4.14  thru  4.37)  take 
the  following  complex  linear  form: 
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J1  +  I12  +  I13  =  IT 

12  +  *21  +  123  =  0 

13  +  X31  +  I32  =  0 


I1Z11+I2Z12+I3Z13-V1A 


K, 


Kl“ff2XC13  + 


XL13 


+  V 


3A 


-K- 


"L13 


I13“V1A 


7K1 

Kl3xC13  ”  -  +  3K3 

XL13 

-K„ 


-  V 


3A 


3*1 

XL13 


=  0 


I1Z31+I2Z32+I3Z33“V1A 


^13 


+  V 


3A 


Kl  + 


K-- 


"LI  3 


~K2XC31 


I31  "  V1A 


-  V 


3  A 


"LI  3 


3K1 

Kl7xc32  ”  +7K3 

XL13 


=  0 


I1Z11+I2Z12+I3Z13-V1B 


K, 


K1-K2XC12+ 


^12 


-K. 


+  V 


2B 


"LI  2 


1 1 3 -V1 B 


3K1 

K1]XC12  “  ~  +  3K3 


LI  2 
-K, 


jxi 


-  V 


I1Z21+I2Z22+I3Z22-V1B  -  +  V: 


2B 


2B 


Kl  + 


"LI  2 

K2 

CL12 


=  0 


~K2XC2lJ 


(4.38) 

(4.39) 

(4.40) 

0  (4.41) 

(4.42) 

=0  (4.43) 

(4.44) 

0  (4.45) 

(4.46) 

=0  (4.47) 


I31  “  V1B  -  "  V2B  K1^XC21  “ 


+  1r,-a  =  0 


XL12 


(4.48) 


Ilz21‘fI2z22+I3z23-v2C  K1-K2XC23  + 


*L23 


f  jKl  jKl 

I13-V2C  k13xC23  “  -  +  3x3  -  V3C  -  -  0 


XL23 


=0  (4.49) 


(4.50) 


XL23 


I1Z31+I2Z32+I3Z33~V2C  3 -  +  V3C  Kl+  -  “K2XC32  0  ^4’51^ 


*L23 


I31  “  V2C  -  “  V3C  KlixC32  “  -  +3K3  -  0 

XL23  XL23  i 


(4.52) 


Therefore  in  the  analysis  case  the  unknown  variables  are 


„  r12 *  Z21 '  I13'  Z31 '  I23*  z32'  V1A'  V3A'  V1B'  V2B'  v2 C' 
V3C'  Ii,  I2,  13 


and  the  known  variables  are 


XC12 '  XC21 '  XC13 '  XC3 1 '  XC23 '  XC32'  IT 
XL12 '  XL13 '  XL23 

The  above  set  of  linear  (since  the  network 
parameters  are  given)  equations  can  be  solved  using  matrix 
techniques  for  solving  a  set  of  linear  simultaneous 
equations  of  the  form  Ax=b.  Here  A  is  a  15x15  complex 
matrix  consisting  of  the  coefficients  of  the  unknown 
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variables,  and  the  complex  matrix  b  consists  of  the  only 

constant  in  the  system  (the  array  input  current  IT) •  The 
PI  network  parameters  are  known  but  they  are  coefficents 
of  the  variable  voltages  at  the  ends  of  the  transmission 
line  sections  connected  to  the  PI  networks  (V1A,  V3C) 

in  the  15  linear  equations  (4.38  thru  3.52). 

A  FORTRAN  subroutine  was  written  which  solves 
equations  4.38  thru  4.52  using  a  linear  equation  solver 
routine  (LEQ2C)  located  in  the  IMSL  subroutine  library 
available  on  the  VAX  11/780.  This  analysis  subroutine 
ANAL_DESG  returns  to  the  user  the  complex  current 
distributions  of  the  array  elements  after  being  supplied 
the  array's  input  current  (IT)  and  the  complex  PI  network 
parameters.  The  source  code  for  this  subroutine  CURR  and 
the  main  program  ANAL_DESG  are  listed  Appendix  B. 


D.  Results 

An  example  run  of  this  subroutine  is  shown  below; 
it  was  supplied  with  constant  magnitudes  for  all  the  PI 
network  parameters  and  a  value  for  the  input  current  IT. 


ilALM.M 


XL23  =  100  XC23  =  -25  XC32 

I-p  (input  current)  =  1.0  +  jO 
Resulting  Magnitude  and  Phase 


=  -5 


Mn  =  1.0 

phu  =  o.o 


mI2  =  0.28  Mi3  =  0.28 

PHj2  =  95.6  PHi3  =  95*58 


The  values  for  the  PI  network  paramters  in  the 
above  example  were  choosen  arbitrarily  and  thus  do  not 
necessarily  result  in  current  distributions  that 
correspond  to  an  useful  radiation  pattern.  The  radiation 
pattern  corresponding  to  the  current  distributions  found 
above  (11,12,  &  13)  is  shown  in  Figure  23. 

Below  is  a  list  of  two  more  examples  showing  the 
given  PI  network  parameters  and  the  resulting  current 
distributions . 


XL13  =  100 

XC1 3  = 

-50 

XC31  = 

-500 

XLi2  =  100 

XC12  = 

-5.67 

XC21 

-50 

xl23  =  100 

XC23  = 

-40 

XC32  = 

-500 

IT  (input  current)  =  1.0  +  jO 


y 

«! 


„  .  .v 

»Vv 


GURE 


Resulting  Magnitude  and  Phase 


Mn 

=  1.0 

Ml  2 

=  0.28 

“13 

=  0.28 

PHI1 

=  0.0 

PHI2 

=  95.48 

PHI3 

=  95.48 

XL13  = 

10 

XC13  = 

-5 

XC31  =  ~5° 

XL12  = 

15 

XC12  ~ 

i 

o 

• 

o 

H* 

XC21  =  -1° 

XL23  = 

20 

XC23  = 

-0.02 

XC32  =  “5 

IT  (input  current)  =  100.0  +  j45 


Resulting 

Magnitude 

and  Phase 

Mn  =  1.0 

MI2 

=  0.11 

mI3 

PHI1  =  °*° 

PHI2 

=  80.67 

phI3 

E.  Summary 


This  chapter  has  demonstrated  that  the  phased 
array  antenna  system  represented  in  Figure  5  can  be 
analyzed  by  solving  the  15  linear  simultaneous  equations 


that  describe  the  system.  The  current  distribution  for 
each  element  of  the  array  are  the  three  variables  of  most 
concern  resulting  from  the  solution  to  equations  4.38  to 
4.52.  However,  these  resulting  current  distributions  are 
as  arbitrary  (useless)  as  the  corresponding  PI  network 
parameters  specified  by  the  user. 

Therefore,  in  order  to  design  a  phased  array 
antenna  system  based  on  the  interconnection  medium 
illustrated  in  Figure  5,  the  24  complex  equations  that 
describe  the  physical  and  electrical  characteristic  of  the 
system  must  be  solved  using  specified  (desired)  current 
distribution.  These  current  distributions  will  correspond 
to  a  desired  radiation  pattern.  This  idea  is  explored  in 
Chapter  Five. 

Presentation  and  analysis  of  the  equations 
describing  the  physical  and  electrical  characteristics  of 
the  three  element  equilateral  triangular  array  was  the 
main  theme  of  Chapter  Four.  The  substitutions  of  the 
mutual  impedance  and  PI  network  equations  into  the 
transmission  line  equations  was  accomplished  and  the 
resulting  set  of  15  equations  was  presented  in  their 
complex  linear  form,  with  the  known  and  unknown  variables 
clearly  indicated.  All  of  the  necessary  equations  and 
theories  has  been  discussed  and  the  stage  is  now  set  for 
developement  and  discussion  of  the  PI  network  design 
procedures  to  be  presented  in  Chapters  Five  and  Seven. 


V.  INTERCONNECTION  NETWORK  DESIGN 


BY  LINEARIZATION:  AN  EXTENSION 


OF  NEWTON'S  METHOD 


This  chapter  will  present  a  detailed  explanation 
of  a  procedure  capable  of  solving  the  equations  that 
describe  the  unique  phased  array  antenna  system  of  Figure 
5  in  the  "design"  situation.  As  mentioned  in  the 
introduction,  the  design  process  pertains  to  the  situation 
where  the  currents  on  the  radiating  elements  are 
classified  as  given  (known  constant  complex  quantities), 
while  the  PI  network  parameters  are  the  unknown  variables 
of  interest.  Therefore,  the  "design"  of  interconnection 

networks  for  the  array  is  just  the  opposite,  where  the 
variables  of  interest  are  concerned,  of  the  "analysis"  of 
interconnection  networks  that  was  presented  at  the  end  of 
Chapter  Four. 

To  a  potential  user  of  the  equilateral  array 
proposed  in  this  thesis,  the  "design"  situation  would  be 
of  most  interest.  The  user  probably  already  has 
calculated  the  current  distributions  on  the  elements  of 
the  particular  array  (discussed  in  Chapter  Two)  and  would 
need  only  the  proper  PI  network  parameters  to  realize  the 
desired  radiation  pattern.  The  design  procedure  of  this 
chapter  provides  these  PI  network  parameters  given  the 
current  distributions  on  the  elements  as  input. 


w: 


The  discussion  will  begin  with  an  explanation  of 
the  substitutions  made  among  the  24  equations  describing 
the  array  (equations  4.14  thru  4.37),  justification  for 
giving  some  variables  a  constant  value,  and  why  one  of  the 
Kirchhoff's  Current  Law  equations  was  disregarded  (not 
used)  as  an  equation  describing  the  array  in  this  design 
procedure.  The  nonlinear  equations  are  presented  in  their 
complex  form  and  a  detailed  discussion  of  the 
linearization  process  performed  on  these  equations 
follows.  The  linearized  set  of  complex  design  equations 
and  the  results  of  an  implementation  of  this  design  method 
on  the  equilateral  triangular  array,  with  Xcl2  and  Xc23 
being  the  unknown  variables,  is  shown  as  an  example. 

In  this  design  method  the  parameters  of  the  PI 
networks  are  unknown  values.  The  desired  current  for  each 
element  in  the  array  will  be  input  as  known  values  from 
the  program  MAX_RAD  (discussed  in  Chapter  Two  and  listed 
in  Appendix  A) .  Reducing  the  amount  of  equations  and  the 
number  of  unknown  variables  by  substitution  of  the  mutual 
impedance  and  PI  network  equations  into  the  transmission 
line  equations  was  done  in  order  to  simplify  the 
computational  procedure. 

Substituting  the  mutual  impedance  equations  into 
the  tansmission  line  equations  for  the  variables  ,  V2, 
and  and  the  PI  network  equations  for  the  variables 

Z13A*  I31A'  Z12B'  Z21B'  Z23C'  and  Z32C  yields  a  set  of  12 


equations.  Together  with  Kirchhoff's  Current  Law  ( KCL) 
equations,  the  number  of  equations  now  describing  the 
system  has  been  dropped  to  15.  The  voltage  at  the  base  of 


each  antenna  (V1#  V2#  V3)  are  constants  for  this  method 
since  the  input  currents  1^  i2,  and  I3  are  given  and  the 

self  and  mutual  impedance  values  Z^,  Z^2,  and  Z13  are 
also  known  values.  Thus,  equations  (4. 7-4.9)  give  t  v2, 
and  V3  directly.  In  order  to  solve  these  15  equations 
(4.38  thru  4.52)  simultaneously  there  must  be  the  same 
number  of  equations  as  there  are  unknowns.  The  current 
flowing  into  the  array  from  the  transmitter  (IT)  is  an 
unknown  variable  that  must  be  disregarded  in  this  design 
procedure.  So  equation  4.38,  which  is  the  KCL  equation  at 
the  base  of  antenna  #1,  is  taken  out  of  the  system  in 
order  to  balance  the  number  of  equations  with  the  number 
of  unknown  variables.  Thus,  equation  4.38  is  assumed 
non-existent  and  is  not  included  in  the  set  of  14 
nonlinear  design  equations  describing  the  array.  There 
remains  14  equations  with  21  unknown  variables. 

In  order  to  make  the  14  equations  solvable,  seven 
of  the  nine  PI  network  parameters  were  assigned  specified 
values  and  so  they  become  knowns,  yielding  14  nonlinear 
complex  equations  in  14  unknowns.  The  PI  network 
parameters  that  were  not  given  a  value  become  the  unknown 
variables  of  interest.  For  the  example  shown  in  this 
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C23 


chapter  X( 


and  X. 


are  the  unknown  network  paramaters. 


1 ; 


v.w'jfrjf'jiv 


A.  Design  Equations 


The  14  nonlinear  complex  equations  have  been 
arranged  in  a  form  that  is  compatible  with  the  linear 
equation  solver  routine  (LEQ2C)  of  the  IMSL  subroutine 
library.  The  design  equations  are  listed  below  in  this 
form: 


I21  +  I23  =  -I2 
I31  +  *32  ~  -I3 


(5.1) 

(5.2) 


V1A  K1  ~  K2XC13  + 


*L13 


+  V- 


=  V1  (5.3) 


T  r  jRi  jRi 

I13_V1A  K1 3XC13 - +  iK3  ~  V3A -  =  0  (5.4) 

XL13  j  XL13 


+  V3A  K1  + 


XL13 


-  k2xc31  =  v3 


'1  iKi 

-  V3A  KxjXc31  - 


131-VlA 


V1B  K1  “  k2xC12  +  -  +  v 


(5.5) 


=  0  (5.6) 


2B  -  =  Vl  (5.7) 

XL12 


jKl  1  3Ki 

I12-V1B  K1^XC12  “  “  +  3K3  “  V2B  “  0  (5*8) 


IB  -  T  V2B  *1  ^  - 

XL12  L  XL12 


^2*021 


v2 


jK, 


1 2 1 _V1 B 


*1,12 


-  V 


2B 


K1^XC21  "  +  3K3 


y,L12 


=  0  ( 


2C 


K1  ~  k2XC23  + 


K- 


eL23 


— K<- 


+  V 


3C 


=  V9  I 


CL23 


I23"V2C 


jKl  , 

Kl3xG23  "  T -  +  3K3 


V 


2C 


-k2 


lL23 


+  V 


3C 


Ki  + 


*1,23 

K2 


jKi 

-  v3c  — -  -  0  ( 


XL23 


XL23 


~  K2XC32 


=  V- 


I32“V2C 


*1,23 


-  V 


3C 


3K1 

K1^XC32  "  +  3K3 


"L23 


=  0  ( 


Since  the  unknown  PI  network  parameter  are  Xd2  and 
the  unknown  variables  of  the  design  equations  listed 
are: 


1 12 '  I21'  I23'  *32 '  *13 '  I31'  V1A'  V3A'  V1B' 
V2C'  V3C'  XC12 '  XC23 

and  the  known  variables  are 


V, 


V- 


'11' 


'12  ‘ 


'13' 


*21' 


5.10) 

[5. 11) 

:  5 . 12) 

:5.i3) 

:5.i4) 

XC23 

above 

V2B' 

!22' 


z23'  z31'  z32 '  z33 '  XL12'  XC21'  XL13'  XC13'  XC31 
XL23 '  XC32 

'The  K^,  K2#  and  K3  terms  are  defined  as  they  were 
in  Chapter  Four  on  page  69* 

Upon  observation  of  the  14  design  equations  one 
can  see  the  nonlinearity  is  due  to  the  occurrence  of 
crossproducts  between  the  unknown  phasor  voltages  and  the 

unknown  PI  network  parameters.  For  instance,  since  Xci2 
is  an  unknown  network  parameter  that  needs  to  be  found,  an 

example  of  a  nonlinear  component  is  the  VibxC12  term  in 
equations  5.7  and  5.8.  Also  the  component  v1AXrq3  is 
nonlinear  in  equations  5.4  and  5.5  if  the  XC13  network 
parameter  was  an  unknown  capacitor  of  interest. 


B.  Taylor's  Series  Expansion 


In  order  to  solve  the  14  design  equations  exactly, 
one  must  use  a  linear  equation  solver  routine  which 
utilizes  matrix  algebra  techniques  (Ax=b).  Therefore  the 
nonlinear  equations  must  be  linearized.  The  method  used 
to  do  this  is  a  generalization  of  Newton's  method  for 
finding  the  roots  of  f(x)=0.  The  generalization 


involves  extension  to  the  multivariable,  complex  variable 
case  presented  by  the  solution  of  equations  5.4  thru  5.15. 


r,»,i 


> jjsVj 


The  extension  is  as  follows. 

Transposing  all  terms  of  equations  5.4  thru  5.15 
to  the  left  of  the  equal  sign,  the  design  relations  may  be 
written  as 


fi<xl'  •*"xn)  -  0  i-1/  ...,n 


(5.16) 


where  f  (xj_ ,  .  . . ,  xn)  is  the  complex  left  hand  side  of  the 
ith  equation  and  0  denotes  the  complex  number  0  +  jO. 
x  =  (x^,...,xn)  denotes  the  complex  (Real  +  jReal)  problem 
unknowns,  whether  voltages,  currents,  or  parameter  values. 

Expanding  f^  in  a  complex  n-dimensional  Taylor’s  series 
about  an  initial  (complex)  guess  Xq  and  retaining  linear 
terms,  one  gets 


n 

f  i(x1,  .  . . ,  x_)=f  i(xlo,  . .  . ,  xno)  +  Y.  i  (x-i  -  x  )  (5.17) 

j=l  r/xj  5o 


Equating  the  right  side  of  equation  5.17  to  O+jO  (complex) 
yields 


fi(*lo' 


,xno)  +  y  (fh  (Xj 

3=1  G^x  j  2Lo 


xjo> 


(5.18) 


tJ  *  i.r 


i 


xx 

( xlo'  •  •  *  *  xno^  +  X 

3=1  r/x.  x- 


3  I  iio 


(5.19) 


X  (xjo) 

3=1  0^xj  xo 


leading  to  a  set  of  n  simultaneous  complex  linear 
equations 


n  n 

X - -  ^  j  )  =X  — _ -  i  X  jo  )  ”^i  (  xlo»  *  •  •  •  xno  )  i=l»***n  (5.20) 


j=l^/x 


3  |±o 


j=l  cjx j  xQ- 


Equation  (5.20)  maybe  expressed  in  complex  vector /matrix 


form  as 


c(Xo)x  =  d{ xQ) 


(5.21) 


where  C(xQ)  is  an  nxn  matrix  of  complex  coefficents 
dependent  upon  x0=(x|0,...,  xno)  aa  the  point  of 
linearization,  where 


Cij(xo>  =  X  Sill  1 '  3=1  * 

3=1  r/xj  2<o 


(5.22) 


and  ^(x^)  is  a  complex  n-vector  whose  elements  are  also 
dependent  on  Xq; 


86 


n 


i=1d: 


fi 


*o 


(xjQ)  -  ft(x 


lo' 


'xno> 


i=l ,  .  . ,  n 


(5.23) 


From  an  initial  complex  vector  guess  Xq  = 

y  )  at  the  root  of  5.16  may  be  solved  via 
complex  matrix  inversion  with 


x=  (x1<...,xn)»C*1  <*o>d(*o> 


(5.24) 


representing  the  approximate  vector  "root"  of  the  original 
set  of  design  equations  5.16. 

Utilizing  the  approximate  "root"  as  a  new 
linearization  point  for  the  Taylor's  series  expancion,  an 
iterative  process  determining  a  sequence  of  approximations 
to  the  solution  of  5.16  is  obtained: 

*k+l  =  d(xk)  (5.25) 

where  the  elements  of  C(x)  and  d(x)  are  as  given  above. 
Essentially,  the  algorithm  is  a  vector/matr icized 
generalization  of  Newton's  method  to  the 
multi (complex) variable  case. 

Any  of  the  14  design  equations  containing  the  two 
unknown  PI  network  parameters  must  be  linearized  using  the 


procedure  discussed  in  the  previous  paragraphs. 


The  abstract  explanation  of  the  linearization 
given  above  is  complemented  by  a  listing  of  the 
linearization  steps  performed  on  the  nonlinear  equation 
(5.7). 

F  =  V1b'k1  +  k2/XL12^  “  V1BXC12K2  “  v2b(k2/xL12^ 

The  Taylor  Series  Expansion  begins  with  the  step  below. 

c/F  dF  dF  s 

Fo  +  ~J77~  (vlB“vlBo)  +  “ -  (xC12-xC12o)  +  -  (V2B“V2Bo) 

IB  0*C12  tfv  2B 

or 

fq  +  a  +  b  +  c  =  v1 

after  the  partial  derivatives  are  performed  FQ>  A,  B,  and 
C  become 


A  ■  <K1  +  X2/XC1'.  -  k2xC12o><v1B  -  vlBo> 

B  *  (-ViBoK2)(Xc  2  -  Xcl2o) 

C  -  (  •'K2/XL12)  (v- B  -  v2Bo) 

Po=V1Bo<K1  +  K2/XL12>  -  V1BoXC12oK2  '  V2Bo<  K2/XL12 » 


Adding  all  the  terms  toge  .her  (with  a  substantial  number 
of  them  cancelling  each  other  out)  the  remaining  terms 
take  the  following  form 


88 


tit 


IB 


K-- 


K1  “  K2XC12  + 


"LI  2 


+  V. 


-K2 


2B 


XL12 


+  XC12(“V1BoK2>  -  vi  - 


C12ovlBo*2 


Upon  completion  of  the  linearization  process  on  the  4 
nonlinear  equations,  the  14  design  equations  take  on  the 
following  form: 


VlA(Kl-K2*XC13+K2/xLl3)  +  V3A< -K2/XL13 >  =  Vx 
I13"V1A(K1 jxC13“iKl/xL13+3K3)“v3A( DKl/xL13*  =  0 
V1A(“k2/xL13)  +  V3A^K1+K2/XL13“K2*XC31)  =  v3 
I31“V1A(  iKl/XL13)_V3A^Kl3xC3l'’;3Kl/XL13  +  3K3)  ~  0 


(5.27) 

(5.28) 

(5.29) 

(5.30) 


'IB 


K, 


K1  ”  K2XC12  + 


-K, 


+  V 


2B 


LI  2  j 


"LI  2 


+  XC12<-V1BoK2>  =  V1  -  XC12oV1BoK2  (5*31) 


*12  ~  V1B 


jK2 

K1  3xC12o - Z -  +jK3 


XL12 


-  V 


jKi 


2B 


"LI  2 


+  XC12^“3kiv1Bo^  “  "3V1BoXC12oK1  ^5-32^ 


vib(-k2/xli2)  +  v2b(k1+k2/xl12-k2*xc21)  =  v2 

121  ~  VlB^Ki/XL12)-V2B(KijXc21-jK1/XL12+jK3)  =  0 


(5.33) 

(5.34) 


V2C  k1-k2*xC23  + 


*1,23 


+  V3C 


+  XC23^ -V2Cok2)  “  v2  "  xC23oV2CoK2  (5-35) 


3xi  jK^ 

I23  “  V2C  Kl3xC23  "  - -  +  3K3  “  V3C 

*L23  XL23 


+  XC23  ( “ 3x1^2Co  ^  =  “3v2CoXC23o 


(5.36) 


V2C(-K2/XL23)  +  V3c(K1+K2/XIj23-K2*xC32)  =  V3  (5.37) 

132  “  V2C( jKi/xL23)  “  V3c(K1jXc32-jK1/XL23+jK3)  =  0  (5.38) 

The  coefficents  of  these  14  design  equations  are 
calculated  using  the  known  values,  and  these  coefficents 
make  up  the  matrix  A  in  the  formula  x  =  A“1b.  Any 
resulting  constants  are  put  into  the  b  matrix. 


C.  Results 


The  design  equations  have  been  incorporated  into  a 
program  named  DESIGN  which  calculates  the  values  for  the 
XC12  and  XC23  PI  networX  design  parameters.  The  program 
has  been  run  with  a  range  of  input  currents  calculated 
using  the  beam  steering  method  of  Chapter  Three  and  is 
listed  in  Appendix  C.  The  input  current  distribution 


il'iVt 


corresponding  to  a  set  of  ©m  and  ©n  (max  and  null) 
directions  were  calulated  using  the  programs  discussed  in 
Chapter  Two  and  stored  in  a  data  file.  The  design  program 
is  able  to  read  in  these  currents  and  calulate  the 

corresponding  XC12  an<^  XC23  va^-ues* 

Figure  24  shows  the  radiation  patten  of  the 

equilateral  array  when  9m  =  0  and  0n  =  90.  The  current 
distributions  corresponding  to  this  pattern  is  used  as 
input  into  this  example  of  the  linearization  design 
method.  The  initial  guesses  for  the  variables  V2Co,  V1Bo, 

xC12o'  and  xC23o  Were 

V2Co  =  1‘°  +  31-0  V1Bo  =  1.0  +  jl.O 

XC23o  =  1,0  +  31-0  XC12o  =  1,0  +  jl,° 

The  magnitude  and  phase  of  the  currents  flowing 
in  each  element  of  the  array,  the  seven  constant  PI 
network  parameters  input  by  the  user,  and  the  resulting 
two  designed  parameters  Xq^2  and  XC23  are  ^iven  below: 

Magnitude  and  Phase  of  Desired  Current  Distributions 


MIi  = 

1.0000 

=  000.0000 

MI2  = 

-0.5600 

PHi2  =  120.0000 

M 


13 


PH 


13 


0.45 


023.000 


Vs  YYVV*. 


FIGURE  24  —  Radiation  Pattern  Desired 


Given  PI  Network  Parameters 


XL13  = 

10.0 

XC13  =  -1 • 0 

XC31 

11 

1 

Ul 

• 

0 

XL12  = 

10.0 

XC21  =  ~1 • 0 

XC12 

=  Unknown 

XL23  = 

10.0 

XC32  =  “5.0 

XC23 

-  Unknown 

Unknown  Parameters  for  Ten  Iterations 


I 

I 

I 

I 

I 

I 

I 
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1 

XC12 

= 

(3456.3+j234.6) 

XC23 

— 

(00234. 3- j485 . 2 ) 

2 

XC12 

= 

( 0008 . 6- j005 . 1 ) 

XC23 

= 

( -0. 0001+ j . 1098 ) 

3 

XC12 

a 

( -. 0018+ jO. 114 ) 

XC23 

a 

(-0.0001-j.1099) 

4 

XC12 

a 

( 0. 0- j .1103 ) 

^23 

a 

(-0.0001-j.1099) 

5 

XC12 

a 

( 0. 0- j .1103) 

XC23 

a 

(-0.0001-j.1099) 

6 

XC12 

= 

( 0. 0- j .1103 ) 

XC23 

a 

(-0.0001-j.1099) 

7 

XC12 

a 

( 0 . 0- j . 1103 ) 

XC23 

a 

(-0.0001-j.1099) 

8 

XC12 

a 

( 0 . 0- j . 1103 ) 

XC23 

= 

(-0.0001-j.1099) 

The  design  program  was  run  on  a  wide  range  of 
current  distributions.  The  following  results  are  the 


designed  2  and  XC23  and  corresponding  Q's  for  0m  and  @n 
values  ranging  from  ©m  =  0  to  360  in  steps  of  90  degress 
and  0n  =  0  to  360  in  steps  of  45  degrees.  Several  of 
these  runs  are  showns  below  with  specific  9 m  and  ^n 
values,  corresponding  current  magnitudes  and  phases,  the 
designed  PI  network  parameters  and  corresponding  Q  values. 


^  Jl  m.  ft  m  ifc  Vm  ^  ^ 


0.0  on  = 

45.0 

Mn  = 

1.00 

PHn  =  21.0 

MI2  = 

4.24 

phI2  =  11,0 

MI3  = 

-4.57 

PHi3  =22.0 

XC13  ~ 

(3.50E-7 

+  jO . 1103 ) 

qc13  =  314658 

XC23  = 

(9.64E-7 

+  jO . 1105 ) 

QC23  =  H41.9 

0.0  0n  = 

225.0 

MI1  - 

1.00 

PHn  =  -3.0 

MI2  = 

0.47 

PHI2  =  "23*° 

MI3  = 

-0.97 

PHi3  =  33.0 

XC13  ~ 

(2.15E-6 

+  jO . 1103 ) 

qc13  =  51158. 

XC23  = 

(5.18E-5 

+  jO . 1104 ) 

Qc23  =  21321 . 

em  =  9o.o  en  =  135.0 


Mn  = 

1. 

.00 

PHI1 

=  7.0 

MI2  = 

-0. 

.71 

PHI2 

=  -3.0 

MI3  = 

-0, 

.35 

PH^ 

=  43.0 

ft 


I 

I 


1  * ' 


NX 


SC 


xd3  =  (8.77E-6  +  jO .  1103  ) 

X  =  ( -4 . 02E-4  +  jO. 1105 ) 


QC13 


=  12574.0 


Qq23  —  275.9 


From  the  above  example  it  is  obvious  that  the  required 
capacitative  reactances  are  complex.  Therefore  the  Q's  of 
the  designed  capacitors  are  non-infinite. 

In  order  to  check  the  results  of  this  example,  the 
given  and  calculated  PI  network  parameters  are  input  into 
the  analysis  program  of  Chapter  Four.  The  current 
distributions  resulting  from  the  analysis  subroutine  are 
used  to  plot  the  radiation  pattern.  As  indicated  in 
Figure  25  the  pattern  produced  by  the  analysis  run  is 
exactly  the  same  of  the  pattern  desired  (Figure  24). 

However,  if  the  resistive  components  of  X^i2  and  XC23  are 
assumed  zero,  the  currents  produced  by  the  analysis 

subroutine  no  longer  result  in  the  exact  same  pattern  (see 
Figure  26).  The  figure  shows  that  the  null  at  en  =  90  is 
not  as  deep  as  it  was  previously. 


D.  Conclusions 


Surprisingly  (given  the  complexity  of  the  physical 
situation  under  discussion),  the  method  usually  converges. 
As  seen  in  the  above  examples,  if  it  is  going  to  do  so,  it 
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converges  to  single  precision  accuracy  in  about  four 
iterations,  yielding  virtually  instantaneous  response  on  a 
multiuser  VAX  11/780. 

The  main  objective  of  th’s»  chapter  was  to  supply 
the  reader  with  the  essential  components  of  the  extended 
Newton's  Method  linearization  design  method  for  solving 
for  the  parameters  of  the  PI  interconnection  networks. 

The  14  design  equations  were  derived  from  the  24  complex 
equations  describing  the  array  and  their  source  of 
nonlinearity  was  clearly  indicated.  The  linearization 
process,  which  resulted  in  an  iterative  Newton's  Method 
type  procedure,  was  examined  quite  closely.  This  method 
converges  for  a  majority  of  the  patterns  it  is  ask  to 
produce. 

This  design  procedure's  biggest  draw  back  is  the 
fact  that  it  only  designs  two  of  the  nine  PI  network 
parameters.  Also,  if  the  given  parameters  are  specified 
with  magnitudes  larger  than  10.0  the  design  routine  will 
not  converge.  However,  the  fact  that  this  method  works 
for  any  PI  network  parameters  is  in  itself  an  achievement. 
This  procedure  has  proven  (at  least  in  theory)  that  it  is 
possible  to  use  the  interconnection  networks  of  this 
thesis  in  a  phased  array  antenna  system. 
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VI.  A  Direct  Search  Optimization  Method: 

Developement  and  Testing 

This  chapter  presents  a  detailed  examination  of  a 
unique  direct  search  optimization  method  developed  by  the 
author.  The  developement  of  this  direct  search 
optimization  method  was  undertaken  as  a  possible 
alternative  technique  for  solving  the  nonlinear  "design" 
equations  of  the  equilateral  array  described  in  Chapter 
Five.  There  is  definite  application  of  this  method  in  a 
proposed  design  technique  based  on  the  constrained 
optimization  of  the  current  distribution  on  the  array’s 
elements  by  adjustment  of  the  PI  network  parameters  in  the 
"analysis"  equations  of  Chapter  Four.  This  design 
technique  will  be  discussed  in  Chapter  Seven. 

A  brief  discussion  of  the  generalities  of  direct 
search  optimization  methods  begins  this  chapter.  A 
detailed  explanation  of  the  various  blocks  of  decision 
control  for  the  direct  search  optimization  method 
developed  here  is  presented  next.  The  operation  of  each 
of  the  separate  procedures  that  are  performed  on  the 
variables  of  the  objective  function  are  clearly  explained. 
As  an  example,  the  various  procedures  of  the  optimization 
method  are  discussed  in  reference  to  a  two  variable 
quadratic  equation.  This  chapter  concludes  with  sample 


optimization  trials  on  several  quadratic  and 
transcendental  multivariable  objective  functions. 


A.  General  Operation  of  the  Hillclimber 

Direct  search  methods  (  hillclimbing  methods)  of 
optimization  base  their  operation  on  a  finite  number  of 
evaluations  of  trial  solutions  to  an  objective  function 
called  the  hill  height.  Comparision  of  subsequent  trial 
solutions  is  the  basis  for  further  evaluation  of  the 
particular  problem.  The  methods  involve  evaulations  of 
the  objective  function  starting  with  a  given  initial  point 
in  the  particular  n-dimensional  space.  This  initial  point 
is  some  arbitrary  value  of  the  variable  vector  x  =  (x^# 

...»  xn)T  containing  the  arguments  of  the  function  under 
consideration.  These  optimization  methods  usually  only 
require  objective  function  evaluation  and  do  not  use 
partial  derivatives  (as  in  the  steepest  decent 
optimization  methods).  Information  accumulated  as  the 
search  proceeds  is  used  by  some  techniques  in  support  of 
directional  movements  in  the  variables. 

A  two  variable  quadratic  function  of  the  form 

2  2 

f(x,,x9)  =  (x-,  -  a)  +  (x9  -  b) 


(6.1) 


! 

100 

will  be  used  to  help  explain  in  detail  the  direct  search 
optimization  method  to  be  discussed  in  this  chapter.  This 
optimization  method  will  be  referred  to  as  the 
hillclimber  throughout  the  following  discussions.  and 

i  x2  define  the  axis  of  the  two  dimensional  x^X2  plane.  The 

1  a  and  b  represent  the  location  of  the  minimum  of  the 

» 

i  simple  quadratic  function  given  in  equation  6.1.  Figure 

!  27  shows  the  x^X2  Plane  and  the  level  curves  of  the 

function  represented  by  equation  6.1. 

J  The  hillclimber ' s  basic  operation,  in  relation  to 

I 

equation  6.1,  is  the  adjustment  of  the  variables  x^  an(j  X2 

(starting  from  some  initial  point  (xi0,X20)  in  tlie  Plane 
until  the  minimum  of  the  objective  function  is  realized. 

;  The  variables  are  adjusted  one  at  a  time,  in  alternating 

i 

i 

sequence,  with  the  magnitude  of  the  objective  function 
j  being  calculated  and  evaluated  after  each  adjustment  of  a 

variable  (see  Figure  28).  The  magnitude  of  the  objective 
function  is  defined  as  the  hill  height  of  the  function. 

For  the  two  variable  function  of  equation  6.1,  with  the 
minimum  defined  at  a=l  and  b=2  at  an  initial  location 
xlo=5  and  X2Q=5,  the  hill  height  equals 

I 

'  2  2 
;  HH  (Hill  Height)  =  (5-1)  +  (5-2) 

;  =  25 

i 

!  The  hillclimber  adjusts  the  variables  of  a 
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function  in  the  following  manner.  The  first  variable  of  a 
multivariable  sequence  (  x1#  Xr  )  is  moved  in  a 

positive  direction  equal  in  magnitude  to  a  present  step 
size  s^.  in  other  words,  the  initial  given  value  for 

the  variable  x^  is  replaced  by  xlo  +  s1*d1  where  si  is  tlie 
step  size  and  d^  is  the  direction  of  the  variable  x^;  d^ 
is  either  +1  or  -1.  The  hillclimber  has  moved  the 
location  of  the  variables  of  the  two  dimensional  function 
of  equation  6.1  from  (xlo,x2o)  to  the  new  location  (xlo  + 
s1*d1,  x2o).  The  hill  height  of  the  function  is 
calculated,  and  by  the  comparision  of  the  new  hill  height 
to  the  previous  hill  height  HHq  (the  hill  height 
calculated  from  the  initial  values  of  the  variables)  the 
hillclimber  routine  is  able  to  make  basic  adjustment 
decisions  as  to  new  values  for  x-^  and  x2  -  hopefully 
closer  to  the  new  location  of  the  minimum. 

At  some  point  in  the  optimization  process 
(sometime  during  the  first  adjustment  sequence  through  all 
of  the  variables  of  the  objective  function)  the  hill 
height  becomes  the  hill  height  of  least  magnitude  to  date. 
The  hill  height  of  least  magnitude  becomes  the  previous 
hill  height  upon  which  the  basic  adjustment  decisions  are 
made.  This  previous  hill  height  will  be  referred  to  as 
the  lowest  hill  height  to  date  (HH^t(j). 
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B.  Hillclimber ' s  Fundamental  Control  Procedures 


Comparision  of  the  present  and  the  lowest  hill 
heights  to  date  logically  can  result  in  three  possible 
situations  as  a  result  of  a  change  in  one  of  the  x^:  the 
hill  height  could  have  increased  (HH-|_  >  it  could 

have  decreased  (HH^  <  HHj_td),  or  it  could  have  stayed  the 
same  (HH^  =  HHltd)*  What  the  hillclimber  does  in  each  of 
the  above  situations  will  be  discussed  next. 

First,  let  us  assume  a  new  hill  height  has  been 
calculated  and  the  magnitude  is  not  lower  than  the  least 
hill  height  to  date  ( HH^  >  HHltd>*  The  decisions  of  the 
hillclimber  are  as  follows:  since  the  value  of  the  hill 
height  is  now  greater  than  it  previously  was  (because  it 
moved  further  away  from  the  desired  value),  the  variable 
xi  shall  be  returned  to  the  previous  value  x^q  and  t*16 
direction  d^  shall  be  changed  to  the  direction  opposite  to 
what  it  was  previously.  These  two  decisions  constitute 
the  first  fundemental  block  of  decision  control  code  and 
are  shown  more  clearly  in  the  notation  below: 


IF  (HHX  >  HHltd)  THEN 


(returned  to  x^) 


When  the  hillheight  H.  is  less  than  the  least  hill 


r5Tr»jr»>r«>T> jn  r: 


height  to  date  HH^td  ( HH^  <  HHltd),  the  variable  is 
allowed  to  retain  the  new  value  (x^+s^*d^)  and  current 
direction  while  the  succeeding  variable  Xj  is  immediately 
adjusted  to  the  new  value  x2+s2*d2'  This  is  shown  below 
in  notation  form  as 

IF  (  HHx  <  )  THEN 

1  *  X  «  —  v  ■  4>  e  •  . 

1  X1  ^  S1  Q1 


2.  x 


(i+l)  *  x( i+l )  +  s(i+l)  d(i+l) 


In  situations  where  the  present  hill  height  equals 
the  current  least  hill  height  ( HHX  =  HHltd)  the  variable 
x1  is  returned  to  the  previous  value  and  the  present 
direction  is  maintained.  The  succeeding  variable  x2  is 
immediately  adjusted  using  the  current  step  size  and 
direction.  In  notational  form  these  decisions  constitute 
the  third  fundemental  block  of  decision  control  code  and 


IF  (  HHX  =  HHltd)  THEN 

1.  xi  =  xi  "  s^*d^ 

2*  x( i+l )  =  x(i+l)o  +  s( i+l ) *d( i+l ) 

3.  Si  =  Si*10.0 


C.  Hill  Height  Not  Changing  (HH.  ,  =  HH- ) 
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Another  block  of  fundamental  decision  code  is  used 
to  handle  situations  where  adjustment  of  the  variables  of 
the  objective  function  does  not  change  the  magnitude  of 
the  hill  height.  In  other  words  the  hill  height  is  not 
being  changed  within  the  word  length  of  the  computer 
making  the  calculations.  This  can  happen  when  the 
variables  being  adjusted  have  reached  values  that  are  too 
large  or  too  small  to  have  any  effect  on  the  objective 
function,  or  whenever  the  surface  is  flat.  This  problem 
is  detected  by  continually  comparing  each  new  hill  height 
(there  will  be  one  hill  height  for  each  iteration  of  the 
process)  with  the  hill  height  immediately  previous. 

An  example  of  the  operation  of  this  block  of 
decision  code  is  shown  by  defining  (I)  to  be  the  iteration 
counter  of  the  optimization  process.  1=0  corresponds  to 
the  hill  height  calculated  from  the  given  initial  values 
for  the  variables.  In  a  two  variable  optimizaton  process, 
1=7  corresponds  to  the  fourth  adjustment  of  the  first 
variable  x^.  Therefore,  HHQ  is  compared  to  HH^  and  ^or 
the  remaining  iterations  is  compared  to  HH^ 

inorder  to  detect  when  the  adjustment  of  the  variables  is 
not  effecting  the  magnitude  of  the  hill  height.  If  this 
condition  is  detected  the  following  action  is  initiated: 


IF  =  HHi)  THEN 


:.  =  x.  -  s-*d.  (returned  to  x4 ) 


1. 
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2*  X( i+1 )  =  x(i+l)  +  s(i+l)*d(i+l) 

3. 

Now  that  the  four  fundamental  blocks  of  variable 
adjustment  decisions  have  been  discussed,  the  two  variable 
optimization  process  referred  to  in  that  discussion  will 
be  continued  in  some  what  more  detail. 

Using  these  four  fundamental  decision  blocks  the 
two  variables  of  equation  6.1  will  continually  be  moved 
towards  the  desired  values  (x1=a,  x2=b),  resulting  in  a 
hill  height  of  zero  magnitude  (see  Figure  28).  The 
optimization  process  will  continue  the  sequence  of 
alternating  adjustments  until  one  of  the  variables  makes  a 
move  in  a  direction  that  is  further  away  from  the  location 
of  the  minimum  ( HHi  >  HHltd) .  Then  the  hillclimber 
reverses  the  variable's  direction  d^  as  shown  in  the 
preceeding  paragraphs. 


D.  Periodic  Pattern  Detection 


Figure  29  shows  this  happening  to  the  x2  variable 
as  it  moves  below  the  minimum  b  to  an  increased  hill 
height  (HHi  >  HHltd).  At  this  point,  the  x2  variable  will 
be  moved  back  to  the  value  it  had  just  prior  to  the  step 
that  increased  the  hill  height  and  the  direction  d^^ 


back  and  forth  sequence  will  consist  of  three  x2  variable 
magnitudes;  the  x ^  magnitude  corresponding  to  the  lowest 
hill  height  (HHltd),  an  x2  magnitude  corresponding  to  a 
hill  height  lower  in  magnitude  than  HHltd,  and  an  x2 
magnitude  corresponding  to  a  hill  height  greater  in 
magnitude  than  the  Detection  of  this  periodic 

sequence  of  steps  above  and  below  the  present  best 


variable  magnitude  constitutes  the  most  important  block  of 
decision  control  within  this  optimization  routine;  the 
reduction  of  step  sizes  upon  detection  of  oscillatory 
behavior  in  the  adjustment  of  the  variable's  directions 
and  magnitudes. 

The  periodic  behavior  described  in  the  above 
paragraph  happens  because  the  size  of  the  steps  being 
taken  in  the  positive  and  negative  directions  are  too 
large.  For  instance,  if  the  current  best  magnitude  found 
for  x0  was  50,  the  desired  or  minima  value  (b)  equals 


routine  would  continually  jump  from  50  to  60,  from  60  back 
to  50,  and  then  from  50  to  40  and  back  to  50  again.  With 
an  S2  of  magnitude  10  maintained,  the  routine  would 


forever  repeat  this  pattern.  This  periodic  pattern  of 
movements  is  shown  in  Figure  30  for  both  and  x2 
arguments.  However,  the  routine  is  able  to  detect  this 
periodic  behavior  in  both  arguments  and  adjust  the  step 
size  accordingly. 

For  the  two  variable  example  being  discussed,  the 
step  size  S2  be  decreased  from  its  initial  value  of 

10  to  a  value  of  1.0  and  the  routine  will  begin  to  move 
from  the  best  magnitude  50  in  the  positive  direction 
toward  54.3.  But.  once  the  x2  variable  reaches  the 
magnitude  of  54,  a  periodic  movement  above  and  below  this 
new  best  magnitude  for  X2  will  begin  again.  s2  will  again 
be  reduced  by  a  factor  of  10  to  equal  0.1.  The  routine 
will  continue  movement  towards  the  minimum  and  will 
quickly  find  the  desired  value  since  the  routine  is  taking 
steps  in  the  decimal  limits  of  the  minimum  (0.1).  This 


periodic  behavior  will  happen  on  the  x^  variable  as  it 
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nears  the  minima  (a)  and  s^  will  also  be  decreased  to  the 
decimal  accuracy  of  the  minimum. 

The  periodic  pattern  detection  procedure  uses  an 
array  containing  the  previous  values  of  the  variables 


routine.  The  length  of  a  periodic  pattern  is  twice  the 
number  of  variables.  If  (I)  is  defined  as  the  iteration 
counter  of  the  optimization  process  then  the  periodic 
pattern  detector  compare"'  the  2xN  (where  N  is  the  number 
of  arguments  in  f(x^))  sequential  values  thru 

xi ( 1-4)  of  a  variable  with  the  next  2xN  set  of  sequential 
values  xi(I_3)  thru  xi(I)  after  each  iteration  of  the 
optimization  routine.  If  there  are  two  variables  in  the 
function  being  optimized  a  periodic  pattern  four 
iterations  in  length  can  occur.  Detection  is  accomplished 
by  comparing  the  the  present  value  of  the  variable  and  its 
previous  seven  values.  If  x1  is  the  variable  being 
observed  the  comparisions  are  as  follows: 


Does 


Does 

Does 


Does 


Xl(I-0)  equal  X1(I_4) 
X1(I-1)  eSual  X1  ( 1-5 ) 
X1  ( I — 2  )  ec3ual  X1  ( 1-6 ) 
X1  ( 1-3  )  e<?ual  X1  ( 1-7  ) 


If  the  four  statements  above  are  true  then 


1.  s^  =  s./lO.O  (step  size  is  reduced) 


2'  xi  xlhhi 


(x^  is  set  to  best  value) 


where  x^^  equals  the  magnitude  of  the  variable 
corresponding  to  the  lowest  hill  height  found  to  date. 


I. 

i. 


If  all  four  comparisons  are  equal,  a  periodic  pattern  has 
been  detected  and  is  oscillating  around  the  best  value 
it  can  find  with  the  present  s^  magnitude. 


E.  Prevention  of  False  Pattern  Detection 


Because  the  variables  are  reset  to  the  magnitude 
corresponding  to  the  lowest  hill  height  found  so  far  in 
the  process  it  is  possible  that  the  pattern  detector 
could,  depending  on  the  present  magnitude  or  even  the 
minimum  of  the  function,  detect  another  periodic  pattern 
on  the  next  adjustment  of  a  variable  This  pattern 

would  be  false  because  the  routine  had  not  been  given  the 
chance  to  adjust  the  particular  variable  using  the  reduced 
step  size.  Therefore  a  new  periodic  pattern  could  not 
exist  one  iteration  past  the  detection  of  a  previous 
pattern.  Preventing  detection  until  2n+l  iterations 
beyond  the  detection  of  a  previous  pattern  insures  the 
variable  step  sizes  do  not  get  reduced  on  every  iteration. 


F.  Step  Size  Enlargement 


Reduction  of  the  step  size  magnitude  is  imperative 
to  finding  the  actual  minimum  of  a  function.  On  the  other 
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hand  step  size  enlargements,  while  not  an  absolute 
necessity  of  the  optimization  routine,  do  increase  the 
speed  of  convergence  by  allowing  the  variables  to  move  in 
larger  steps  towards  the  desired  results.  The  step  size 
enlargement  module  depends  on  a  user  input  constant  Mdist, 
which  specifies  the  length  of  continuous  movement  in  the 
same  direction.  The  word  "length"  means  the  number  of 
movements  of  a  particular  variable  in  the  same  direction. 
The  Mdist  variable  should  usually  be  specified  to  be  10 
steps  in  the  same  direction,  but  has  been  given  other 
values  in  test  runs  to  see  if  there  is  any  optimal  length 
the  variables  should  be  allowed  to  move  in  the  same 
direction  before  the  step  size  is  enlarged.  For  Lhc 
objective  functions  tested,  values  of  Mdist  anywhere 
between  5  and  13  seem  to  perform  equally  well.  A  length 
of  10  was  chosen  for  all  the  example^  in  this  thesis 
because  all  other  step  size  evaluation  procedures  base 
their  movements,  reductions,  and  comparisons  on  the  factor 
of  10.  This  length  variable  can  be  input  to  a  users  main 
program  and  passed  into  the  optimization  routine  for  the 
purpose  of  tuning  the  hillclimber  to  the  users  particular 
problem. 


G.  Step  Size  Boundaries 


m 


There  are  two  procedures  in  this  direct  search 
optimization  routine  which  prevent  the  step  size  of  the 
variables  from  becoming  too  small  or  from  becoming  too 
large.  These  procedures  are  refered  to  as  the  step  size 
maximum  and  minimum.  The  user  is  able  to  specify  the  step 
size  range  for  all  variables  and  the  step  sizes  are 
prevented  from  going  beyond  or  below  these  limits.  The 
step  size  maximum  procedure  checks  the  step  sizes  against 
the  maximum  magnitude  allowed  and  resets  them  to  a 
magnitude  of  ten  below  the  limit  specified  by  the  user. 
Immediately  after  this  step  size  maximum  procedure,  the 
step  size  minimum  module  checks  the  step  size  of  all  the 
variables  against  the  minimum  allowed.  If  one  of  the 
variable's  step  sizes  is  at  the  minimum  magnitude  it  will 
be  reset  to  a  value  that  is  a  factor  ten  above  the  limit 
and  the  particular  variable  in  question  is  reset  to  the 
best  value  found  (to  date).  The  step  size  can  be 
continually  increased  if  the  variable  continues  to  move  in 
one  direction.  This  can  happen  if  the  movements  of  the 
variable  have  only  a  very  small  effect  on  the  hill  height. 
The  step  size  maximum  procedure  prevents  this  from 
happening . 


H.  Basic  Procedures  Reviewed 


There  are  five  distinct  control  modules  for  this 


optimization  routine  and  they  are  as  follows: 

1.  Direction  Reversal  (due  to  increase  in  hill  height) 

2.  Step  Size  Reduction  (due  to  periodic  pattern  detection) 

3.  Step  Size  Enlargement  (to  increase  convergence  rate) 

4.  Step  Size  Maximum  Limit  (to  prevent  divergence) 

5.  Step  Size  Minimum  Limit  (to  prevent  divergence) 

This  optimization  routine  is  in  the  form  of  a 
subroutine  which  is  called  by  a  main  program.  The  main 
program  is  nothing  more  than  a  medium  for  input ing  the 
initial  values  x^Q  and  step  sizes  s^0  and  transferring 
them  to  the  hillclimber.  The  hillclimber  will  call  a  user 
created  subroutine  containing  the  function  to  be 
optimized.  The  hillclimber  passes  the  current  values  of 
the  variables  to  the  subroutine  FUNC  where  the  hill  height 
is  computed.  Appendix  E  contains  a  listing  of  the 
hillclimber  subroutine  HCSUBR,  the  objective  function 
subroutine  FUNC,  and  the  main  program  MAINHC. 


I.  Results 


The  hillclimber  was  tested  on  several  quadratic 


and  transcendental  functions  of  which  a  two  will  be  shown 
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S  below.  The  examples  will  present  the  form  of  the 

objective  function,  the  initial  values  of  the  variables 
and  their  step  sizes,  the  step  size  maximum  and  minimum 
limits,  along  with  the  an^  tlie  corresponding 

optimized  values  of  the  variables. 

The  first  function  tested  has  the  form, 

F<xl)  =  (xx2  -  78.7)2 


with  the  hill  height  equal  to 

HHi  =  (x,2  -  78.7)2 


I 

* 

! 

! 


The  initial  values  and  other  parameters  of  interest  are 


x 


lo 


10.0 


slo 


1.0 


smax 

smin 


1000.0 

0.0001 


and  the  results  after  67  iterations  are 


xllhh  =  8,8713  HHltd  0,00 

The  table  of  data  shown  below  contains  the  value 
for  the  argument  x and  the  resulting  hill  height  for  the 
first  20  iterations  of  the  hillclimber  on  this  first 


* 


example. 
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I 

X1 

F(x-l ) 

1 

10.0 

453.69 

2 

11.0 

1789.29 

3 

10.0 

453.69 

4 

9.0 

5.29 

5 

9.0 

5.29 

6 

8.0 

216.09 

7 

9.0 

5.29 

8 

10.0 

453.69 

9 

9.0 

5.29 

10 

8.0 

216.09 

11 

9.0 

5.29 

12 

10.0 

453.69 

13 

8.9 

0.26 

14 

8.8 

1.59 

15 

8.9 

0.26 

16 

9.0 

5.29 

17 

8.9 

0.26 

18 

8.8 

1.59 

19 

8.9 

0.26 

20 

9.0 

5.29 

The  second  function  used  to  test  this  hillclimber 
is  of  the  form 


$1 


$ 

% 


I 


SSI 


F(xl,x2)  =  [x22  -  20sin(0.05x1 )  Y  +  O.lCx^  +  x^] 
The  initial  values  and  other  parameters  of  interest  are 


lo  =100.  0 

x2o  =  100.0 

lo  =  10*° 

s2o  =  10.0 

s  .  =  0.0001 

smxn 


and  the  results  after  200  iterations  are 


xlhh  “  °-10  x2hh  “  °*11  HHltd  *  2*3E"2 


J.  Conclusions 


For  the  relatively  simple  functions  tested,  this 
unique  hillclimbing  routine  seems  to  work  quite  well. 

The  initial  guesses  can  usually  be  varied  by  the  user 
until  an  optimal  solution  to  the  function  or  functions  is 
reached.  The  step  size  enlargement  procedure  contains  a 
discretionary  variable  Mdist.  This  variable  can  be 
specified  by  the  user  if  desired  and  therefore  creates 
some  uncertainity  in  this  procedure.  The  author  used  a 
value  of  ten  for  this  variable.  The  four  fundamental 
blocks  of  decision  control  code  and  the  periodic  pattern 
detector  are  the  most  important  parts  of  this  routine  and 
the  logic  behind  them  is  quite  sound. 


it 


Since  developement  of  this  routine  was  only  a  part 
of  my  entire  thesis,  further  testing  and  refinement  of  the 
algorithm  is  in  order.  This  routine  was  use  in  the  beam 
steering  procedure  of  Chapter  Two  and  the  results  were 
quite  satisfactory.  The  radiation  patterns  always  had  the 
main  lobe  fairly  close  to  the  desired  direction.  The 
hillclimber  is  also  implemented  in  a  phased  array  design 
procedure  to  be  discussed  in  Chapter  Seven. 


VII.  Interconnection  Network  Design: 

Direct  Search  Optimization 
of  Constrained  Analysis  Equations 

In  Chapter  Four  quantitative  analysis  of  the  24 
complex  equations  describing  the  physical  and  electrical 
characteristics  of  the  three  element  equilateral  array  was 
performed.  The  analysis  of  the  equilateral  array,  given 
values  for  all  of  the  PI  network  parameters,  results  in 
the  current  distribution  on  the  array.  The  hillclimbing 
technique  presented  in  Chapter  Six  was  developed 
specifically  for  application  in  design  of  interconnection 
networks  for  the  phased  array  system  of  Figure  5.  In  this 
chapter,  the  hillclimber  will  be  used  to  adjust  the  PI 
network  parameters  until  a  desired  current  distribution  is 
realized. 

The  discussion  begins  with  an  explanation  of  how 
the  calculated  and  desired  current  distributions  are 
defined.  How  the  hill  height  is  defined  and  calculated  is 
discussed  next.  All  of  the  procedures  of  this  design 
method  have  been  programmed  in  FORTRAN.  An  explanation  of 
the  structure  and  interaction  of  the  various  subroutines 
is  presented.  The  subroutines  of  this  design  method  are 
pulled  together  in  a  main  program  which  accepts  inputs  and 
provides  the  resulting  output.  The  output  of  this  design 


method  is  the  best  current  distribution  found  after  a 


finite  number  of  iterations  of  the  hillclimber.  Sample 
runs  of  this  design  method  showing  the  resulting  current 
distributions  and  PI  network  parameters  are  presented  at 
the  end  of  this  chapter. 

The  design  method  of  this  chapter  uses  the 
hillclimber  discussed  in  Chapter  Six  to  adjust  the 
parameters  of  the  PI  networks  to  whatever  values  will 
result  in  the  desired  currents  flowing  in  the  elements  of 
the  equilateral  phased  array  of  Figure  5.  The 
quantitative  analysis  of  this  system  of  antennas  has  been 
placed  in  a  subroutine  CURR  which  is  called  by  the 
hillclimber  subroutine  CAP6HC.  The  subroutine  CURR  returns 
the  complex  values  of  the  current  distributions  on  the 
elements  for  the  PI  network  parameters  provided  by  the 
hillclimber . 


A.  Hill  Height  Defined 


The  current  distributions  desired  by  a  user 
(corresponding  to  the  desired  pattern)  are  provided  as 
input  to  this  design  method  and  are  part  of  the  hill 
height.  These  desired  currents  ( I ^ ,  i2,  and  Ij)  are 
normalized  and  used  as  the  minima  of  the  solution  to  the 
linear  equations  describing  the  array.  The  desired 


current  distributions  of  elements  #2  and  #3  are  normalized 
to  the  value  that  corresponds  to  a  current  distribution  on 
element  #1  of  amplitude  1.0  and  phase  0.0  and  take  the 
following  form: 


The  currents  being  calculated  after  each  iteration 
of  the  hillclimber  (b1#  b2,  and  b3)  are  also  normalized  to 
values  with  respect  to  the  current  distribution  on  element 
#1  and  are  shown  below. 


The  PI  network  parameters  are  loaded  with  the 
initial  guess  values  (x^Q)  and  passed  into  the  analysis 
subroutine.  After  the  15  complex  linear  equations  are 
solved,  the  analysis  subroutine  CURR  passes  the  computed 
current  distributions  back  into  the  main  program  where 
they  are  normalized  as  shown  above.  The  hill  height  is 
defined  to  be  the  squared  difference  between  the 
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normalized  desired  currents  (Iln,  i2n,  and  l3n)  and  the 
normalized  calculated  currents  (bln,  b2n,  and  b3n)  and  is 
shown  in  equation  7.1.  The  hill  height  is  calculated 
after  each  adjustment  of  a  PI  network  parameter  during  the 
progression  of  the  optimization  routine.  The  real  and 
imaginary  components  of  these  complex  current 
distributions  are  subtracted  and  the  difference  is  then 
squared.  The  resulting  squared  differences  between  the 
currents  real  and  imaginary  components  are  then  added 
together  and  the  sum  is  the  hill  height  upon  which  the 
hillclimber  will  base  decisions.  The  hill  height  takes 
the  following  form: 

2  2 

HH  =  [Re(Iln)-Re(bln)]  +  [Re( I2n) -Re(b2n) ] 

2  2 

+  CRe(I3n)-Re(b3n)3  +  [ Im( Iln) -Im(bln) 3 

2  2 

+  Clm(I2n)-Im(b2n)]  +  [Im(I3n)-Im(b3nn  (7.1) 

Using  input  statements  in  the  main  program  the 
variable  vector  the  iteration  count  I «  the  number  of 

variables  in  the  objective  function  (Nvar),  the  variable's 
intial  step  sizes  s^  »  and  the  range  of  the  step  sizes 

smax  an<^  smin-  are  a^  inPut  by  user.  With  these  inputs 
from  the  main  program  the  hillclimber  is  able  to  evaluate 


the  hill  height  magnitude  and  use  any  periodic  behavior  in 
the  variables  in  order  to  make  the  decisions  necessary  for 
adjustment  of  the  variables. 

Upon  completion  of  the  specified  number  of 
iterations  the  hillclimber  returns  the  best  calculated 
current  distributions  along  with  the  values  of  the  PI 
network  parameters  that  produce  these  currents.  These 
best  calculated  currents  are  normalized  and  used  as  inputs 
to  a  subroutine  which  will  plot  the  resulting  radiation 
pattern  for  the  triangular  geometry  of  the  array.  Also 
the  hillclimber  can  be  conveyed  variables  that  allow  the 
user  to  observe  the  pattern  being  produced  by  the  best 
current  distributions  found  to  date,  at  whatever  iteration 
count,  during  the  hillclimbers  progression.  Listings  of 
the  program  used  in  this  design  procedure  are  given  in 
Appendix  E. 


B.  Results 

With  only  minor  changes  in  the  programs  listed  in 
appendix  F  the  PI  network  parameters  can  be  adjusted  in 
any  combination  and  sequence.  The  programs  listed  in 
appendix  F  are  designed  for  adjustment  of  the  six 
capacitive  reactances  while  the  three  inductive  reactances 
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are  held  constant.  The  results  to  be  presented  in  this 
chapter  will  also  include  the  problems  of  adjustment  of 
the  three  inductive  reactances  with  the  six  capacitive 
reactances  held  constant  and  adjustment  of  all  nine  PI 
network  parameters.  The  specific  data  items  to  be  listed 
for  the  three  design  problems  being  presented  are  the 
desired  radiation  pattern  and  corresponding  current 
distribution,  the  reactances  of  the  Pi  network  parameters 
held  constant,  the  initial  magnitude  and  step  size  of  the 
parameters  being  adjusted  by  the  hillclimber,  the  best 
current  distribution  resulting  from  a  finite  number  of 
iterations  of  the  hillclimber  and  corresponding  designed 
PI  network  parameters,  and  the  resulting  radiation 
pattern. 

The  current  distributions  corresponding  a 
radiation  pattern  with  @m  =  60  and  0n  =  180  will  be  used 
as  input  for  the  three  design  problems.  This  desired 
radiation  pattern  is  shown  in  Figure  31  and  can  be  used  to 
compared  with  the  radiation  patterns  produced  by  the  three 
design  problems  to  be  presented  below. 

The  f  rst  example  uses  the  hillclimber  to  adjust 
the  three  inductive  reactances  of  the  three  PI  networks, 
given  constant  values  for  the  six  capacitative  reactances, 
until  the  desired  pattern  is  realized.  The  resulting 
current  distributions  and  corresponding  designed 
inductances  are  presented  next  and  the  radiation  pattern 
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is  shown  in  Figure  32. 


Desired  Current  Distribution  in  Magnitude  and  Phase 


MIX  =  1.000 

PHn  - 

0.00 

mi2  =  0.8480 

PH^  = 

-148.05 

Mi3  =  0.2593 

PH13  = 

-42.15 

Given  Capacitive  Reactances 


XC13 

-100 

227 

PF 

Y 

C3 1 

= 

-200 

113 

pF 

XC12 
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-300 

75 

pF 

XC21 
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-400 

57 

pF 

XC23 

= 

-500 

45 

pH 

XC32 

s 

-600 

38 

pF 

Initial  Values  and  Step  Sizes 


XL13  =  xlo  = 

100.0 

Slo  =  10.0 

XL12  =  x3  = 

134.0 

s3o  "  10*° 

XL23  =  X5o  = 

23.0 

s5o  =  1.0 

Hillclimber  Parameters 


FIGURE 


1  =  JUU  °max  =  10000.0  Smin  =  0.00001 
Initial  Value  of  Hill  Height 

HHX  =  2.331952 

Calculated  Currrent  Distribution 


%ln  =  1  •  0000 

PHbin  =  0.000 

Mb2n  =  l-l239 

pHb2n  =  _164-7 

%3n  =  0.2011 

PHb3n  *  23*08 

Resulting  Inductive  Reactances  and  Inductance  at  7  MHz 


<L13  =  291925 

6.64 

mH 

XL12  =  -®5-375 

266 

pF 

«l23  =  1650515 

37.5 

mH 

The  second  design  problem  involves  adjustment  of 
the  six  capacitative  reactances,  while  the  inductances  are 
held  constant.  After  500  iterations  of  the  hillclimber 
the  resulting  radiation  pattern  is  plotted  and  is  shown  in 
Figure  33.  The  following  list  of  data  is  in  the  same 
format  as  the  first  design  problem. 


Desired  Current  Distribution  in  Magnitude  and  Phase 

Mji  =  1.000  PHn  =  0.0^ 

Mi2  =  0.8480  PHX2  =  "148.05 
mi3  =  0.2593  PHi3  =  -42.15 

Given  Inductive  Reactances 


LI  3 
*L12 


=  100  2.27  uH 

=  5  0.1136  uH 

=  20  0.4547  uH 


Initial  Values  and  Step  Sizes 


XC13 

s 

xlo 

= 

12.0 

slo 

= 

10.0 

XC31 

s 

X2o 

= 

34.0 

s2o 

= 

10.0 

XC12 

= 

x3o 

= 

50.0 

s3o 

= 

10.0 

XC21 

= 

x4o 

= 

9.0 

s4o 

= 

1.0 

XC23 

s 

x5o 

= 

78.0 

s5o 

= 

10.0 

XC32 

3 

x6o 

= 

134.0 

s6o 

= 

10.0 

MK 

Sv 

ft 

> 

I 

$ 

f 


B! 

I 


Ki 

At 

I 


Hillclimber  Parameters 


SSWwWJSCiW 


1  -  600  smax  *  10000.0  Smln  =  0.000001 
Initial  Value  of  Hill  Height 

HHi  =  2.294484 

Calculated  Currrent  Distribution 


^bln 

=  1.0000 

PHbln  =  °-00° 

Mb2n 

=  0.9172 

PHb2n  =  -143.9 

MbSn 

=  0.1226 

PHb3n  =  "58-9 

Resulting  Capacitive  Reactances  and  Capacitance  at  7  MHz 


*C13 

s 

154022 

0.147 

PF 

XC31 

= 

0.2181 

104 

nF 

XC12 

= 

9. 59E-2 

237 

nF 

XC21 

= 

1.49E-2 

1.52 

uF 

XC23 

= 

4.5428 

5.00 

nF 

XC32 

= 

0.1136 

200 

nF 

The  final  designed  problem  to  be  presented  uses 


the  hillclimber  to  adjust  all  nine  variables  of  the  PI 
networks  until  the  best  approximation  of  the  desired 
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radiation  pattern  is  realized.  The  resulting  radiation 
pattern  is  shown  in  Figure  34. 


Desired  Current  Distribution  in  Magnitude  and  Phase 


Mj 3  —  1 . 000 

Mi2  =  0.8480 
Mi3  =  0.2593 


PHX1  =  °'00 
PHI2  =  -148.05 
PHi3  =  "42.15 


There  are  no  Given  Reactances 


Initial  Values  and  Step  Sizes 


XC13 

= 

xlo  = 

23.0 

slo  - 

10.0 

XC31 

s 

x2o  = 

34.0 

s2o  = 

10.0 

XC12 

s 

x3o  = 

13.0 

S3o  - 

10.0 

XC21 

s 

x4o  = 

45.0 

s4o  ~ 

10.0 

XC23 

s 

x5o 

56.0 

S5o  = 

10.0 

XC32 

= 

x6o  - 

78.0 

s6o  = 

10.0 

XL13 

= 

x7o  = 

7.0 

s7o  = 

1.0 

XL12 

s 

x8o  = 

132.0 

s8o  = 

10.0 

XL23 

a 

X9o  = 

58.0 

S9o  = 

1.0 

Hillclimber  Parameters 


1  =  900  Smax  =  10000.0  Smin  =  0.000001 


Initial  Value  of  Hill  Height 

=  2.295662 

Calculated  Currrent  Distribution 

%ln  =  1-0000  PHbln  =  0.000 

Mb2n  «  0.8582  PHb2n  =  “148*6 
%3n  =  0.2918  PHb3n  =  -4.127 

Resulting  Capacitive  Reactances  and  Capacitance  at  7  MHz 


C13 

s 

-2.36E-2 

536 

pH 

C31 

SI 

0.18658 

121 

pH 

C12 

= 

-2.93E-3 

666 

pH 

C21 

= 

0.54461 

41.7 

nF 

C23 

= 

2.47199 

9.19 

nF 

C32 

= 

611.242 

37.2 

pF 

LI  3 

= 

-32.495 

699 

pF 

L12 

= 

6.09429 

138 

nH 

T  O  "5 

55 

698.966 

15.9 

uH 

www'IwW®WSn5,i 


Upon  comparing  the  desired  radiations  pattern 
(Figure  31)  with  the  calculated  radiation  patterns 
(Figures  32,33,  and  34)  one  can  see  that  for  this 
particular  set  of  current  distributions  a  fairly  good 
pattern  was  generated  from  all  three  types  of  problems. 

The  adjustment  of  the  three  inductive  reactances  resulted 
in  a  radiation  pattern  with  a  null  direction  that  is  about 
5  degrees  away  form  the  desired  direction.  The  main  lobe 
direction  is  quite  close  to  the  desired  direction  and  the 
large  beamwidth  of  the  main  lobe  is  quite  evident. 

The  problem  which  holds  the  inductive  reactances 
constant  while  adjusting  the  capacitative  reactances 
performed  better  than  the  first  problem.  The  null 
direction  was  within  0.5  degrees  of  the  desired  directions 
and  the  main  lobe  was  clearly  in  the  60  degree  direction. 
The  final  problem's  (adjustment  of  all  nine  PI  network 
parameters)  results  were  nearly  as  good  as  the  second 
problems  results  (Figure  33)  except  that  the  null  depth 
was  not  as  deep  (see  Figure  34). 


C.  Conclusions 


The  design  problems  above  were  all  examined  using 
other  input  current  distributions  corresponding  to  other 


W.‘;V7i 


radiation  patterns.  The  results  were  generally  very  good. 
Excellent  approximation  to  the  desired  pattern  was 
obsevered  in  about  50%  of  the  patterns  desired,  rough 
approximation  was  realized  in  about  30%,  and  about  20%  of 
the  patterns  desired  were  not  able  to  be  realized  at  all. 
The  complexity  of  these-  designed  problems,  the  use  of  the 
unique  optimization  routine  of  Chapter  Six,  and  the  small 
real  components  resulting  from  the  design  method  of 
Chapter  Five  leads  one  to  believe  that  finding  the  desired 
current  distribution  using  purely  reactive  interconnection 
networks  is  indeed  a  difficult  problem. 
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VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  beam  steering  procedure  presented  in  Chapter 
Two  worked  quite  well  on  all  the  radiation  patterns  the 
author  desired.  The  program  MAX_RAD  listed  in  Appendix  A 
was  run  inside  two  Do  Loops  with  one  loop  incrementing  the 
0m  direction  and  the  other  incrementing  the  0n  direction. 
Since  the  null  direction  is  realized  by  solving  a  set  of 
linear  equations  for  the  current  amplitudes  it  will  always 
be  in  the  direction  desired.  However,  additional  nulls 
can  be  produced  in  unspecified  directions,  n  the  other 
hand,  the  current  phases  provide  a  means  of  controlling 
the  direction  of  maximum  radiation. 

The  hillclimber  of  Chapter  Six  performed  a 
maximization  procedure  on  the  current  phases  of  the  three 
elements  in  the  antenna  system  of  Figure  5.  Since  there 
was  just  three  elements  in  the  array  the  potential  degree 
of  maximization  was  minimal  because  the  null  is  already 
defined  and  the  remaining  radiation  always  had  a  very 
broad  beamwidth  and  was  usually  already  fairly  close  to 
the  desired  direction. 

This  beam  steering  procedure  can  be  modified  in 
order  to  handle  as  many  elements  as  desired.  This 
procedure  provided  the  author  with  the  current 
distributions  needed  for  designing  the  interconnection 


medium  of  this  thesis. 


The  analysis  of  the  phased  array  antenna  system 
proposed  in  this  thesis  was  demonstrated  in  Chapter  Four. 
The  reader  was  provided  with  the  electrical  and  physical 
characteristic  to  be  used  by  the  author  in  his 
investigations.  As  was  shown  in  the  results,  after 
assuming  some  reasonable  constant  values  for  the  PI 
network  parameter  and  solving  the  linear  complex  equations 
describing  the  system,  the  resulting  current  distributions 
on  the  array  elements  are  of  no  use  to  a  potential  user 
(see  Figure  24).  However,  the  idea  of  adjusting  the  PI 
network  parameters  until  a  useful  current  distribution  is 
realized  makes  the  analysis  equations  of  Chapter  Four  the 
basis  for  the  design  technique  of  Chapter  Seven. 

In  the  design  of  the  phased  array  antenna  system 
of  Figure  5  one  assumes  the  desired  current  distributions 
is  known  through  some  form  of  beam  steering  technique  as 
demonstrated  in  Chapter  Two.  Therefore,  the  values  of  the 
PI  network  parameters  are  the  unknown  variables  desired  by 
the  designer  of  the  system.  Chapter  Five  discusses  the 
nonlinearity  of  the  systems'  equations  in  this  situation 
and  demonstrated  a  linearization  procedure  and  technique 
for  solving  these  linearized  complex  equations  (see 
equations  5.27-5.40).  The  resulting  PI  network  parameters 
from  this  technique  are  not  very  practical.  The  designed 
parameters  contain  very  small  real  components  and  the 
reactive  components  are  also  quite  small.  The  capacitors 


that  would  correspond  to  the  ohm  values  of  these 
reactances  are  very  large  at  7-Mhz. 

Another  draw  back  to  the  design  technique  is  the 
fact  that  one  must  chose  seven  of  the  nine  PI  network 
parameters.  This  situation  immensely  limits  this 
techniques  fexibility.  Another  problem  with  this 
technique  is  that  the  seven  PI  network  parameters  that 
must  be  specified  by  the  user  can  be  no  greater  in 
magnitude  than  10  ohms  reactive.  Any  time  this  technique 
is  tried  with  parameters  larger  than  10  ohms  in  value,  it 
diverges  (  no  solution  is  obtained). 

However,  of  all  the  radiation  patterns  this 
technique  was  run  with,  using  parameters  of  10  ohms  or 
less,  a  90%  convergence  rate  was  achieved  in  less  than  5 
iterations  of  the  iterative  matrix  solver  subroutine 
DESIGN  listed  in  Appendix  C. 

Chapter  Six  was  completely  devoted  to  the 
development  and  demonstration  of  a  direct  search 
optimization  routine  developed  by  the  author  for 
application  in  a  design  technique  discussed  in  Chapter 
Seven.  The  hillclimber  subroutine  HCS JBR  in  Appendix  D 
was  developed  to  the  point  where  it  now  functions  as  a 
n-variable  optimization  routine  requiring  user  input  of 
the  initial  guess  and  step  size  for  each  variable  along 
with  the  range  of  the  variables  step  sizes  and  the  number 
of  iterations  desired.  Also  a  user  must  create  a 


subroutine  containing  the  functions,  inequalities,  or 
system  to  be  optimized.  An  example  of  this  subroutine's 
form  is  shown  in  Appendix  D  under  the  name  FUNC. 

It  is  the  author's  opinion  that  this  hillclimber 
performs  very  well,  and  with  some  refinement  and  continued 
study,  could  become  even  more  competitive  with  some  of  the 
more  traditional  direct  search  optimization  routines  ( 
Simplex  and  Rosenbrock).  This  routine  has  been  used  on 
quadratic  and  transcendental  functions  of  several 
variables  as  well  as  in  the  beam  steering  procedure  of 
Chapter  Two  and  the  design  technique  of  Chapter  Seven. 

From  observation  of  the  hillclimber ' s  performance  in  the 
design  technique  of  Chapter  Seven  and  on  some  of  the  more 
complex  nonlinear  quardratic  functions,  the  routine  seems 
to  be  dependent  on  a  good  initial  guess  for  the  magnitude 
of  the  variables  in  order  to  reach  an  acceptable  solution. 

Chapter  Seven  developed  and  demonstrated  a 
promising  technique  for  the  design  of  interconnection 
networks  for  the  three  element  phased  array  antenna  system 
around  which  this  thesis  has  evolved.  This  technique  uses 
the  15  complex  linear  equations  of  Chapter  Four  as  the 
system  of  equations  that  when  solved  result  in  the  current 
distributions  on  the  three  elements  for  a  particular  set 
of  PI  network  parameters.  These  current  distributions 
will  be  optimized  to  some  desired  magnitude  and  phase 
using  the  hillclimber  of  Chapter  Six.  Basically  the  user 


inputs  an  initial  guess  for  the  PI  network  parameters  and 
the  hillclimber  adjusts  these  parameters  until  the  desired 
current  distributions  result  from  the  solution  of  the  15 
analysis  equations. 

The  first  form  of  this  design  technique  to  be 
investigated  was  the  six  variable  optimization  of  the  six 
capacitive  reactances  of  the  PI  network  parameters.  The 
idea  behind  this  problem  was  the  practical  considerations 
of  using  adjustable  capacitors  in  conjunction  with  three 
constant  inductors  if  this  proposed  array  was  ever 
constucted.  The  results  from  the  investigation  of  this 
problem  were  mixed.  Given  enough  iteration,  a  reasonable 
approximation  of  the  desired  pattern  was  achieved  for 
about  50%  of  the  radiation  patterns  this  procedure  was 
tried  on.  The  author  observed  that  if  the  0m  and  0n 
directions  were  within  45  degrees  of  one  another,  the 
hillclimber  had  problems  trying  to  produce  a  reasonable 
pattern.  Also  the  hillclimber  tended  to  adjust  the 
capacitative  reactances  to  negative  values  which  indicated 
the  system  wanted  an  inductor  at  that  point  instead  of  a 
capacitor.  The  author  then  put  constraints  on  the  range 
of  the  variables  to  try  and  prevent  this  occurrence.  But, 
this  caused  problems  for  some  examples  in  the  number  of 
iterations  needed  for  convergence.  In  other  problems 
these  magnitude  constraints  even  prevented  some  examples 
that  were  converging  before  the  constraints  from  producing 


usable  patterns  at  all. 


The  design  technique  was  tried  on  several 
combinations  of  PI  network  parameters,  besides  the  three 
problems  discussed  in  Chapter  Seven,  and  the  results  were 
quite  similar  to  the  results  already  shown.  The  design 
problem  which  uses  all  nine  PI  network  Parameters  seems  to 
be  of  most  promise  from  a  point  of  flexiblitiy  in  the 
systems  choice  of  network  configuration.  Even  though  more 
iterations  of  the  hillclimber  are  needed  because  of  the 
number  of  variables  in  the  problem,  the  added  flexibility 
outweighs  this  consideration. 

The  author  suggests  further  investigation  of  the 
design  of  interconnection  networks  for  phased  array 
antenna  system  in  the  following  areas:  formulating 
analysis  equations  for  larger  arrays  to  see  how  lossless 
reactive  networks  function  as  the  interconnection  medium 
on  a  larger  scale;  investigation  of  a  more  suitable 
interconnection  medium,  possibly  some  form  of  reactive 
network  where  the  parameters  can  function  as  both  inductor 
and  capacitor;  and  the  method  of  feeding  this  unqiue  array 
must  be  explored  in  further  detail  because  how  the  system 
is  powered  is  a  very  important  consideration. 
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Program  Array  Pattern  Plotter  (H  it  PAT"  j 
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Version  2.1 

written  by:  Joseph  D. Downs  II 


Date:May  1986 


This  program  allows  a  user  to  input (Nantj  antenna  elements 
at  any  locaion  in  the  XY  plane.  The  spacing  between  each  antenna 
is  also  a  variable  that  can  be  input  by  the  user.  The  user  is 
prompted  for  the  excitation  amplitude  and  phase  for  each  element 
and  the  resulting  radiain  pattern  is  plotted  using  a  call  to  the 
module  PPLOT.FOR. 

dimension  amp <  50) , ph (50), x(50),y(50),xy (360,2), er(50),ei(50),Etot (360) 
real  k,d,p, length, Emax,db 
integer  nant 

open(unit*7,f ile='hilpat' ,status='new’ ) 

Emax=0 . 0 

pi=2.0*asin(1.0) 

k=2.0*pi 

write(5,*)'  Input  length  factor  between  elements  (length)' 
read(6, *) length 

write(5,*) ' Input  #  of  Elements' 
read(6, *)nant 

do  J=l,nant 

write(5,*) ' Input  position  of  each  element  #',J 
read( 6, *)x( J) ,y( J ) 

end  do 
do  J=l,nant 

write(5, *)' Input  amplitude  &  Phase  (deg.)  of  element  #’,J 
read(6,*)amp(J) ,Ph(J) 

Ph(J)=Ph(J)*pi/180.0 

end  do 
do  1*1, 360 

Theta=float( I)*pi/180.0 

do  J=l,nant 

if (X(J) .eq.0.0)then 

fphi=sign(pi/2.0,Y(J) ) 
goto  50 

end  if 

fphi=atan(Y(j)/X(J) ) 
if(X(J) .lt.0.0)fphi=fphi+pi 
D=sqrt(X(J)**2.+Y(J)**2. ) *cos ( f phi -Theta ) 

P=(b/ieny th) *K 
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i 


A 
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wr  i  te ( 7 , * ) ' bxl  =  '  ,  x(  1  )  ,  '  bx2=',x(2),’  bx3=',x(3) 
wr i te( 5 , * ) ' ampl  = ' , amp(  1  ) ,  ‘  amp2= ' , amp( 2 ) , '  amp3= ' , amp( 3 ) 
wr  i te ( 5 , * ) ' bxl  = '  , x(  I ) ,  '  bx2=',x(2),'  bx3=',x(3) 


ph( 1 )=x(1 ) *pi/1 80 . 0 
ph(2)=x(2)*p i/180.0 
ph(3)=x(3)*p i/180.0 


wr i te ( 7 , * ) ' END  — Convergence' 

wr i te ( 7 , * ) ' oerr= ' , oerr , '  phi* ’ ,ph( 1 ) , '  ph2=< ,ph(2) , '  ph3 
write(7,*)'  * 

write(7,*)'  The  Real  &  Imag  currents  equal' 


do  1=1 ,3 


end  do 


curl  23 ( I )=amp( I )*cos(ph( I ) )+cj*ainp( I)*sin(ph( I) ) 
write(7,*)'  The  curr  equals ', curl  23 ( I ) 


write(7,*) '  ' 

call  patt3 (amp, ph, kl , pi , length) 
CLOSE ( 7 ) 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

C  PROGRAM  MAX_RAD 

C 

C  THIS  MAIN  PROGRAM  CALLS  THE  HILCLIMBER  SUBROUTINE  TO 

C  TO  ADJUST  THE  ELECTRICAL  PHASE  OF  THE  ELEMENTS  OF  THE  THREE 

C  ELEMENT  ARRAY  IN  FIGURE  1  UNTIL  THE  RADIATION  ON  THE  DIRECTION 

C  Om  IS  MAXIMIZED.  THE  HILLCLIMBER  SUBROUTINE  MAX_HC  CALLS  THE 

C  THE  SUBROUTINE  FUNC_MAX  WHICH  IN  TURN  CALLS  THE  NULL  AND  FFMAG 

C  ROUTINES.  THE  NULL  MODULE  CALCULATES  THE  AMPLITUDES  OF  THE 

C  CURRENT  CORRESPONDING  TO  THE  RADIATION  PATTERN  WITH  THE  DESIRED 

C  NULL  DIRECTION  On. 


REAL  x(3),dincrx(3),pi ,THmax,THnull , acur , bampl , bamp2 ,bamp3 
REAL  amp( 3 ) ,ph( 3 ) 

complex  cur123(3),cj 

integer  time,IHOW 

open (unit *7, f ile='max_rad' , status=’new ) 

wr ite( 5, * ) 'How  long  will  iterations  for  this  run  (Int)' 
read(6, *) time 


cj*(0.0, 1 .0) 
pi*2. 0*asin( 1.0) 

write(5,*) 'Please  input  guesses  for  Phase  and  step’ 


do  ir=1,3 


read(6, *)x( ir ) ,dincrx( ir ) 

wr ite( 7 , * ) ' Ini t  ial  Guess  and  Step' ,x( ir) ,dincrx( ir) 


end  do 


write( 5, *) ' Input  main  lobe  direction(THmax) ' 
read(6,*)THmax 

write(7,*) 'The  main  lobe  direction(THmax)  is  ' ,Thmax 
THmax=THmax*pi/1 80 . 0 

write(5,*) ' Input  null  direction  (THnull)’ 
read(6,*)THnull 

write(7, *) 'The  null  direction  (THnull)  is',Thnull 
THnull=THnull*p i/1 80.0 

write(5,*)'  when  do  you  want  output  of  iterations  to  begin?' 
read ( 6 , * ) I HOW 


CALL  MAX_HC ( THNULL , THMAX , X , AMP , TIME, DINCRX , IHOW) 

wr i te ( 7 , * ) ' ampl = ' , amp( 1 ) , '  amp2= ' , amp( 2 ) , '  amp3= ' , amp( 3 ) 
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er { J) =amp( J ) *cos(ph( J )+P) 
ei  ( J )=amp( J )*sin(ph(J)+P) 


end  do 


do  J=2,nant 

er(l)=er(l)+er(J) 
ei ( l)=ei ( 1) +ei ( J) 

end  do 


Etot ( I )=Sqrt (er (1) **2.+ei ( 1) **2. ) 
if (Etot ( I ) . gt . Emax )Emax=Etot ( I ) 

end  do 

vrite(5,*)'  Emax=',Emax 


do  1  =  1 , 360 


Theta=float( I ) *pi/180 . 0 
xy( I , 1) =ET0T( I ) *C0S( theta ) 
xy ( I , 2 ) =ETOT( I)*SIN(theta) 


end  do 


call  plot (xy ) 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  SUBROUTINE  MAX  HC 

c 

C  THIS  SUBROUTINE  IS  THE  N- VAR I ABLE  VERSION  OF  THE 

C  THE  HILL  CLIMBER  OPTIMIZATION  ROUTINE  DEVELOPED  BY  THE  THESIS 

C  STUDENT.  THIS  IS  THE  MOST  GENERAL  PURPOSE  ROUTINE  AND  IT 

C  IS  USED  FOR  TESTING  THE  OPTIMIZATION  PROCESS  ON  GIVEN 

C  QUADRATIC,  LOGRYTHMIC,  AND  TRANSCENDENTAL 

C  FUNCTIONS.  IT  IS  USED  FOR  MAXIMIZATION  INSTEAD  OF  MINIMIZATION. 


SUBROUTINE  MAX  HC (THnul 1 , THMAX , x , AMP , t ime , dincrx , ihov) 


REAL  OERR, ERR, what ,acur ,acurlO 

real  mlimit,mlimitl0,Merr(40000) ,Ediff 

real  ar rx (10, 40000), x(10),dx (10), dincrx (10) 

real  xsml(lO) , xlrg(10) 

real  bestx( 10 ), THnul 1 , THMAX, AMP 


integer  ch,vait,vait2,nvar,nvar2, ixdir(10) , idch(10) 
integer  te,td,time,mdist, ihow 


INITIALIZATION  OF  COUNTERS 

ACUR=0 . 001 

MLIMIT=1000.0 

NVAR-3 

MDIST«10*NVAR 

ch=l 

oerr=1.0E-10 

1=0 

K=0 

J=0 

L=0 

M=0 


VARIABLES  USED  IN  STEP  SIZE  RESTRAINT  CODE 


mlimitl0=mlimit*10 . 0 
acurl0=acur/10 . 0 


VARIABLES  OF  THE  ROUTINE  THAT  ARE  BASED  ON  THE  NUMBER 
OF  VARIABLES  BEING  ADJUSTED  BY  THE  PROCESS. 


nvar 2=nvar *2 
wai t=nvar*4 
vait2=nvar*2+l 


TO 


YVfV'rv  ‘fvwtrsxM  irjuT.'T.n  it. 


INITIALIZATION  OF  A  COUNTERS  AND  DIRECTION  INDICATORS  THAT  DESCRIBE 
THE  PARTICULAR  STATE  OF  EACH  VARIABLE  DURING  THE  OPTIMIZATION. 


This  begins  the  Hclimber  by  calling  fuction  to  be  optimized 


Call  Func_MAX(THNULL , THMAX , er r , x , nvar , AMP) 


This  is  iteration  counter  used  in  Hclimber 
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THE  ARRAY  MERR(I)  CONTAINS  ALL  THE  VALUES  OF  THE  HILL  HEIGHT. 
EDIFF  IS  THE  DIFFERENCE  BETWEEN  THE  PRESENT  HILL  HEIGHT  AND  THE 
PREVIOUS  HILL  HEIGHT  AND  OEDIFF  IS  THE  DIFFERENCE  BETWEEN  THE 
PRESENT  HILL  HEIGHT  AND  THE  HILL  HEIGHT  OF  LOWEST  MAGNITUDE  FOUND 
TO  DATE. 


Merr ( I )=err 

Ed if f=abs(Merr( I )-Merr( 1-1) ) 
OEdif f=abs(err-oerr) 


Printing  option  that  prints  out  error, I,  and  Network  prameters. 

i  f ( I . ge . Ihow) then 
wr ite{ 7 , * ) 1  ’ 

WRITE(7,15) I ,err,ch,edif f 

Format(/'  1=  1  , 1 5 , 3x , ' Error  =  ' , E16 . 6 , 3x, ' ch= ' , 12 , 2x, ' Edif f ■ ' , E16 . 6 ) 

write ( 7 , * ) '  ' 

vrite(7,17)x(l),x(2),x(3) 

FormatC  Xl=  '  ,F16. 6, 3x  ,  ’  x2=  '  ,  F16. 6, 3x,  ’  x3*  ',F16.6) 
wr  ite ( 7 , * ) '  ' 

write(7,*)'  Step  sizes  . ' 

write(7, *) '  ' 

write(7,16)dincrx(l),dincrx(2)fdincrx(3) 


I 
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16  Format ( ’ i 1= ’ , F10 . 4 , 2x , ■ i 2  = ’ , F10 . 7 , 2x , * i 3= 1 , F10 . 7 ) 


wr  ite( 7 , * ) ’  1 

write(7,*)*  Direction 
wr  ite( 7 , * ) 1  ’ 


wr i te( 7 , 19)dx(l) ,dx(2) ,dx(3) 

19  Format. ’dl*’ , F4 . 1 , lx , ' d2= ' , F4 . 1 , lx, ‘d3*' , F4 . 1 ) 


wr  i  te ( 7 , * ) '  ' 


endi  f 


This  loop  loads  the  storage  array  with  all  values  of  variables  through 
out  the  iteration  sequence. 


do  ix=l,nvar 


arrx( ix, I )=x( ix) 


end  do 


These  two  do  loops  determine  determine  when  a  variable  needs  to  be 
reduced  in  size  when  it  has  gotten  relatively  close  to  an  answer 
and  begins  jumping  around  it. 


IF( I .GE.wait ) then 


do  ia=l,nvar 


J , K , L ,M  AND  MM  ARE  USED  TO  COUNT  THE  LENGTH  OF 
A  PERIODIC  PATTERN. 


3  =  1  +  1 

k=l  +  wait2 

m=0 

mm=0 


do  ib=l,nvar2 


THIS  STATEMENT  DETECTS  THE  PERIODIC  PATTERN. 


if(arrx(ia,j).eq.arrx(ia,k) )m=m+l 
3  =  j  +  l 
k  =  k+ 1 


end  do 


THIS  STATEMENT  PREVENTS  PERIODIC  PATTERN  FROM  BEING 
DETECTED  ONE  RIGHT  AFTER  ANOTHER. 
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if ( 1 . le. idch( ia) )goto  1 


IF  PATTERN  IS  DETECTED  THEN  THE  STEP  SIZE  OF  THE  PERIODIC  VARIABLE 

IS  REDUCED  IN  MAGNITUDE  BY  A  FACTOR  TEN  AND  THE  MAGNITUDE  OF 

THE  VARIABLE  IS  RESET  TO  THE  BEST  VALUE  FOUND  SO  FAR  IN  THE  PROCESS. 


if (m.eq.nvar2)then 

idch(ia)=l  +  wait2 

x( ia)=bestx( ia) 

dincrx( ia) =dincrx( ia)/10 . 0 

write(7, *) ' *********»*--var iable  is'  ,  ia 

vrite(7,*) 'Decreased  increment  of  X,dincrx= ' ,dincrx( ia) 

vrite(7,*) 'x=  ' ,X( ia) 


1  end  do 


endif 


This  loop  keeps  the  step  size  of  the  variables  from  going  out  of  limit. 


do  if*l,nvar 

if (dincrx( if) . Ie.acurl0)then 
dincrx( if )=acur 

vrite(7,*)  * - ' 

vr ite ( 7 , * ) 'dincrxl  has  gone  to  low,  increased  it  to'.dincrxl 

endif 
end  do 


This  block  of  code  determines  when  to  increases  the  size  (mag.)  of  the 
variable  based  on  how  far  it  has  climber  in  one  direction. 


do  ic=l,nvar 


THIS  STATEMENT  CHECKS  THE  VARIABLES  DIRECTION  INDICATOR 
AND  COUNTS  HOW  MANY  ITERTION  THEY  HAVE  MOVED  IN  THE  SAME 
DIRECTION. 

if(  ixdir(  ic)  .eq.Dthen 

xlrg( ic)=xlrg( ic)+l 

else 


end  i  f 


xlrg( ic)=0 


THIS  STATEMENT  DOES  THE  SAME  AS  THE  ABOVE  FOR  THE 
OPPOSITE  DIRECTION. 


if(ixdir(ic).eq.O) then 

xsml( ic)=xsml( ic)+l 

else 

xsml ( ic ) =0 

end  if 


THE  NEXT  TWO  IF-THEN' S  INCREASE  THE  STEP  SIZE  OF  THE 
VARIABLES  IF  THEY  HAVW  MOVED  CONTINUOUSLY  IN  THE 
SAME  DIRECTION. 


if (xlrg( ic) . eq.mdist)then 

dincrx(ic)*dincrx(ic)*10.0 

Write(7, *) ‘ ####**###### — variable  is' ,ic 

Write(7,*) ' Increased  increment  of  X,dincrx=' ,dincrx( ic) 


end  if 


xlrg( ic)=0 


if  (xsml( ic) .eq.mdist ) then 

dincrx( ic)*dincrx(ic)*10.0 

vrite(7f*) '############  —variable  is’,ic 

write(7,») ' Increaesd  increment  of  X,dincrx=' ,dincrx( ic) 


end  if 


xsml( ic)=0 


enddo 


IF  THE  STEP  SIZE  HAS  GONE  ABOVE  THE  LIMIT  SPECIFIED  BY  THE  USER 
THIS  LOOP  WILL  REDUCE  THE  STEP  SIZE  BY  A  FACTOR  OF  TEN. 


do  if=l,nvar 


if(dincrx(if).ge.mlimitlO)then 


end  if 
end  do 


dincrx( i f ) =ml imit 

write (7,*) 'STEP  SIZE  GONE  ABOVE  LIMIT’ 
wr ite( 7 , * ) '  ' 


This  block  takes  the  first  initial  step  in  a  positive  incrmental 
direction  for  X(l) . 


I F ( I .EQ. 1)THEN 
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END  IF 


x(l)=x(l)  +  dx(l)*dincrx(l) 
wr ite ( 7 , * ) '  • 

WRITE(7,*) ’ Initial  HH  =',err 
if (err . It .oerr ) then 
oerr=err 

endif 
GOTO  100 


c 

c 

c 

c 

c 

c 


C 

C 

C 

C 

C 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 


This  block  of  code  changes  the  variable  back  to  its  previous  value 
if  changing  it  caused  the  hill  height  to  increase.  This  block  also 
adjusts  the  variables  based  on  the  direction  indicators  and  step  size 
magnitudes  by  evaluating  the  magnitude  of  the  hill  height. 


do  id=l,nvar 


THIS  OUTER  IF-THEN  STATEMENT  PERFORMS  THE  VARIABLE  ADJUSTMENTS 
WHEN  THE  ROUTINE  IS  OPERATING  ON  THE  LAST  VARIABLE  IN  THE 
FUNCTION. 


i f (ch.eq.nvar ) then 


THIS  IF-THEN  REDUCES  THE  LAST  VARIABLE  BACK  TO  IT'S 
PREVIOUS  VALUE  IF  THE  PRESENT  HILL  HEIGHT  EQUALS  THE 
PREVIOUS  HILL  HEIGHT.  IT  ALSO  IMMEDIATELY  ADJUSTS 
THE  FIRST  VARIABLE  TO  ITS  NEW  POSITION.  THE  STEP 
SIZE  OF  THE  LAST  VARIABLE  IS  INCREASED  BY  A  FACTOR 
OF  TEN. 


IF(EDIFF.eq.O.O) THEN 

x ( nvar ) «x ( nvar ) -dx ( nvar ) *dincrx ( nvar ) 
x(l)=x(l)+dx(l)*dincrx(l) 

DINCRX(NVAR) *DINCRX(NVAR) *10 . 0 

WRITE(7,*) 'ERR  -  PREVIOUS  ERR  LT  IE-6 , Bx= ' , X ( nvar ) 
WRITE (7,*)'i=',i,'  id*', id,'  step*' ,dincrx(nvar ) 

goto  70 
ENDIF 


THE  SAME  AS  ABOVE  IS  PERFORMED  WHEN  THE  PRESENT 
HILL  HEIGHT  EQUALS  THE  BEST  HILL  HEIGHT  FOUND 
TO  DATE. 


I F ( OEDI FF . eq . 0 . 0 )THEN 

x(nvar)=x(nvar) -dx ( nvar ) *d i ncrx ( nvar ) 
x(l)=x(l)+dx(l)*dincrx(l) 

DINCRX ( NVAR )=DINCRX( NVAR) *10.0 
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c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 


WRITE( 7 ,  * ) ' Oerr  and  Err  are  =  to  e-6' 
wr i te ( 7 , * ) 1 EEEERRRRRRRRRRR  =  OOOOeeerrrr' 

WRITE (7 , * ) ' 1= ' ,  I , '  id=',id,'  step= 1 ,dincrx(nvar ) 

goto  70 


END  IF 


IF  THE  HILL  HEIGHT  IS  LOWER  THEN  THE  LAST  VARIABLE 
IS  KEPT  AT  IT'S  NEW  MAGNITUDE  AND  THE  FIST  VARIABLE 
IS  ADJUSTED. 


if (err.gt.oerr)x(l)=x(l)+dx(l)*dincrx(l) 


IF  THE  HILL  HEIGHT  IS  GREATER  THE  LAST  VARIABLE 
IS  REDUCED  AND  FIRST  VARIABLE  IS  ADJUSTED  IMMEDIATELY. 


i.f  (err .  It  .oerr )  then 

x { nvar ) «x ( nva  r ) -dx ( nvar ) *d i nc  r x ( nva  r ) 
x( 1 ) »x ( I ) +dx(  1 ) *dincrx( 1 ) 

endi  f 

goto  70 
endi  f 


THE  FOLLOWING  IF-THEN  STATEMENTS  ARE  FOR  ALL  OF  THE  ITERATIONS 
BESIDES  ADJUSTMENT  OF  THE  LAST  VARIABLE  (  IN  OTHER  WORDS 

MOVEMENT  OF  THE  VARIABLES  IN  SEQUENCE  1,2,3, .  N,  BUT 

NOT  N.  MOVEMENT  OF  THE  LAST  VARIABLE  (N)  IS  HANDLED  ABOVE. 

CH  TRACKS  WHICH  VARIABLE  IS  CURRENTLY  BEING  ADJUSTED. 


I F ( CH . EQ . I D  .AND.  EDI FF. eq. 0 . 0 )THEN 
TE=ID+1 

x( id)*x( id)-dx( id) *dincrx( id) 
x(te)»x(te)+dx( te)*dincrx(te) 

DINCRX ( ID) -DINCRX ( ID) *10.0 

WRITE(7,*) ’ERR  -  PREVIOUS  ERR  LT  IE-6 , Bx= ' , x( id ) 
WRITE (7, *)'!  =  ', I,'  id=',id,'  step=' ,dincrx(nvar) 


goto  70 
END  IF 


IF(CH. EQ. ID  .AND.  OEDI FF . eq . 0 . 0 ) THEN 


TE= ID+ 1 

x ( id ) =x( id) -dx( id ) *dincrx ( id) 
x( te) =x( te) +dx( te ) *dincrx ( te ) 

DINCRX ( ID)=DINCRX( ID) *10.0 

WRITE ( 7 ,*)' Oerr  and  Err  are  =  to  e-6' 
vrite(7,*) 'EEEERRRRRRRRRRR  =  OOOOeeerrrr' 
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END  IF 


WRITE( 7 , *) *  1  = ' , I , '  id=',id,'  step=’ ,dincrx(nvar) 
goto  70 


IF(ch.eq.id  .and.  err.gt.oerrlthen 
td=id+l 

if ( id.eq.nvar)td=l 
x(td)=x(td)+dx(td)*dincrx(td) 

endif 


If(ch.eq.id  .and.  err . It .oerr ) then 
td=id+l 

x( id)=x( id)-dx( id)*dincrx( id) 
x( td)=x( td)+dx( td) *dincrx( td) 

endif 

end  do 


c 

c  This  loop  changes  the  variables  direction  based  on  the  magnitude  of 

c  the  hill  height.  If  the  variables  adjustment  causes  the  hill  height 

c  not  to  be  reduced  below  the  current  best  hill  height  then  the 

c  particular  variables  direction  indicators  are  reversed  so  the  variable 

c  can  move  in  the  other  direction, 

c 


70  do  ie=l,nvar 

if(ch.eq.ie  .and.  ixdir( ie) .eq.l  .and.  err.lt.oerr)then 
ixdi r ( ie ) =0 
dx{ ie)=-1.0 
goto  20 

endif 

if(ch.eq.ie  .and.  ixdiri ie) .eq.O  .and.  err . It .oerr ) then 
ixdir( ie)=l 
dx(  ie)*=1.0 
goto  20 

endif 


end  do 


This  block  keeps  track  of  which  variable  is  to  be  changed  next, 
c 


i f ( ch. Ie . nvar ) then 

i f { ch . eq . nvar ) ch=0 
ch=ch+l 


3 


20 


endif 


This  block  records  value  of  the  best  variables  and  current  values 
as  the  Hclimber  proceeds  thru  the  climb. 


if (err ,gt .oerr)then 

write(5,*) 'err=' ,err, '  at  I=',I 
oerr=err 

do  ip=l,nvar 

bestx( ip)=arrx( ip, I ) 

end  do 

end  if 


This  tells  machine  to  do  another  iteration 


IF( I .LT. time) GOTO  100 
vrite(7,*) '  ' 

wr i te ( 7 , *) ' HH  (hill  height)  at  end  of  run  «',oerr 

DO  IH-1,3 

X(IH)«BESTX(IH) 

ENDDO 


return 

end 
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SUBROUT I KE  FUNC_MAX ( THNULL , THMAX , H I LLH , X , NVAR , AMP ) 

REAL  X ( 3 ) , p i , THmax, H i 1 1H , THnu 1 1 , Fmax , acur , amp ( 3 ) , ph ( 3 ) 
integer  nvar 
complex  cj,curl23(3) 


c  j  =  ( 0 . 0 , 1 . 0 ) 

pi=2.0*asin(1.0) 

k=2.0*pi 


ph( 1 ) =x{ 1 ) *pi/180 . 0 
ph( 2 ) =x( 2 ) *pi/180 . 0 
ph( 3 ) =x( 3 ) *pi/180 . 0 


call  null (pi ,ph,THnull, amp) 
call  FFmag(pi , amp, ph, THmax, Fmax) 
HillH=Fmax 


RETURN 

END 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

subroutine  null (pi ,PH,THnul 1 , amp) 

real  x ( 3 ) , y ( 3 ) ,ph( 3 ) , amp( 3 ) ,p ( 3 ) 

real  kl ,THnull ,c1 ,c2,x11,x12,pi,x21, x22,denom 

X(1)=0.0 
X(  2 ) =0 . 5 
X(  3 )  =  1 .0 
Y( 1 ) =0 . 0 
Y(2)=0. 8660254 
Y( 3 ) *0 . 0 

kl=  2 . 0*pi/4 . 0 

do  20  J*1 ,3 

if (x( J) .eq.O .0  .and.  y ( J ) . eq . 0 . 0 )goto  20 

P(J)=kl*sqrt(x(J)**2.+y(J)**2. ) *cos ( atan2 ( y ( J) ,x(J) )-THnull) 

20  continue 
P(1 ) =0 . 0 

Cl  =-cos(ph( 1 ) +p( 1 ) ) 

XI 1=cos(ph(2)+p(2) ) 

X12=cos(ph(3)+p(3) ) 

C2«-sin(ph( 1 )+p( 1 ) ) 

X21=sin(ph(2)+p(2) ) 

X22=s in (ph ( 3 ) +p( 3) ) 

denom=X1 1 ‘X22-X21 *X1 2 


if (denom. eq.O. 0) then 

write(5,*)‘  Denom  equals  0.0’ 


end  if 


goto  39 


amp ( 1 ) = 1 . 0 

amp( 2 ) =(c1 *X22-X1 2*c2 ) /denom 
amp( 3 ) = (XI 1 *c2-C1 *X2 1 ) /denom 

39  return 
end 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc ccc 

subrout ine  FFmag (pi , amp , ph,THmax , Fmax ) 

dimension  x(3),y(3),ph(3), amp( 3),p(3),er(3),ei(3) 
real  kl , THmax , Fmax , ptotal , p i 


X(  1  )  =  0.0 
X(  2  )  =0 . 5 
X(3)=1 .0 
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Y  ( 1  )=0.0 
Y(2)=0. 8660254 
Y( 3 ) =0 . 0 

kl  =  2 . 0*pi /4 . 0 

do  30  J=1  ,3 

if (x(J) .eq.O .0  .and.  y ( J ) . eq. 0 . 0 )goto  29 

P(J)=kl*sqrt(x(J)**2.+y(J)**2.) *cos(atan2(y(J ) ,x( J) )-THmax) 

29  p(  1 )  =0 . 0 


er(J)=amp(J)*cos(ph(J)+p(J) ) 
ei ( J) =amp( J ) *sin (ph( J ) +p( J ) ) 

30  continue 


do  J=2 , 3 

er  ( 1  )=er'(  1 )  +er  ( J ) 
ei ( 1 ) =ei ( 1 ) +ei ( J) 

end  do 


Fmax=sqrt(er( 1 )**2.0  ♦  ei(1)**2.0) 

call  Power  in (amp,ph,ptotal) 

Fmax*=Fmax/ptotal 

return 

end 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

subrout i ne  Powe r I n ( amp , ph , pto t a 1 ) 
dimension  amp(3),ph(3) 

complex  cj, 11,12, I3,V1,V2,V3,Zs,Z14,z1,z2,z3 
real  ptotal ,p1 ,p2 ,p3 

cj=(0.0, 1 .0) 

I1*amp( 1 )*cos(Ph( 1 ) )  +  cj *amp( 1 ) *sin(ph( 1 ) ) 

I2=amp( 2 ) *cos ( Ph( 2 ) )  +  cj*amp(2 )*sin(ph(2) ) 

I3=amp< 3 ) *cos (Ph( 3 ) )  +  cj *amp( 3 ) *s in(ph( 3) ) 
c 
c 

Z1 4* (20. 4, -14. 18) 

Zs=  ( 36.5,21 .0) 
c 
c 

V1  =  I1*Zs  +  (  12  +  13) *21 4 
V2=l2*Zs  +  ( 1 1 +1 3 ) *Z1 4 
V3  =  I 3*Zs  +  (  1 1 +12 ) *Z1 4 
c 
c 


pi =real (  conjg(I1)*v1) 
p2-real(conjg( I2)*v2) 
p3= real ( con jg( I3)*v3) 

c  write(7,*),p1=',p1, '  p2=’,p 2,’  p3*',p3 

c 

Ptotal«sqrt( (pi +p2+p3 ) /36 . 5 ) 
c 

c  21 *V1 /1 1 

c  Z2*V2/I2 

c  Z3-V3/I3 

c  write! 7, *)' Driving  point  impedances  equal 

c  vrite(7f *)z1 ,z2,z3 


return 

end 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

PROGRAM  ANAL_DESG 

THIS  PROGRAM  SOLVES  THE  15  LINEAR  EQUATIONS  DESCRIBING 
THE  PHASED  ARRAY  ANTENNA  SYSYTEM  OF  FIGURE  5.  THE  VALUES  FOR 
THE  PI  NETWORK  PARAMETERS  AND  THE  INPUT  CURRENT  ARE  PASSED 
TO  THE  ANALYSIS  SUBROUTINE  CURR  WHICH  RETURNS  TO  THIS  PROGRAM 
THE  CALCULATED  CURRENT  DISTRIBUTIONS  FOR  EACH  ELEMENT  OF  THE 
ARRAY. 


complex  cktl (3),ckt2(3),ckt3(3),B(15),It,ITT,curr1(3), Itdesign 
complex  bd( 14) ,crr1 ( 3 ) 

real  amp<3) ,ph(3) 

open  (unit=7,file=:  a.n.al_desg  '  .status- '  new' ) 

pi  =  2 . 0*asin( 1 . ) 

CKT1 (1 >*(0.0,1 .0) 

CKT1 <  2 )  =  ( 0 . 0,50.0) 

CKT1 (3)* (0.0 ,560.0) 

CKT2 ( 1  )  =  ( 0  . 0,10.0) 

CKT2 ( 3 ) * ( 0 . 0,50.0) 

CKT3(1 )= (0.0, 100.0) 

CKT3(3)=(0. 0,5.0) 

write(5,*) ' Input  value  for  XC12  and  XC13,  complex' 
read(6,*)ckt2(2) ,ckt3(2) 

write(7,*) ’ Input  capacitors  XC12  &  XC23  equal ’ , ckt2 ( 2 ) , ckt3( 2) 
WRITE ( 7 , * ) 1 

write(7,*)'  PI  Network  Parameters' 
write(7,*)  '  ' 

wri te( 7 , 75)aimag(ckt1 ( 1 ) ) , aimag (cktl (2) ) , aimag (ckt 1 ( 3 ) ) 

75  format ( 3x , '  XL13  « ' , 2X,F6 . 2 , '  XC13  =’,2X,F6.2,'  XC31  =',2X,F6.2) 

WRITE( 7 , * ) '  ' 

wr ite( 7 , 76 ) aimag ( ckt 2 ( 1 ) ) , aimag (ckt 2 ( 2 ) ) , aimag (ckt 2 ( 3 ) ) 

76  format ( 3x, '  XL1 2  - ' , 2X,F6 . 2 , '  XC12  =’,2X,F6.2,'  XC21  =',2X,F6.2) 

WRITEt  7 , *  )  '  ' 

wr ite( 7 , 77 ) aimag (ckt 3 ( 1 ) ) , aimag ( ckt  3(2)) , aimag (ckt 3 ( 3 ) ) 

77  format( 3x, '  XL23  =',2X,F6.2,’  XC23  =',2X,F6.2,'  XC32  *',2X,F6.2) 

WRITE ( 7 , * )  '  ' 


wr i te( 5 , * ) ' Input  the  current  at  the  input  of  antenna  1' 
read( 6, * ) It 

wr  i  te ( 7 , * ) ' The  input  current  It  =  ',It 
WRITE ( 7 , * ) ' 
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call  curr (CKT1 ,CKT2 ,CKT3 , I t , B ) 


write(5,*)b( 1 ) ,b(2) ,b(3) 

write(7,*)'  The  Resulting  Current  Distributions  (complex)' 

vr ite( 7 , *  )  '  ' 

wr ite ( 7 , * ) 1 1 1  = ' ,b( 1 ) 

write(7, *) ' 12= 1 , b( 2 ) 

wr ite( 7 , * ) 1 13=' ,b(3) 

WRITE ( 7 , * )  1 

write(7,  *) 'The  Corresponding  Magnitudes  and  Phases  (degrees)' 
wr i te ( 7 , *  )  '  ' 

err  1 ( 1 ) =b( 1 ) 
err  1 ( 2 ) =b( 2 ) 
cr r 1 ( 3 ) =b( 3 ) 

call  mag_ph_norm.(  3  ,  crrl  ,  amp,  ph) 
do  i=1,3 

ph(  i )=ph( i  )*180.0/pi 

end  do 

wr ite( 7 , 80 )amp( 1 ) , amp( 2 ) , amp( 3 ) 

80  format ( 3x, '  [II]  =',2X,F7.2,'  [12]  =',2X,F7.2,'  [ 1 3 ]  =',2X,F7.2) 
WRITE (7, 81  )ph( 1 ) ,ph( 2 ) ,ph( 3 ) 

81  f ormat ( 3x ,  '  PHS1  =',2X,F7.2,'  PHS2  =’,2X,F7.2,’  PHS3  =',2X,F7.2) 


ITT=B( 1 )+B(4)+B(6) 

write(5, *) ' It  =11  +  112  +113  is  equal  to', ITT 
WRITE (7,*)'  ' 


call  design(crr1 ,bd, Itdesign) 

wr i te( 5 , * ) ' XC1 2= ’ , bd( 1 3 ) , '  XC2 3= ' , bd( 1 4 ) 
write (5,*) ' Itdesign=' , Itdesign 

WRITE(7,*) 'Designed  Pi  Network  Parameters' 

write (7,*) 'XC12=' ,bd( 13) , '  XC23= ' , bd ( 1 4 ) 
wr i te ( 7  ,  *  )  ' 

wr i te( 7 , * )  '  1 1  (Input  Current  resulting  from  Design)' 
wr ite( 7 , * ) ' ltdesign= '.Itdesign 

stop 

end 

c 

c 

c 


c 


subroutine  CURR(CKT1 ,CKT2,CKT3, I t , B ) 


nn  nnn  nn  nnn  nnn  non  nnnnnn  nnn 
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This  program  is  the  two-port  analysis  of  a  triangular  three 
element  antenna  system. 

************************************************************* 


declare  all  variables- 


INTEGER  N, IA,M, IB, IJOB, IER 

COMPLEX  A(15,15),B(15), WA (255) ,CKT1 (3) ,CKT2( 3) ,CKT3( 3) ,XL1 3 ,XC13 
COMPLEX  XC3 1 ,Z1 1 , Z1 2 , Z1 3 , IT , C46 , C02 , PI  I , XI , X2 , X3 , X4 , X5 , X6 
REAL  wk< 15) 

- matrix  parameters - 

I  A=  1  5 
IB=  1  5 
N=  1  5 
M=  1 


- mutual  impedances- 

Z1 1  =  ( 36 . 5 , 2 1  . 0  ) 

Z13=(20. 4,-14. 18) 

Z1 2=Z1 3 


-constants  generated  from  separation  distances- 


- change  above  reals  into  complex  #'s  with  zero  reals 

C46=(0. 0,46. 19397663) 

C02=(0.0, .018477591 ) 

PII=( .382683432,0.0) 

- set  constant  matrix  b  back  to  zero - 

do  5  j  =  1  ,15 

B(j)=(0. 0,0.0) 

continue 


-DECLARE  INPUT  CURRENT  IT- 


B  (  1  )  =  IT 


■INSERT  GIVEN  PI  NETWORK  VALUES  OF  REACT ANCE- 


-BEGIN  CALCULATING  THE  MATRIX  "A"  COEFFS- 


XL1 3=CKT1 ( 1 ) 

XC1 3=CKT1 (2) 

XC3 1 =CKT1 ( 3 ) 

CALL  SUBROUTINE  COEFF  FOR  FIRST  TIME 

CALL  COEFF (C4 6 ,C02 ,PI I , XL1 3 , XC1 3 , XC3 1 , XI , X2 , X3 , X4 , X5 , X6 ! 
- INSERT  COEFFS  INTO  MATRIX  "A" - 


nnn 


c 

c 


c 

c 


c 

c 


c 

c 


A( 1 , 1 )  =  ( 1  .  ,0.  ) 
A(1,4)=(1.,0.) 
A( 1 ,6)  =  ( 1 . ,0.  ) 
A(2,2)  =  ( 1  .  ,0. ) 
A(2,5)=( 1  .  ,0.  ) 
A ( 2 , 8 )  =  ( 1 . ,0.  ) 
A(3,3)=(1.,0.' 
A(3,7)  =  ( 1 .  ,0.  ) 
A(3,9)  =  ( 1 . ,0. ) 


A ( 4 , 1 ) =Z1 1 
A( 4 , 2 ) =Z1 2 
A( 4 , 3 ) =Z1 3 
A ( 4 , 1 0 ) =X1 
A( 4 , 1 1 )=X2 


A(5,6)  =  ( 1 . ,0.  ) 
A( 5 , 1 0 ) =X3 
A( 5, 11 )=X4 


A(  6, 1 )=Z12 
A( 6 , 2 ) =Z1 3 
A ( 6 , 3 ) =Z1 1 
A(  6 , 1 0 ) =X2 
A  (  6 , 1 1 ) =X5 


A(7,7)«(1.,0.) 

A( 7 , 1 0) *X4 
A( 7, 1 1 ) »X6 

- SECOND  SET  OF  EQUATIONS 

XL1 3=CKT2 ( 1  ) 

XC1 3*CKT2 ( 2 ) 

XC3 1 =CKT2 ( 3 ) 


CALL  COEFF (C46,C02,PI1 ,XL13, XC1 3 , XC3 1 , XI , X2 , X3 , X4 , X5 , X6 ) 

- LOAD  COEFFS  INTO  MATRIX  "A" - 

A(  8 , 1 )-Z1 1 
A( 8 , 2 ) =Z1 2 
A( 8 , 3 ) =Z1 3 
A(  8 , 1  2 ) =X1 
A ( 8 , 1 3 ) =X2 


A ( 9 , 4 )  =  ( 1  .  ,0.  ) 
A( 9 , 1 2 ) =X3 
A ( 9 , 13) =X4 


a( 10, 1 ) »Z1 2 
A ( 1 0 , 2 ) =Z 1 1 
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C 

C 

C 


C 

C 

c 

c 

c 


c 

c 


c 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 


A( 1 0 , 3 ) =Z1 3 
A( 10, 12)=X2 
A ( 1 0 , 1 3 ) =X5 


A(  11 , 5)  =  ( 1 . , 0.  ) 

A( 1 1 , 1 2) =X4 
A( 1 1 , 1 3 ) =X6 

- THIRD  SET  OF  EQUATIONS 

XL1 3=CKT3 ( 1 ) 

XC13=CKT3(2) 

XC31=CKT3(3) 


CALL  COEFF (C46,C02,PII , XL1 3 , XC1 3 , XC3 1 , X 1 , X2 , X3 , X4 , X5 , X6 ) 

A  ( 1 2 , 1 ) =Z1 2 
A  { 1 2 , 2 ) =Z1 1 
A( 1 2 , 3 )*Z1 3 
A( 12, 1 4 ) =X1 
A( 1 2 , 1 5 ) =X2 


A( 13,8)“( 1 . ,0.  ) 
A( 1 3 , 1 4 ) =X3 
A  ( 1 3 , 1 5) =X4 


A(  1 4 , 1 ) =Z 1 2 
A  ( 1 4 , 2 ) =Z1 3 
A( 1 4 , 3 ) =Z1 1 
A( 1 4 , 1 4 ) =X2 
A ( 1 4 , 1 5) *X5 


A(15,9)  =  (1. ,0.  ) 
A( 1 5 , 1 4 ) =X4 
A( 1 5 , 1 5) *X6 


- SET  FUNCTION  DEFINITION 

I JOB=0 


C 

c  - CALL  MATRIX  INVERSION  IMSL  ROUTINE 

C 

CALL  LEQ2C(A,N, IA,B,M, IB, I JOB , WA , WK , I ER ) 

C 

c  - OUTPUT  CALCULATED  VALUES - 

C 

wr i te ( 5 , * ) '  II  =  ' ,b( 1 ) 
write (5,*/' 12  =  ' , b ( 2 ) 
wr i te( 5 , * ) '13  =  ' ,b(3) 


VLV’.VA'WA’  W  k~.  s.^  CT  v.^  ».'  5OTCT  7?.  k^x."  VI X  TOTETSJV  tc  nvrninvmv.'v-v.v.  -i#<jv 
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RETURN 

END 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


***************************************************** 

THIS  SUBROUTINE  CALCULATES  THE  MATRIX  COEFFS 
***************************************************** 


SUBROUTINE  COEFF (C46 , C02 , PI  I , XL1 3 , XC1 3 , XC3 1 , XI , X2 , X3 , X4 , X5 , X6 ) 

- DECLARE  COMPLEX  VARIABLES - 

COMPLEX  T1 ,T2,X1 . X2 , X3 , X4 , X5 , X6 , C46 , C02 , PI  I , XL1 3 ,XC1 3 , XC3 1 

- ACTUAL  COMPLEX  ALGEBRA  EQUATIONS - 

T1 = (PI 1+ (C46/XL1 3 ) ) 

T2= (C02+ (PI I/XL1 3 ) ) 

XI =- (T1 + (C46*XC1 3 ) ) 

X2= (C46/XL1 3 ) 

X3=-(T2+(PII *XC1 3 ) ) 

X4= (PI I/XL1 3 ) 

X5«-(T1+(C46*XC31 ) ) 

X6=- (T2+ (PI  I *XC3 1 ) ) 


RETURN 

END 


V.VJA.W.  TIW1 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

■C 

C  DESIGN. FOR 

c 

C  THIS  PROGRAM  USES  THE  DESIGN  TECHNIQUE  DISCUSSED  IN 

C  CHAPTER  5  TO  DESIGN  TWO  OF  THE  PI  NETWORK  PARAMETERS  OF  THE 

C  INTERCONNECTION  MEDIUM  IN  FIGURE  5.  THE  USER  IS  PROMPTED  FOR 

C  THE  DESIRED  CURRENT  DISTRIBUTIONS  AND  THE  OTHER  PI  NETWORK 

C  PARAMETERS  ARE  ASSIGNED  SOME  CONSTANT  MADNITUDE.  THE  SUBROUTINE 

C  DESIGN  IS  CALLED  AND  THE  RESULTING  PARAMETERS  ARE  LOADED  INTO 

C  AN  OUTPUT  FILE.  THE  DESIGN  SUBROUTINE  PERFORMS  THE  ITERATIVE 

C  PROCESS  ON  THE  DESIGN  MATRIX  UNTIL  INSTRUCTED  TO  STOP. 


dimension  CKT 1(3) ,CKT2 ( 2 ) , CKT3 ( 2 ) 
real  amp( 3 ) ,ph(  3 ) 

complex  11,12,13, curr ( 3 ) , I t , bd ( 1 4 ) , cckt 1 ( 3 ) , cckt2 ( 3) , cckt3( 3 ) 
open(unit=7,f ile='design’ , status= ' new’ ) 

write(5,*) ’ Input  magnitude  and  phase  of  desired  current  distribution' 
write(7,*) ' Input  magnitude  and  phase  of  desired  current  distribution* 


do  i = 1 , 3 


end  do 


write(5,*) ' Input  amp  and  phase  for  element  #',i 
read ( 6 , * ) amp ( i ) , ph ( i ) 

vrite(7, *) ’For  element  #’,i,'  Mag  i  Phase® ', amp( i ) ,ph( i ) 


call  cmpn(3,amp,ph,curr) 


CKT1 ( 1 )=10.0 
CKT" (^)=1 .0 
CKT1 ( 3 ) =5 . 0 
CKT2( 1 )*CKT1 ( 1 ) 
CKT2 ( 2 ) =5 . 0 
CKT3 ( 1 ) =CKT1 ( 1 ) 
CKT3 ( 2 ) =7 . 5 


call  design (CKT1 ,CKT2,CKT3,curr,bd) 
writet  7, *) '  ' 

vrite(7, *) '  . NOW  PERFORMING  ANALYSIS. 

write (7,*)' 

It® curr ( 1 ) +bd ( 1 ) +bd ( 3 ) 
wr i te ( 7 , * ) ’  lt=’,It 

ccktl (1 ) =cmplx( 0 . 0 , cktl ( 1 ) ) 
cckt 1 ( 2 ) =cmplx ( 0.0, cktl  (2) ) 
ccktl ( 3 ) =cmplx( 0 . 0 , ckt 1 ( 3 ) ) 
cckt 2 ( 1 ) =cmpl x (0.0,ckt2(1)) 


v;,-;.- v.'.-.v.  .v.v;.- 
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176 


T1=aimag(bd( 13)) 

T2=real(bd(13) ) 

cckt2 ( 2 ) =cmplx( 1 1 , t2  ) 

cckt2<3)=(0.0,7.0E-6) 

cckt3 ( 1 )=cmplx( 0 . 0 , ckt3 ( 1 ) ) 

cckt3(2)=cmplx(aimag(bd( 14) ) ,real(bd( 14) ) ) 

cckt3(3)=(0.0,4.0E-7) 

write! 5, *) ' cckt2= ■ , cckt2 ( 1 ) , cckt2 ( 2) , cckt2 ( 3  ) 
write (5,*)  ' cckt3=' ,cckt3(1),cckt3(2),cckt3(3) 
write(5,*)'cckt1=',cckt1(1),cckt1(2) ,cckt1 ( 3 ) 
call  curr 1 (ccktl , cckt2 , cckt3 , 1 1 ) 
c 
c 

c  CKT1 ( 1 ) =CKT1 ( 1 )  +  plus2 

c  i f (CKT1 ( 1 ) .LE. 600.0) go to  10 

c  CKT1 ( 1 )=startC1 

c  CKT1 (2)=CKT1 (2)  +plus1 

c  if (CKT1 (2) .GE.-60.0)goto  10 

c  CKT1 ( 2 ) =startC2 

c  CKT1 (3)=CKT1  (3)  +plus1 

c  I F ( CKT 1(3). GE .-£0.0) GOTO  10 

C 
C 

close(7) 

stop 

end 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  SUBROUTINE  DESIGN 

C 

C  THIS  ROUTINE  IS  AN  ITERATIVE  DESIGN  PROCEDURE 

C  WHICH  USES  GIVEN  CURRENT  DISTRIBUTIONS  FOR  THE  ARRAY  ELEMENTS 

C  AS  INPUT.  IT  CALCULATES  TWO  USER  DESIGNATED  PI  NETWORK 

C  PARAMETERS  THAT  SOLVE  THE  LINEARIZED  SET  OF  14  DESIGN 

C  EQUATIONS. 

C 

C  INPUTS:  CKT 1 , CKT 2, CKT 3, CURR 

C  OUTPUTS :  B 

C 
C 

subroutine  design(ckt1 fckt2,ckt3,curr,bd) 

complex  a(14,14),curr(3),vt(3),I1,I2,I3,b(14),y1,y2,y3,y4 
complex  cl , c2 , c3 , c4 , V2C0 , VI  BO ,GV1 B,GV2C,wa( 224 ) ,bd(14) 
integer  n, ia,m, ib, i job, ier ,Q 

real  wk ( 1 4 ) ,XL1 3 , XC1 3 , XC3 1 , Gxcl 2 , Gxc23 ,x1 ,x2,x3,x4,x5,x6,c46,c02,pii 
dimension  CKT1 ( 3 ) , CKT2 ( 2 ! , CKT3 ( 2 ) 

C 

C  SETS  DIMENSIONS  OF  MATRIX 


i  a  =  1  4 
ib=  1  4 
n=  1  4 
m=  1 
Q=0 


►V*  »  <-*  *. ■  *.*  *. I 


SETS  INITIAL  GUESS  VALUE  FOR  PI  NETWORK  PARAMETERS  AND 
VOLTAGE  VARIABLES. 


Gxcl 2=1  . 

GV1 B= ( 1 . , 1 . ) 
Gxc23=1 . 
GV2C= ( 1 . , 1 . ) 


OUTPUTS  GIVEN  PI  NETWORK  PARAMETERS  TO  DATA  FILE. 

wr i te ( 7 , * ) ' THE  given  capacitor  values  are  XC1 3 , XC3 1 ,XC21 , XC32 ' 
wr  i  te ( 7 , * )CKT1 (2) , CKT1 ( 3 ) , CKT2 ( 2 ) ,CKT3(2) 

WRITE( 7 , * ) 1  The  given  inductor  values  are  equal  to',CKT1(1) 
INITIALIZED  THE  MATRIX'S  OF  COEFFICENTS  AND  CONSTANTS  (A  &  B) . 


do  i = 1 , 14 

do  j=l , 14 

a(  i ,  j  )  =  (0. ,0. ) 

end  do 

end  do 

do  i  * 1 ,14 

b(i)-(0.,0.) 

end  do 


THE  CALL  TO  CMMUL  PERFORMS  THE  COMPLEX  MATRIX  MULTIPLICATION 
BETWEEN  THE  CURRENT  DISTRIBUTIONS  AND  THE  MUTUAL  IMPEDANCES  VALUES 
WHICH  RESULTS  IN  THE  VOLTAGES  AT  THE  BASE  OF  EACH  ARRAY  ELEMENT. 


call  cmmul(curr,vt) 


LOADS  PARAMETERS  OF  THE  FIRST  PI  NETWORK  INTO  VARIABLES  THAT  ARE 
USED  TO  REPRESENT  ALL  OF  THE  PARAMETERS  IN  THE  CALCULATIONS  OF  THE 
COEFFICENTS  OF  THE  UNKNOWN  VARIABLES.  C 


XL1 3=CKT1 ( 1 ) 
XC1 3=CKT1 (2) 
XC31=CKT1 (3) 


CALL  TO  ROUTINE  THAT  PERFORMS  THE  COEFFICENT  CALCULATIONS  FOR  THE 
FIRST  PI  NETWORK. 


call  numb( xl ,x2,x3,x4,x5,x6,c46,c02,pii,XL13,XC13,XC31) 


LOADIND  OF  MATRIX  A  WITH  CALCULATED  COEFFICENTS. 


ft 


555s 


«.■  vr  v 
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a(1 , 2 )  =  ( 1 .0,0.0) 
a(1,5)=(1.,0.) 
a(2,4)=(1.,0.) 
a(2,6)=(1.,0.) 
a ( 3 , 7 ) =cmplx(x1 , 0 . ) 
a(3,8)=cmplx(x2,0. . 
a  ( 4 , 3 )  =  ( 1 . ,  0 . ) 
a(4,7)=cmplx(0. ,x3) 
a{4,b)“cmplx(0. ,x4) 
a( 5 , 7 ) «cmplx( x2 , 0 . ) 
a(5,8)=cmplx(x5,0. ) 
a( 6 , 4 )  =  ( 1 .  ,0.  > 
a( 6 , 7 ) =cmplx( 0 . , x4 ) 
a( 6 , 8 ) =cmplx( 0 . , x6 ) 

C 

C  LOADS  PARAMETERS  OF  THE  PI  NETWORK  CONTAIN I ND  ONE  OF  THE 

C  UNKNOWN  PARAMETERS  INTO  THE  VARIABLES  THAT  ARE  BEING 

C  USED  TO  REPRESENT  ALL  OF  THE  PARAMETERS  IN  THE  CALCULATIONS 

C  OF  THE  COEFFICENTS  OF  SOME  OF  THE  UNKNOWN  VARIABLES. 

C 

XL1 3=CKT2 ( 1 ) 

XC1 3=Gxc1 2 
XC3 1 =CKT2 ( 2 ) 

V1B0=GV1B 

C 

C  CALL  TO  ROUTINE  THAT  PERFORMS  THE  COEFFICENT  CALCULATIONS  FOR  THE 

C  FIRST  PI  NETWORK. 

C 

call  numb(x1 , x2 , x3 , x4 , x5, x6 , c46 , c02 , pi i ,XL1 3 ,XC1 3 ,XC31 ) 

C 

C  CALCULATION  OF  MORE  COEFFICENTS 

C 


yl *- ( VI B0*cmplx(C46 ,0.0)) 
y2=- ( V1B0*cmplx ( 0 . 0 , pi i ) ) 
cl 1=c46*XC1 3 

c1=vt(  1  )-(cmplx(c1 1 , 0.0)*V180) 
c22=pi i*XC1 3 

c2=-(cmplx(0.0,c22)*V1B0) 


C  LOADIND  OF  MATRIX  A  WITH  CALCULATED  COEFFICENTS. 

C 


a ( 7 , 9 ) =cmplx( xl , 0 . ) 
a(7 , 10)=cmplx(x2,0.  ) 
a( 7 , 1 3 )=y1 
a(8 , 1 )  =  ( 1 . , 0 . ) 
a( 8 , 9 ) =cmplx( 0 . , x3 ) 
a (8 , 1 0 ) =cmplx ( 0 . , x4 ) 
a(8,13)«y2 
a(9,9)=cmplx(x2,0. ) 
a( 9 , 1 0 ) =cmplx( x5 , 0  .  ) 


l. 
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a(10,2)=(1.,0.) 
a(  10,9)=cmplx(0. ,x4) 
a  ( 1 0 , 1 0 ) =craplx ( 0 . , x6 ) 


LOADS  PARAMETERS  OF  THE  PI  NETWORK  CONTAIN IND  ONE  OF  THE 
UNKNOWN  PARAMETERS  INTO  THE  VARIABLES  THAT  ARE  BEING 
USED  TO  REPRESENT  ALL  OF  THE  PARAMETERS  IN  THE  CALCULATIONS 
OF  THE  COEFFICENTS  OF  SOME  OF  THE  UNKNOWN  VARIABLES. 


XL1 3=CKT3 ( 1 ) 
XC1 3=Gxc23 
XC3 1 =CKT3 ( 2 ) 
V2C0=GV2C 


CALL  TO  ROUTINE  THAT  PERFORMS  THE  COEFFICENT  CALCULATIONS  FOR  THE 
FIRST  PI  NETWORK. 


call  numb( xl , x2,x3,x4,x5,x6,c46,c02,pii,XL13,XC13,XC31) 


CALCULATION  OF  MORE  COEFFICENTS 


y3=-(V2C0*cmplx(c46,0.0) ) 
y4»- (V2C0*cmplx( 0 . 0 ,pi i ) ) 
c33«c46*XC13 

c3=vt ( 2 ) - (cmplx( c33 , 0.0) *V2C0 ) 
c44=pi i*XC1 3 

c4=- ( cmplx ( 0 . 0 , c44 ) *V2C0 ) 


LOADIND  OF  MATRIX  A  WITH  CALCULATED  COEFFICENTS. 


,11 )=cmplx(x1 , 0. ) 
, 1 2 )=cmplx(x2 , 0 . ) 
, 1 4 ) =y3 
, 5 ) * ( 1 . ,0. ) 

,14) =y4 

,11) =cmplx(0. , x3 ) 
,12) =cmplx( 0, , x4 ) 

1.11 ) =cmplx(x2, 0. ) 

1. 12) =cmplx(x5,0. ) 

l,6)»(1.,0.) 

,,11  )»cmplx(0. ,x4) 
i,  12)=cmplx(0. ,x6) 


CONSTANT  MATRIX  B  FORMATION 

b(  1 )=-curr(2) 
b(2)*-curr(3) 
b(  3  )»vt ( 1 ) 
b(4)=(0. ,0. ) 
b(5)=vt(3) 
b{6)=(0.0,0.) 


rV>JT»jT'->r«>  -J«  >-J< -jt  r 


Dl/|=Cl 

b(  8 ) =c2 
b(  9  )  =vt ( 2 ) 
b(  1  0 )  =  ( 0 . , 0 . ) 
b(  1 1  )=c3 
b(  1  2 ) =c4 
b(  1  3 ) =vt ( 3 ) 
b(  14)  =  (0. ,0.  ) 


CALL  TO  LINEAR  EQUATION  SOLVER  IN  IMSL  SUBROUTINE  LIBRARY, 
i job=0 

call  leq2c(a,n, ia,b,m,ib,i3ob,va,vk,ier) 

Q=Q+1 

write(7,*)'I  (number  of  iterations)  =  ’,Q 
wr ite ( 7 , * ) ' 

wr i te ( 7 , 90 ) AIMAG(B( 13)) ,REAL(B(13) ) ,aimag(b( 14 > ) ,real(b(14) ) 
FORMAT ( ' XC1 2= ' , 2X , F9 . 4 , 1  +j',F9.4,'  XC1 2= 1 , 2X, F9 . 4  , '  +j',F9.4) 

LOADING  THE  NEW  FOUND  VALUES  INTO  NEXT  GUESS  VARIABLES  FOR  NEXT 
ITERATION. 


GV1B*=b(9) 

Gxc12-real(b(13) ) 

GV2C-b(l1 ) 

Gxc23=real(b(14) ) 

if (Q.le. 10)goto  100 

vrite(7,*) 'End  of  iterations' 

do  i  j  =  1  , 1 4 

bd( i j )*b( i j ) 

write(7,*) 'For  ij  Bd( i j )=' ,Bd( i j ) 

end  do 

return 

end 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  SUBROUTINE  NUMB 

c 

C  THIS  ROUTINE  CALCULATES  THE  COEFFICENTS  OF  THE  UNKNOWN 

C  VARIABLES  BASED  ON  THE  INPUT  CONSTANTS.  THE  INPUT  CONSTANTS  ARE 

C  BASED  ON  THE  MAGNITUDE  OF  THE  GIVEN  PI  NETWORK  PARAMETERS. 

C 

C  INPUT:  XL1 3 , XC1 3 , XC3 1 

C  OUTPUT:  XI ,X2,X3,X4,X5,X6 


........ .J 
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c 

subroutine  numb( xl ,x2,x3,x4,x5,x6,c46,c02,pii,XL13,XC13,XC31) 
real  xl ,x2,x3,x4,x5,x6,pi,pii,c46,cQ2,rad,t1 ,  1 1 1 , t2 , t22 , coni 

C 

C  CALCULATIONS  BASED  ON  THE  TRANSMISSION  LINE  LENGTH 

C 

pi=2.0*asin( 1.0) 
rad=3 . 0*pi/8 . 0 
con1=sin(rad) 
pi i=cos( rad) 
c46=50.0*con1 
c02=con1 /50 . 0 


C  CALCULATIONS  BASED  ON  THE  NETWORK  PARAMETERS 

C 

tl =c46/XL1 3 

tl 1 =pi i+tl 

x1=(t1 1-(c46*XC13) ) 

x2  — tl 

t2=pi i/XLI 3 

t22=c02-t2 

x3=-(t22+(pi i*XC1 3 ) ) 
x4=-t2 

x5=(t1 1-(c46*XC31 ) ) 
x6=-U22*(pii*XC31>) 

return 

end 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  SUBROUTINE  CMMUL 

C 

C  THIS  ROUTINE  PERFORMS  THE  COMPLEX  MATRIX  MULTIPLICATION 

C  OF  A  3X3  AND  A  3X1  MATRIX.  IN  PARTICULAR  THIS  ROUTINE  CALCULATES 

C  THE  VOLTAGE  AT  THE  BASE  OF  EACH  ARRAY  ELEMENT.  THE  INPUT  IS  THE 

C  CURRENT  FLOWING  IN  EACH  ELEMENT. 

C 

C  INPUT:  CURR 

C  OUTPUT:  VT 

C 
C 


subroutine  cmmul ( cur r , vt ) 
complex  vt(3),curr(3), impe( 3 , 3 ) , a 
C 

C  MUTUAL  IMPEDANCE  VALUES  FOR  ELEMENT  OF  THE  EQUILATERAL  TRIANGULAR 

C  ARRAY. 

C 


nnn 


c 

impet 1 , 1 )  =  ( 36.5,21 .0) 
imped ,2)»(20.4,-14.18) 
impe( 1 ,3)=(20.4,-14.18) 
impe(2, 1 )=(20.4,-14. 18) 
impe (2,2)=(36.5,21 .0) 
impe (2,3)=(20.4,-14.18) 
impe( 3 , 1 ) = (20 .4 ,-14 . 18 ) 
impe( 3,2)=(20.4,-14.18) 
impe (3,3)— (36.5,21 .0} 


THESE  TWO  DO  LOOPS  PERFORM  THE  MUL I PL I CAT I ON. 


do  k=1  ,  3 

a=(0. 0,0.0) 

do  m=1 , 3 

a=a-*-impe(k,m)  *curr  (m) 
.end  do 

vt(k)*a 

end  do 

return 

end 


a 


K 

$ 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  PROGRAM  MAINHC 

C 

C  THIS  PROGRAM  IS  USED  TO  TEST  THE  HILL  CLIMBER 

C  OPTIMIZATION  ROUTINE  ON  A  VARIETY  OF  FUNCTION  AS  SPECIFIED 

C  BY  THE  USER  IN  THE  SUBROUTINE  FUNC .  THIS  PROGRAM  CALLS  THE 

C  MOST  GENERAL  VERSION  OF  THE  OPTIMIZATION  ROUTINE  AND  PASSES 

C  TO  IT  ANY  VARIABLES  THAT  ARE  INPUT  BY  THE  USER.  THESE  VARIABLES 

C  ARE  USUALLY  THE  NUMBER  OF  VARIABLES  IN  THE  FUNCTION  (NVAR) ,  MAX 

C  AND  MINIMUM  STEP  SIZE  (MLIMIT, ACUR) ,  NUMBER  OF  ITERATIONS  FOR 

C  THE  PARTICULAR  RUN  (TIME),  AND  THE  INITIAL  GUESSES  FOR  THE 

C  VARIABLES  AND  THEIR  INITIAL  STEP  SIZE. 


REAL  X(10) ,Dincrx(10) ,mlimit,mlimitlO,acur, acur 10 

Integer  md ist , nvar , ihow, t ime 

open( un i t=7 ,file='mainhc' , status= 'new' ) 

wr i te ( 5 , * ) 1  input  #  of  variables  in  this  calculation' 
read(6, *)nvar 

wr i te( 5 , * ) *  Input  distance  of  climb  in  one  direction  (mdist)' 
read( 6 , * )mdist 

mdist=10*nvar 

write(5,*) ' Input  starting  values,  and  increments' 
write( 7 ,*)' Initial  Guess  for  Magnitude  and  Step  Size' 
wr i te ( 7 , * )  '  ’ 

do  io=l,nvar 

write(5, *) 'Please  input  guesses  for  variable  #',io 
read( 6 , * ) x ( io ) 

wr i te ( 5 , * ) ' input  starting  increments  for  x  #’,io 
read(6,*)dincrx( io) 

wr i te( 7 , * ) ' Xi= ' , x{ io) , ' Si= ' ,dincrx( io) , '  for  var.  #',io 


end  do 


Input  how  many  iterations  for  the  particular  run 


wr i te ( 5 , * ) ' How  many  iterations' 
read ( 6 , * ) t ime 

wr i te ( 7 ,*)' Number  of  I terat ions= ' , t  ime 
wr i te ( 7 , * ) '  1 

Input  the  acurracy  limit  of  this  calculation 


write(5,*)'  Input  maximun  step  size' 
read( 6, * )ml imit 


c 

c 

c 


wr  ite( 7 ,*)' Maximum  step  size  =',mlimit 
wr i te( 7 , * ) '  ' 

wr i te ( 5 , * ) ' Input  the  minimun  step  size' 

Read (6 , *) ac ur 

wr i te ( 7 , * ) ' Min imun  step  size  =',acur 

wr i te( 7 , * ) ' Var iables  move  in  same  direction  (mdist ) ' , mdist 

wr i te ( b , * ) ' Beg  in  outputing  HC  data  when?' 
read ( 6 , * ) I  how 

CALL  HC( nvar , x , t ime , diner x , mdi st , i how , ml imi t , acur ) 
c lose ( 7 ) 


».« ki.» K«j.l  |.t  1.1 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  SUBROUTINE  HC 

c 

C  THIS  SUBROUTINE  IS  THE  N-VARI ABLE  VERSION  OF  THE 

C  THE  IiLL  CLIMBER  OPTIMIZATION  ROUTINE  DEVELOPED  BY  THE  THESIS 

C  STUDENT.  THIS  IS  THE  MOST  GENERAL  PURPOSE  ROUTINE  AND  IT 

C  IS  USED  FOR  TESTING  THE  OPTIMIZATION  PROCESS  ON  GIVEN 

C  QUADRATIC,  LOGRYTHMIC,  AND  TRANSCENDENTAL 

C  FUNCTIONS. 


SUBROUTINE  HC( nvar , x,time,dincrx,mdist, ihow ,  ml imit,acur) 


REAL  OERR, ERR, what, acur,acurlO 
real  mlimit,mlimitl0,Merr(40000),Ediff 
real  arrx(10, 40000), x(10),dx (10), dincrx ( 10 ) 
real  xsml ( 10 ) , xlrg ( 10 ) 
real  bestx(10) 

integer  ch , wait ,wait2 , nvar , nvar 2 , ixdir(10) , idch(10) 
integer  te, td, time, mdist, ihow 


INITIALIZATION  OF  COUNTERS 


ch  =  l 

oerr=l . 0E30 
1-0 
K=0 
J  =  0 
L=0 
M=0 


VARIABLES  USED  IN  STEP  SIZE  RESTRAINT  CODE 


mlimitl0=mlimit*10.0 
acurl0=acur/10 . 0 


VARIABLES  OF  THE  ROUTINE  THAT  ARE  BASED  ON  THE  NUMBER 
OF  VARIABLES  BEING  ADJUSTED  BY  THE  PROCESS. 


nvar2=nvar*2 
wa i t=nvar *4 
wait2=nvar*2+l 


INITIALIZATION  OF  A  COUNTERS  AND  DIRECTION  INDICATORS  THAT  DESCRIBE 
THE  PARTICULAR  STATE  OF  EACH  VARIABLE  DURING  THE  OPTIMIZATION. 


MSBWBCBPI 
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do  io=l,nvar 


dx( io ) =1 . 0 
ixdir( io)=l 
xsml ( io ) =0 . 0 
xlrg(  io)=0.0 
idch! io)=0 


end  do 


This  begins  the  Hclimber  by  calling  fuction  to  be  optimized 


Call  Func(err,x,nvar) 


This  is  iteration  counter  used  in  Hclimber 


1  =  1  +  1 


i 

i 

% 


$ 

I 

is 


THE  ARRAY  MERRt I )  CONTAINS  ALL.  THE  VALUES  OF  THE  HILL  HEIGHT. 
EDIFF  IS  THE  DIFFERENCE  BETWEEN  THE  PRESENT  HILL  HEIGHT  AND  THE 
PREVIOUS  HILL  HEIGHT  AND  OEDIFF  IS  THE  DIFFERENCE  BETWEEN  THE 
PRESENT  HILL  HEIGHT  AND  THE  KILL  HEIGHT  OF  LOWEST  MAGNITUDE  FOUND 
TO  DATE. 


Merr( I )=err 

Ed if f=abs(Merr( I )-Merr(  I  —  1 ) ) 
OEdif f=abs(err-oerr) 


Printing  option  that  prints  out  error, I,  and  Network  prameters. 

if ( I .ge. Ihowithen 
write( 7 , * )  ’  ’ 

WRlTE(7,15)I,err,ch,Ediff 

Format  < / ’  1=  \I5,3xr'HH=  ’  ,  E12 . 6 , 3x , ’Ch= ' , 12 , 2x, ' HHdi f =  '  ,  E13 . 6 ) 
wr ite( 7  ,  * )  '  ' 

wr i te( 7 , * ) 'Magni tudes  . ' 

write(7,17)x(l),x(2),x(3),x(4) 

Format  (6x,  '  Xl=  '  ,  F12 . 6 , 2x,  '  x2=  '  , F12 . 6 , 2x, ' x3= ' , F12 . 6 , 2x, ' x4= ' , F12 . 6  ) 
wr i te( 7 , * )  '  ' 

write! 7 , 18 ) x( 5 ) ,x(6) ,x(7) ,x(8) 

Format!'  X5=  '  , F16 . 6 , 3x , ' x6= ' , F16 . 6 , 3x , ' x7=  ' , F16 . 6 , 3x , ' x8=  ’,F16.6) 
wr i te ( 7 , * ) x ( 9 ) 
write! ^  *)’  ’ 

write(7,*)'  Step  sizes  . ' 

write!7,16)dincrx(l),dincrx(2),dincrx(3),dincrx(4) 

Format (6x, ' sl= ' , FI  2 . 6 , 2x , ' s2= ' , FI  2 . 6 , 2x , ' s3= ' , F12 . 6 , 2x , ' s4= ' , FI  2 . 6 ) 


I 

I 

to 
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c  write(7,21)dincrx(5),dincrx(6),dincrx(7),dincrx(8) 

c  21  Format ( ' ixl= ' , F14 . 7 , lx, ' ix2= ' , F14 . 7 , lx, ’ ix3= ' , F14 . 7 , lx 
c  $  ,  ' ix4=' ,F14.7) 

wr i te( 7 , *  )  '  ' 

write(7,*)'  Direction  . ' 

wr ite(7 , 19)dx( 1 ) ,dx(2) ,dx(3) ,dx(4) 

19  Format (6x, ' dl= ’ , F4 . 1 , lx, ' d2= ' , F4 . 1 , lx , ' d3= ’ , F4 . 1 , lx , ’ x4= ' , F4 . 1 ) 
wr ite( 7  ,  *  )  1  1 

'  end if 


c 

c 

c 

c 


c 

c 

c 

c 

c 


This  loop  loads  the  storage  array  with  all  values  of  variables  through 
out  the  iteration  sequence. 


do  ix=l,nvar 

arrx ( ix , I ) =x ( ix ) 


end  do 


These  two  do  loops  determine  determine  when  a  variable 
reduced  in  size  when  it  has  gotten  relatively  close  to 
and  begins  jumping  around  it. 


needs  to  be 
an  answer 


IF( I .GE.wait ) then 

do  ia=l,nvar 


J,K,L,M  AND  MM  ARE  USED  TO  COUNT  THE  LENGTH  OF 
A  PERIODIC  PATTERN. 


j  =  l*l 

k=l  +  wait2 

m=0 

mm-0 

do  ib=l,nvar2 


THIS  STATEMENT  DETECTS  THE  PERIODIC  PATTERN. 


if(arrx(ia,j).eq.arrx(ia,k))m=m+l 

j  =  j  +  l 

k=k+ 1 

end  do 


THIS  STATEMENT  PREVENTS  PERIODIC  PATTERN  FROM  BEING 
DETECTED  ONE  RIGHT  AFTER  ANOTHER. 


nnnnn  noon  nn  nnn  nnn  nnnnnn 
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if (l.le.idch(ia) )goto  1 


IF  PATTERN  IS  DETECTED  THEN  THE  STEP  SIZE  OF  THE  PERIODIC  VARIABLE 

IS  REDUCED  IN  MAGNITUDE  BY  A  FACTOR  TEN  AND  THE  MAGNITUDE  OF 

THE  VARIABLE  IS  RESET  TO  THE  BEST  VALUE  FOUND  SO  FAR  IN  THE  PROCESS. 


i f (m.eq. nvar2  >  then 

idch (ia)=l  +  wait2 

x( ia) =bestx<  ia  ) 

dincrx( ia)=dincrx( ia  )/10 . 0 

wr i te ( 7 , * ) ******»*****--variable  is' ,  ia 

wr ite ( 7 ,*)' Decreased  increment  of  X , d i ncrx= ' , d i ncrx ( ia ) 
write(7,*) ’X=  ' , X ( ia) 

endi  f 


1  end  do 


1=1  +  1 


endi  f 


This  loop  keeps  the  step  size  of  the  variables  from  going  out  of  limit, 
do  if=l,nvar 

if(dincrx(if) . le. acurlO ) then 
dincrx( if )=acur 

wr  i  t  e  (  7 ,  * ) ' - ' 

write(7,*) 'dincrxl  has  gone  to  low,  increased  it  to' ,dincrxl 

end  if 
end  do 


This  block  of  code  determines  when  to  increases  the  size  (mag.)  of  the 
variable  based  on  how  far  it  has  climber  in  one  direction. 


do  ic=l,nvar 


THIS  STATEMENT  CHECKS  THE  VARIABLES  DIRECTION  INDICATOR 
AND  COUNTS  HOW  MANY  ITERTION  THEY  HAVE  MOVED  IN  THE  SAME 
DIRECTION. 

i  f  ( i  xd  i  r  ( i  c ) .  eq .  1 )  t  hen 

xlrg( ic)=xlrg(  ic)  +  l 

else 


end  i  f 


xl rg ( ic ) =0 


L*aL'(I  «il  Wl  *rt_ ‘il  i'l,i  |>_ 


enddo 


THIS  STATEMENT  DOES  THE  SAME  AS  THE  ABOVE  FOR  THE 
OPPOSITE  DIRECTION. 


if (ixdir( ic) . eq.0)then 

xsml (ic)=xsml(ic)+l 

else 

xsml ( ic ) =0 

endi  f 


THE  NEXT  TWO  IF-THEN' S  INCREASE  THE  STEP  SIZE  OF  THE 
VARIABLES  IF  THEY  HAVW  MOVED  CONTINUOUSLY  IN  THE 
SAME  DIRECTION. 


if (xlrg( ic) . eq .mdist ) then 

dincrx(ic)=dincrx(  ic)*10.0 

Write (7,*) * ### ######### -var i able  is’  ,  ic 

Wr i te ( 7 ,*)' Increased  increment  of  X,dincrx= ' ,dincrx( ic) 


endi  f 


xlrg( ic)=0 


i  f ( xsml ( ic) .eq.mdi st ) then 

dincrx(ic)=dincrx(ic)*10.0 

write<7,*)  '############  --variable  is',ic 

wr i te( 7 , * ) ' Increaesd  increment  of  X,dincrx= ' ,di ncrx ( ic ) 


xsml ( ic ) =0 


endi  f 


IF  THE  STEP  SIZE  HAS  GONE  ABOVE  THE  LIMIT  SPECIFIED  BY  THE  USER 
THIS  LOOP  WILL  REDUCE  THE  STEP  SIZE  BY  A  FACTOR  OF  TEN. 


do  if=l,nvar 


if(dincrx(if).ge.mlimitlO)then 


endi  f 


end  do 


dincrx ( i f ) =ml imi t 

vrite(7, *) 'STEP  SIZE  GONE  ABOVE  LIMIT' 
wr ite(7  ,  *  )  '  ' 


This  block  takes  the  first  initial  step  in  a  positive  incrmental 
direction  for  Xil). 


I F ( I .EQ. 1 (THEN 


mm 

•wyTMKsc 


x(l)=x(l)  +  dx( 1) *dincrx ( 1 ) 
if(err.lt.oerr)oerr=err 


GOTO  100 


END  IF 


This  block  of  code  changes  the  variable  back  to  its  previous  value 
if  changing  it  caused  the  hill  height  to  increase.  This  block  also 
adjusts  the  variables  based  on  the  direction  indicators  and  step  size 
magnitudes  by  evaluating  the  magnitude  of  the  hill  height. 


do  id=l,nvar 


THIS  OUTER  IF-THEN  STATEMENT  PERFORMS  THE  VARIABLE  ADJUSTMENTS 
WHEN  THE  ROUTINE  IS  OPERATING  ON  THE  LAST  VARIABLE  IN  THE 
FUNCTION. 


i  f ( ch . eq. nvar ) then 


THIS  IF-THEN  REDUCES  THE  LAST  VARIABLE  BACK  TO  IT'S 
PREVIOUS  VALUE  IF  THE  PRESENT  HILL  HEIGHT  EQUALS  THE 
PREVIOUS  HILL  HEIGHT.  IT  ALSO  IMMEDIATELY  ADJUSTS 
THE  FIRST  VARIABLE  TO  ITS  NEW  POSITION.  THE  STEP 
SIZE  OF  THE  LAST  VARIABLE  IS  INCREASED  BY  A  FACTOR 
OF  TEN. 


IF(EDIFF. eq. 0 . 0 )THEN 

x(nvar) =x( nvar ) -dx ( nvar ) *dincrx ( nvar ) 
x( 1 )=x( 1 )  -t-dx ( 1 )  *dincrx ( 1 ) 

D I NCRX ( NV AR ) =  D I NCRX ( NV AR ) * 1 0 . 0 

WRITE(7,*) 'ERR  -  PREVIOUS  ERR  LT  IE-6 , Bx= ’ r x ( nvar ) 
WRITE (1 ,*) ' l-' , l , '  id=',id,'  step= ' ,dincrx ( nvar ) 

goto  70 
END  IF 


THE  SAME  AS  ABOVE  IS  PERFORMED  WHEN  THE  PRESENT 
HILL  HEIGHT  EQUALS  THE  BEST  HILL  HEIGHT  FOUND 
TO  DATE. 


I F ( OED I FF . eq . 0 . 0 ) THEN 

x( nvar ) =x ( nvar )-dx(nvar)*dincrx(nvar) 
x  ( 1 )  =  x ( l)+dx(l)*dincrx(l) 

D I NCRX ( NVAR )=D I NCRX (NVAR) *10. 0 

WRITE ( 7 , * ) ' Oer r  and  Err  arc  =  to  e-6' 
wr i t  e ( 7 , * ) ' EEEERRRRRRRRRRR  =  OOOOeeerrrr' 


f 
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wKirtsw  ,  - )  •  i=  •  ,  i ,  •  ia=  ’  ,  id,  •  step=  '  .dincrxlnvar  ) 
goto  70 


IF  THE  HILL  HEIGHT  IS  LOWER  THEN  THE  LAST  VARIABLE 
IS  KEPT  AT  IT'S  NEW  MAGNITUDE  AND  THE  FIST  VARIABLE 
IS  ADJUSTED. 


i f ( err . It .oerr ) x( 1 ) =x( 1 ) +dx( 1 ) *dincrx( 1 ) 


IF  THE  HILL  HEIGHT  IS  GREATER  THE  LAST  VARIABLE 
IS  REDUCED  AND  FIRST  VARIABLE  IS  ADJUSTED  IMMEDIATELY. 


if(err.gt.oerr)then 

x(nvar)=x(nvar)-dx(nvar)*dincrx( nvar ) 
x( 1 )*x( 1 ) +dx( 1 ) *dincrx ( 1 ) 

endi  f 


goto  70 
endi  f 


THE  FOLLOWING  IF-THEN  STATEMENTS  ARE  FOR  ALL  OF  THE  ITERATIONS 
BESIDES  ADJUSTMENT  OF  THE  LAST  VARIABLE  (  IN  OTHER  WORDS 

MOVEMENT  OF  THE  VARIABLES  IN  SEQUENCE  1,2,3, .  N,  BUT 

NOT  N.  MOVEMENT  OF  THE  LAST  VARIABLE  (N)  IS  HANDLED  ABOVE. 

CH  TRACKS  WHICH  VARIABLE  IS  CURRENTLY  BEING  ADJUSTED. 


IF(CH.EQ. ID  .AND.  EDI FF . eq . 0 . 0 ) THEN 
TE= I D+ 1 

x ( id ) =x ( id) -dx ( id ) *d incrx ( id ) 
x(te)=x(te)+dx(te)*dincrx(te) 

DINCRXi ID) =DINCRX( ID) *10.0 

WRITE< 7 , * ) ' ERR  -  PREVIOUS  ERR  LT  IE-6 , Bx= ' , x ( id ) 
WRITE( 7, *)'!=', I,'  id=',id,'  step= ', dincrx ( nvar ) 

goto  70 


I F ( CH . EQ . I D  .AND.  OEDI FF . eq . 0 . 0 ) THEN 
TE= I D+ 1 

x(  id)=x(  id ) -dx ( id)*dincrx(  id) 
x(te)=x(te)+dx(te)*dincrx(te) 

DINCRX ( I D ) =DI NCRX ( ID) *10.0 

WRITE(7, *) 'Oerr  and  Err  are  =  to  e-6' 
vr i te( 7 , * ) ’ EEEERRRRRRRRRRR  =  OOOOeeerrrr' 

WRITE ( 7 , * ) ' I = ' , I , '  id=',id,'  step= ’ , d i ncr x ( nvar ) 
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END  IF 


goto  70 


IF(-;..eq.id  .and.  err .  1 1 .  oerr )  then 
td=id+l 

i  f ( id . eq . nvar ) td=l 
x(td)=x(td)+dx(td)*dincrx(td) 

endi  f 


If(ch.eq.id  .and.  err . gt . oerr ) then 
td= id+ 1 

x(  id )  =  x ( id)-dx( id)*dincrx( id) 
x(td)=x(td)+dx(td)*dincrx(td) 

endi  f 


end  do 


This  loop  changes  the  variables  direction  based  on  the  magnitude  of 
the  hill  height.  If  the  variables  adjustment  causes  the  hill  height 
not  to  be  reduced  below  the  current  best  hill  height  then  the 
particular  variables  direction  indicators  are  reversed  so  the  variable 
can  move  in  the  other  direction. 


do  ie=l,nvar 

if(ch.eq.ie  .and.  ixdir( ie) .eq. 1  .and.  er r . gt . oerr ) then 
ixdi r ( ie) =0 
dx( ie)=~1.0 
goto  20 

endi  f 

if(ch.eq.ie  .and.  i  xd i r  ( ie )  .  eq .  0  .and.  er r . gt . oerr ) then 
ixdir ( ie)=l 
dx(  ie)  =  1.0 
goto  20 

endi  f 

end  do 


This  block  keeps  track  of  which  variable  is  to  be  changed  next. 


i  f ( ch . Ie . nvar ) then 

i  f ( ch . eq. nvar ) ch=0 
ch=ch+l 

endi  f 


V 


This  block  records  value  of  the  best  variables  and  current  values 
as  the  Hclimber  proceeds  thru  the  climb. 


if(err.lt.oerr)then 

wr  i  te ( 5 , * ) ' err= ' , err , *  at  I = ’ , I 
oer r=er r 

do  ip=l,nvar 

bestx(ip)=arrx(ip,I) 

end  do 


end  i  f 


This  tells  machine  to  do  another  iteration 


IF< I . LT . t ime ) GOTO  100 
write (7,*) ’  END  OF  RUN' 
wr  i  te ( 7 , * ) '  ' 

write( 7 ,*) 'Values  of  the  Variables  corresponding  to  HHltd’ 
write(7r*)'Xl=',bestx(l),,X2=',bestx(2),,B3=,,bestx(3) 
write(7,*) 'Hill  Height  =',oerr 

return 

end 
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SUBROUTINE  FUNC 

THIS  SUBROUTINE  13  CALLED  BY  THE  HILL  CLIMBER  OPTIMIZATION 
ROUTINE.  IT  CALCULATES  THE  MAGNITUDE  OF  THE  HILL  HEIGHT  BASED 
ON  WHATREVER  FUNCTION  IT  CONTAINS.  THE  FUNCTIONS  CONTAINED 
IN  THIS  SUBROUTINE  ARE  USE  TO  TEST  THE  PERFORMANCE  OF  THE 
THE  OPTIMIZATION  ROUTINE  DEVELOPED  BY  THE  THESIS  STUDENT. 


subroutine  func(err,x,nvar) 

integer  nvar 
real  x(nvar),err 
pi-3.141593 
q-x(l) 

err-  ( x( 2 )-20 . 0*SIN( 0 . 05*x{ 1 ) ) ) *  *  2  +  0.1*(x(2)**2  +  x(l)**2) 
c  err- ( exp ( -x ( 1 ) )  —  1 0 . O*cos( q) +0 . 5*sin(q )-X( 2 ) **  2 . +exp( -x( 2) ) ) *  *  2 . 

c  err- ( q*  *  5-4 . 5*q**4+4 . 55*q**3+2.675*q**2~3. 3*q~l . 4375-x(2) )  *  *  2 

c  err-100* (x( 2)  -  x(l)**2)**2  +  (l-x(l))**2 

c  err=(x(l)**2  -  78.7)**2 

c  ERR  -  ERR1  +  ERR2  +  ERR3  +  ERR4  +  ERR 5  +  ERR6 

return 

end 
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PROGRAM  V6PATT 

THIS  PROGRAM  PERFORMS  THE  HILL  CLIMBER  OPTIMIZATION 
ON  THE  SIX  CAPACITIVE  PARAMETERS  OF  THE  PI  NETWORKS  IN  AN  ATT '"KPT 
TO  ENFORCE  THE  DESIRED  CURRENT  DISTRIBUTION  ON  THE  ARRAY  ELEMENTS. 
THIS  PROGRAM  CALLS  CAP6HC . FOR  WHICH  IS  AN  ADAPTED  HILL  CLIMBER 
OPTIMIZATION  SUBROUTINE  THAT  HAS  BE  MODIFIED  TO  ANALYSIS  THE 
EQUILATERAL  ARRAY  DESIGN  PROBLEMS.  THE  NECESSARY  INPUTS  FROM 
THE  USER  ARE  ENTER  IN  THIS  PROGRAM  AND  PLOTTING  OF  THE  RESULTING 
RADIATION  PATTERN  CAN  BE  DONE  BY  CALLING  SUBROUTINE  PATT3. 


REAL  Tight ,OERR, ERR, length , k 1 , what , SWR , rt ,Mi t ,phi t , Pin ,pi , acur , acurlO 

real  scale, indl , ind2, ind3 ,ml imi t ,ml imi tlO , max, null 

real  x( 10) ,dincrx( 10 ) 

real  FAMP(4) ,FPH(4) ,Cap(10) 

real  bestx( 10 ) , amp( 3 ) , ph ( 3 ) , ampC( 3 ) ,phC( 3 ) , LI , L2 , L3 


complex  Zll,Z12,Z13,It,Il,I2,l3, cktl  ( 3) , ckt2 ( 3 ) , ckt 3 ( 3 ) , B ( 1 5 ) 
complex  SI (4 ) , w( 3) , voltl ,Zin,Rfcoe£f , cj , 1 123 ( 3) 

COMPLEX  B123 ( 3 ) 

integer  start , of  ten, ch, wait , wai t2,nvar,nvar2, ixdir ( 10) , idch( 10) , numl 
integer  te, td, t ime .mdist , ihov, icount,dir(10, 40000) 

open(unit»7,file=’v6PATT_one' ,status«'nev' ) 


c 

c  Self  and  Mutual  impedances  values  for  the  3  elements 

c 

zll-(36.5,21.0) 

zl2«(20.4,-14.18) 

zl3=(20.4,-14.18) 

pi«2.0*asin(l.Q) 

cj«(0.Q,1.0) 

c 

c  These  vars.  define  the  dimensions  of  the  3  element  triangular  array 

c 

c  if (Tight .eq.O . 125) length-4 . 0 

c  if (Tight .eq. 0. 25) length-2 . 0 

c  if (Tight. eq.O. 5) length-1.0 

length-4.0 

C  vrite(5,*) ' Input  lenght  of  T.L.' 

C  read( 6 , * )Tlght 

TLGHT-. 125 


C 

C 


vrite(5,*)'  Given  Inductor  Values  (ohms).' 
read (6,*) indl, ind2,  ind3 
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IND1=10 . 0 
IND2-15.0 
IND3-22.0 

vrite(7,*) 'The  Inductive  Reactances  (ohms).’ 

vr  i  te(  7  ,  * ) ’ XL13* ’ , IND1 , 1  XL12= ’ , IND2 , ’  XL23=’,IND3 

write(7,*) *  ’ 

cktl(l)*cmplx(0.0, indl) 
ckt2(l)«cmplx(0.0, ind2) 
ckt3(l)*cmplx(0.0, ind3) 

L1=IND1/(2.0*PI*7.0E6) 

L2*IND2/( 2 . 0*PI *7 . 0E6 ) 

L3«IND3/(2.0*PI*7.0E6) 

WRITE( 7 , * ) ’ At  7  MH2  the  Inductor  Values  are  (henries),’ 

WRITE(7,*) ’Ll-’ ,L1,  ’  L2=  ’  ,  L2 ,  ’  L3*=’,L3 

write(7,*) 

write(7,*) ’ Initial  guess  and  step  size  for  the  Caps  * ,  ip 
write(7,*)’  ’ 

do  ip»l,6 

vrite(5,*) ’ Initial  guess  and  step  size  for  Cap  #’,ip 
read(6,*)X( ip) ,dincrx( ip) 

write (7,*) ’For  Cap#’ , ip, ’Xi** ,x( ip) , ’Si*’ ,dincrx( ip) 
end  do 

vrite(7,*)’  • 

WRITE(5,*) ’ INPUT  LENGTH ( #  OF  ITERATIONS)?’ 

READ(6,*)TIME 

WRITE (7,*) ’Number  of  Iterations  «’,TIME 

WRITE (5,*) ’MAX  STEP  SIZE?’ 

READ ( 6 , *  JMLIMIT 
write(7, *)  ’  ’ 

vrite(7,*) ’MAXIMUM  STEP  SIZE?’ ,mlimit 

vrite(5,*) ’MINIMUM  STEP  SIZE  (acurracy)?’ 

Read(6, *)acur 
WRITE(7 , * ) ’  * 

WRITE(7,*) ’MINIMUM  STEP  IS  ’ ,ACUR 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c  These  are  the  input/given  current  distribution  on  3  elements 

c 

vrite(7,*)  ’  ’ 

wr i t  e ( 7 ,*)’!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!  1  !  !  !  !  I  ’ 
do  ir*l,3 


end  do 


vrite(5,*) ’ Input  Magnitude  and  Phase’ 
read ( 6 , * )amp( ir ) , ph( i r ) 
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write(7, *)' Input  Om  and  On' 
read (6,*) max, null 


WRITE ( 5 , * ) '  MAX* ' ,MAX, '  NULL  ' ,NULL 
write(7,*) '  ' 

Write(7,*) 'Current  Distribution  Corresponds  to  a  Om  and  On  of' 
WRITE ( 7 , * ) ' Om  * ' , MAX, '  On  -’.NULL 


Change  amp  and  phase  to  real  and  imaginary  component 


CALL  CMPN ( 3 , AMP , PH , W ) 


Change  complex  current  values  to  normal  amplitudes  and  phase 


CALL  MAG_PH_NORM(3,W,AMP,PH> 


CALL  CMPN ( 3 , AMP , PH , W ) 


wr i te( 7 , * )  '  ' 

write(7,*) ' Amps  Normalized  * ' , amp( 1 ) , amp( 2 ) , amp( 3 ) 

write(7,*) 'Phase  (Deg.)  *' ,ph(l)*180. 0/pi, pb(2)*180. 0/pi, ph(3)*180. /pi 


Calculate  It  (input  current  into  array)  using  normalized  given  curr. 
write(7,*)'  ' 

It«(w(l)*Zll  +  v( 2 ) *Z12  +  w(3)*Z13)/(50. 0,0.0) 
write(7,*) ' It  (Desired  input  current)  *',It 


CALL  CAP6HC(cktl,ckt2,ckt3,PI,W, IT,X,DINCRX,TIME,MLIMIT,ACUR,S1 ) 


CALL  MAG_PH_NORM ( 4 , SI , FAMP , FPH ) 

rt*180.0/pi 

write (5,*) '  I* ' ,  I 
write(7,*) '  I-' , I 
write(7,*) ’  ’ 

write (7,*) ’AMPS*' ,Famp(l) ,Famp(2) ,Famp<3) , ' ITM* ' ,FAMP<4) 

write (7,*) ' PHASE*’ ,Fph( 1 ) *rt , Fph( 2) *rt , Fph( 3)*rt , ’ ITPH-’ ,FPH(4)*RT 

write(7,*) ’  ' 

write(7,*) 'Best  Capacitative  Reactances' 

write(7,*)x(l),x(2),x(3),x(4) 
wr  ite(7 , * ) x( 5 ) ,x(6) 

WRITE( 7 , * ) '  ' 

write(7, *) 'Capacitance  of  elements  at  7  MHz' 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c  Calculates  the  Capacitance  of  L.C.  Netwirk  Parameters 

c 


do  ip=l,nva.- 

cap( ip)=1.0/(x( ip)*2.0*pi*7.0E6) 

write(7,*)'  Cap.  of  ',ip,'  is ' , cap( ip) , '  Farads' 

end  do 


voltl=Sl(l)*zll  +  SI { 2 ) *zl2  +  SI ( 3 ) *zl3 
Zin=voltl/Sl(4) 

Rfcoef f=(Zin-( 50.0, 0.0 ) )/(Zin+< 50. 0,0.0) ) 
SWR=(l+cabs(Rfcoeff ) )/(1.0  -  cabs(rfcoef f ) ) 

write(7, *) '  ' 

write (7,*) 'Zin=' ,Zin, '  SWR-'.SWR 
write (7,*) '  Itbest“' ,S1(4) , 'Vl«' ,voltl 

wr i te ( 5 ,  * ) ' Zin  found*=',Zin 
write(5,*) '  SWR»',SWR 

Calculate  values  for  It  for  singlr  dipole  problem 


Pin«( ( (CABS(S1(4) ) )**2.0)*REAL(ZIN) ) 

scale«sqrt(1000.0/Pin) 

mit«sqrt( ( 1000 . 0 ) /36 . 5 ) 

phit *0.0 


numl*=l 

call  patt3 (scale, Famp, Fph, pi , length, numl , mi t, phit ) 

numl-=2 

Famp(l)=1.0 

do  ip=l,3 

Famp( ip) «Famp( ip) *scale 

end  do 

call  patt 3 ( scale, Famp, Fph, pi , length, numl , mi t, phit ) 


CLOSE ( 7 ) 

STOP 

END 


INCLUDE  'CAP6HC.FOR' 
INCLUDE  ' FUNC6C. FOR' 
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c 

c 

c 

c 

c 

SUBROUTINE  CAP6HC(cktl , ckt2,ckt3,pi,v, It ,x, diner*, time, mlimit , acur ,S1 ) 

REAL  OERR, ERR, what, acur, acur 10 
real  mlimit,mlimitlO,Merr(40000),Ediff 
real  arrx(10,40000),x(10),dx(10),dincrx(10) 
real  xsml(lO) ,xlrg(10) 
real  bestx(lO) ,pi,BSTIN(3) 

complex  B(15),It,w(3),Sl(4) ,cktl ( 3 ) , ckt2( 3 ) ,ckt3 ( 3 ) 
integer  ch,wait,wait2,nvar,nvar2, ixdir(lO) , idch(10) 
integer  te,td,time,mdist, ihow 
IHOW=time-l 


ch=l 

oerr=l . 0E30 
1=0 
K=0 
J=0 
L=0 
M=0 

mlimitl0=mlimit*10 . 0 
acurl0=acur/10. 0 


nvar=6 

mdist=nvar*10 

nvar2=nvar*2 

vait=nvar*4 

vait2=nvar*2+l 


Initialize  points  for  Hclimber  vars.  and  increments 


do  io=l,nvar 

dx( io)=1.0 
ixdir ( io) =1 
xsml( io)=0.0 
xlrg ( io) =0 . 0 
idch( io)=0 

end  do 


$1 


c 

c 


This  begins  the  Hclimber  by  calling  fuction  to  be  optimized 


100  Call  Func6c(pi,v,x, i t , b, cktl ,ckt2 , ckt3 , err ,S1 ) 


c  This  is  iteration  counter  used  in  Hclimber 


1*1  +  1 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c  This  block  takes  the  first  initial  step  in  a  positive  incrmental 

c  direction  for  X(l). 

c 

I F ( I . EQ . 1 ) THEN 

x(l)=x(l)  ♦  dx(l)*dincrx(l) 
vrite(7,*)'  ' 

WRITE(7,*) ' Initial  Hill  Height  =',err 
if (err.lt.oerr)then 
oerr=err 

end  if 
GOTO  100 

END  IF 

cccccccccccccccccccccccccccccccccccccccccccccccccc 

Merr ( I ) *err 

Edit f*abs(Merr( I )-Merr ( 1-1) ) 

OEdif f*abs(err-oerr ) 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c  Printing  option  that  prints  out  error, I,  and  Ntwk.  elements  as  desired 

c 

if ( I .gt , Ihow) then 
vrite(7,*)'  ' 

vrite(7,*) 'Number  of  Iterations, HHi, Variable  Seq. ,  &  HH( i)-HH( i-1) ' 
vrite(7,*) 

WRITE(7 , 15)I,err,ch,ediff 

15  Formati/'  I*  ',I5,2x,'HH*  ’ , E12 . 6 , 2x , 'Ch= ' , 12 , 2x, ' HHdi f- ' , E12 . 6 ) 
wr i te{ 7 , * ) ' Var ibles  present  value 
write(7, *) '  ' 

vrite(7,17)x(l),x<2),x(3),x(4) 

17  Format (6x, 'XI*' ,F12.6,2x, ' x2= ' , F12 . 6 , 2x, ' x3« ' ,F12 . 6 , 2x, ' x4* ' , F12 . 6 ) 

write(7,») '  ' 

write(7,18)x(5) ,x(6) 

18  Format (6x, 'X5=' ,F12.6,2x, 'x6*' ,F12.6) 

write(7,*) '  ' 

vrite(7,*)'  Step  sizes  . ' 

wr i te( 7 , * ) ' 

write(7,16)dincrx( 1) ,dincrx(2) ,dincrx(3) ,dincrx(4) 

Format ( 'si* ' , F12. 6, lx, ' s2= ' , F12 , 6, lx, ' s3« ’ , F12. 6, lx 
$  , 's4«' ,F12.6) 
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write(7,21)dincrx(5),dincrx(6) 

21  Format ( 's5=',F12.6,lx, 1 s6=* ,F12.6) 

write(7,*) 1  ’ 

vrite(7,*)'  Direction  . ' 

write(7,*) ' 

write(7,19)dx(l) ,dx(2) ,dx(3) ,dx(4) 
lb  Format ( 'dl=' , F4 . 1 , lx, ' d2= ' , F4 . 1 , lx , 'd3= ’ , F4 . 1 , lx , ’ x4= ' , F4 . 1 ) 

write(7,*)'  ' 

endi  f 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c  This  loop  loads  the  storage  array  with  all  values  of  variables  through 

c  out  the  iteration  sequence, 

c 

do  ix=l,nvar 


arrx( ix, I )=x( ix) 


end  do 


These  two  do  loops  determine  determine  when  a  variable  needs  to  be 
reduced  in  size  when  it  has  gotten  relatively  close  to  an  answer 
and  begins  jumping  around  it. 


IF( I .GE.wait)then 


do  1  ia«l,nvar 
j  =  l+l 

k=l  +  wait2 

m-0 

mm=0 

do  2  ib*l,nvar2 


if (arrx(ia, j ) .eq.arrx( ia,k) )m«m+l 

j»j  +  l 

k*k+l 

continue 

if (l.le. idch( ia) )goto  1 

if (m.eq. nvar2 (then 

idch(ia)=l  +  wait2 

x( ia) «bestx( ia) 

dincrx(ia)»dincrx(ia)/10.0 

wr ite( 7 , * ) ' *********** — variable  is' ,  ia 

write{7,*) 'Decreased  increment  of  X,dincrx*’ ,dincrx(  ia) 

write(7, *) 'X=  ' ,X( ia) 


204 


end  if 


continue 
1*1  +  1 
endi  f 


This  loop  keeps  the  acurracy  of  the  variables  from  going  out  01  limit 


do  if=l,nvar 

if (dincrx( if ) . le.acurlO) then 
dincrx( if )=acur 

write(7,*) 1 - ’ 

write(7,») ’dincrxl  has  gone  to  low,  increased  it  to’ .dincrxl 

endif 
end  do 

This  block  of  code  determines  when  to  increases  the  size  (mag.)  of  the 
variable  based  on  how  far  it  has  climber  in  one  direction. 


do  3  ic»l,nvar 

if (ixdir(ic) . eq. 1) then 

xlrg( ic)-xlrg( ic)+l 

else 

xlrg( ic)-0 

endif 


if(ixdir(ic) .eq.OJthen 

xsml( ic)*xsml( ic)+l 

else 

xsml( ic)=0 

endif 

if (xlrg(ic) . eq.mdist)then 

diner x( ic ) -diner x( ic)*10.0 

Write(7, *) '############ — variable  is' ,  ic 

Write( 7 ,*)' Increased  increment  of  X,dincrx=' ,dincrx( ic) 


endif 


xlrg( ic) *0 


if (xsml ( ic) . eq.mdist ) then 

dincrx( ic) -dincrx ( ic)* 10.0 

write(7,*) '############  --variable  is',ic 

wr i te( 7 , * ) ' Increaesd  increment  of  X, dincrx- ' , dincrx ( ic ) 
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xsml ( ic ) =0 

end  if 

3  continue 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCD-CCCCCCCCCCC 

do  if=l,nvar 

if (dincrx( if ) .ge.ml imi tlO ) then 
dincrx( if )=mlimit 

c  vr ite( 7 , * } ' STEP  SIZE  GONE  ABOVE  LIMIT’ 

c  vrite(7,*)'  ' 

endi  f 

end  do 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c  This  block  of  code  changes  the  variable  back  to  its  previous  value 

c  if  changing  it  caused  the  error  to  increase.  This  block  also  increases 

c  or  decreases  the  variable  based  on  the  current  direction  of  the  var. ,  by 

c  the  value  in  incrx(?). 

c 


do  4  id=l,nvar 

if (ch.eq.nvar)then 

I F { ED I FF . eq . 0 . 0 ) THEN 

x(nvar )»x(nvar )-dx(nvar ) *dincrx(nvar ) 
x  ( 1 ) =x ( 1 ) +dx ( 1 ) *di ncr x ( 1 ) 

D I NCRX (NVAR) *D I NCRX ( NVAR) *10.0 

c  WRITE(7,*) 'ERR  -  PREVIOUS  ERR  LT  IE-6 , Bx= ’ , x ( nvar ) 

c  WRITE { 7 , * ) ' 1= ' , I , '  id=',id,'  step*' .dincrx(nvar) 

goto  70 
END  IF 


I F ( OED I FF . eq . 0 , 0 ) THEN 

x { nvar ) »x ( nvar ) -dx ( nva r ) *di ncr x( nvar) 
x ( 1 ) =x( 1 )+dx( 1 ) *dincrx( 1 ) 

D I NCRX ( NVAR ) =D I NCRX ( NVAR ) *  1 0 . 0 

c  WRITE(7,*) 'Oerr  and  Err  are  =  to  e-6’ 

c  wr i t  e ( 7 , * ) ' EEEERRRRRRRRRRR  =  OOOOeeerrrr’ 

c  WRITE( 7 , * ) ' 1  =  ' , I , '  id=',id,'  step*’ , dincrx ( nvar ) 

goto  70 


END  IF 


u  u 


if (err.lt.oerr)x(l)=x(l)+dx(l)*dincrx(l) 
if(err.gt.oerr)then 

x(nvar )=x(nvar)-dx(nvar )*d incrx (nvar) 
x ( 1 ) =x ( _ ) +dx ( 1 ) *d incrx ( 1 ) 

end  if 

goto  70 
end  if 

IF(CH.EQ.ID  .AND.  EDI FF . eq . 0 . 0 )THEN 
TE=ID+1 

x( id) *x( id) ~dx( id) *dincrx( id) 
x( te ) =x( te )+dx(te)*d incrx (te) 

DINCRX ( ID)=DINCRX( ID) *10.0 

WRITE(7, *) 'ERR  -  PREVIOUS  ERR  LT  IE-6 r Bx= ’ , x { id) 
WRITE(7, *) ' I*' , I , '  id= ' , id, '  step*’ ,dincrx(nvar) 

goto  70 
END  IF 


I F ( CH . EQ .ID  .AND.  OEDIFF. eq. 0 . 0 )THEN 
TE-ID+1 

x( id)-x( id)~dx( id)*dincrx( id) 
x(te)*x( te)+dx(te)*dincrx(te) 

DINCRX( ID) “DINCRX ( ID)*10.0 

WRITE(7,*) 'Oerr  and  Err  are  *  to  e-6’ 
wr i t e ( 7 , * ) 1 EEEERRRRRRRRRRR  *  OOOOeeerrrr’ 

WRITE ( 7 ,  * ) 1 1  = ' ,  I , '  id*', id,'  step* ' ,dincrx(nvar) 

goto  70  ■ 

END  IF 


IF(ch.eq. id  .and.  err. lt.oerr)then 
td-id+1 

c  if ( id.eq.nvar)td=l 

x(td)*x( td)+dx( td)*d incrx (td) 

endif 


I f (ch. eq. id  .and.  err.gt.oerr)then 
td=id+l 

x( id)=x( id)-dx( id) *d incrx ( id) 
x(  td)  =x{  td)+dx( td) *d incrx ( td) 

endif 


«/  »/  » / 
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4  continue 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c  This  loop  changes  the  variables  direction  basd  on  the  value  of  error 

c 


70  do  5  ie=l,nvar 

if(ch.eq.ie  .and.  ixdir( ie) .eq. 1  .and.  err.gt.oerr)then 
ixdir ( ie)*0 
dx( ie)=-1.0 
goto  20 

end  if 
c 
c 

iffch.eq.ie  .and.  ixdir( ie) .eq.O  .and.  err.gt.oerr)then 
ixdir ( ie) =1 
dx( ie)“l,0 
goto  20 

end  if 

5  continue 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c  This  block  keeps  track  of  which  variable  is  to  be  changed  next, 

c 

20  if(ch.le.nvar)then 

if (ch.eq.nvar )ch=0 
ch=ch+l 

end  if 


c 

c  This  block  records  value  of  the  best  variables  and  current  values 

c  as  the  Hclimber  proceeds  thru  the  climb, 

c 

if(err.lt.oerr)then 

oerr«err 

do  ip«l,nvar 

bestx(ip)«arrx(ip,I) 

end  do 

endi  f 
c 

c  This  tells  machine  to  do  another  iteration 

c 

I F ( I .LT. time) GOTO  100 
do  iz»l,nvar 

x( iz)=bestx( iz) 
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c 

C  SUBROUTINE  FUNC6C 

c 

C  THIS  SUBROUTINE  PERFORMS  THE  ANALYSIS  OF  THE  EQUILATERAL 

C  ARRAY  BY  CALLING  THE  ANLYSUBR  SUBROUTINE.  THE  SIX  CAPACITIVE 

C  REACTANCES  BEING  ADJUSTED  BY  THE  HILL  CLIMBER  OPTIMIZATION 

C  ROUTINE  ARE  SENT  TO  ANLYSUBR  FOR  CALCULATION  OF  THE  RESULTING 

C  CURRENT  DISTRIBUTION.  THE  HILL  HEIGHT  TS  THEN  CALCULATED  AND 

C  RETURNED  TO  THE  OPTIMIZATION  ROUTINE  FOR  EVALUATION.  THE  DESIRED 

C  CURRENTS  ARE  INPUT  AS  VARIABLE  W  AND  THE  PI  NETWORK  PARAMETERS 

C  ARE  IN  VECTORS  CKT1 , CKT2 , CKT3 . 

C 

C  INPUTS:  PI , W, X , IT , CKT1 ( 1 ) , CKT2 ( 1 ) , CKT3 ( 1 ) 

C  OUTPUT:  ERR, SI , B 


SUBROUT I NE  FUNC 6C (PI,W,X,IT,B, CKT 1 , CKT 2,CKT3,ERR,S1) 


REAL  err , x ( 10 ) ,Pi , AMPCN(4 ) ,PHCN( 4 ), Tight 

COMPLEX  cktl ( 3 ) , Ckt2 ( 3 ) , Ckt3 (3) ,B(15) ,  Si ( 4 ) ,v(3) , It,B123(4) 


THIS  VARIABLE  IS  THE  LENGTH  OF  THE  TRANSMISSION  LINE 
BETWEEN  THE  ARRAY  ELEMENTS  AND  THE  PI  NETWORKS. 


TLGHT= .125 


THE  PRESENT  VALUES  FOR  THE  SIX  CAPACITIVE  REACTANCES  ARE  LOADED 
INTO  THE  VARIABLES  TO  BE  PASSED  INTO  THE  ANALYSIS  ROUTINE. 


cktl(2)»cmplx(0.0,x(l) ) 
ckt2(2)»cmplx(0.0,x(2) ) 
ckt3(2)=cmplx(0.0,x(3) ) 

cktl(3)»cmplx(0.0,x(4) ) 
ckt2 ( 3 ) -cmplx ( 0 . 0 , x ( 5 ) ) 
ckt3(3)=cmplx(0.0,x(6) ) 


This  is  a  call  to  subr.  curr  in  file  anlysubr.for  which  solves 
the  analysis  of  the  array. 


call  curr (pi .Tight , cktl , ckt2 , ckt3 , i t , b) 


This  loops  detects  when  the  variables  have  gone  below  zero  and 
adds  the  square  of  the  variables  magnitude  to  the  hill  height. 


err8=0 . 0 
do  ip=l,nvar 


n  n 
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c  if (x( ip) .lt.O.O)then 

c 

c  err8*err8  +  x( ip) **2 

c  write(7,*) ' x=' , ip, '  gone  LT.  O.O'.x(ip) 

c 

c  end if 

c  end  do 

c 

c  This  loop  is  a  warning  block  of  code  which  detects  when  the 

c  currents  for  the  array  elements  found  by  the  analysis  routine 

c  are  zero  (undetermined).  This  means  that  there  is  no  solution 

c  for  the  present  values  of  the  network  parameters, 

c 
c 

do  ip*l,3 

if (real(b( ip) ) .eq.0.0  .and.  aimag(b( ip) ) .eq.0.0)then 
•write(7, *) ' 

wr i te ( 7 ,*)' WARNING  CURR' , ip, '  ....  to  ZEROES' 

wr i te ( 5 , * ) ' WARNING  CURR', ip,'  _  to  ZEROES' 

writeO,*)1 

endif 

bl23 ( ip)=b( ip) 

end  do 

bl23 ( 4 ) «b( 1 )+b( 4 ) +b( 6 ) 
c 

c  Changes  complex  current  distributions  in  normalized  magnitude 

c  and  phase  vectors, 

c 

CALL  MAG_PH_N0RM ( 4  ,  B 1 2  3  , AMPCN , PKCN ) 

Changes  notmaiized  magnitude  and  phase  vectors  into  complex  values, 
c 

CALL  CMPN ( 4 , AMPCN, PHCN, SI ) 
c 

c  This  is  the  calculation  of  the  hill  height,  which  is  the 

c  difference  between  the  desired  and  calculated  values, 

c 

errl=( real(w( 1) )-real(sl( 1) ) ) **2.  +  ( real(w(2) )-real(sl<  2) ) ) **2 . 
err2=(real(w(3))-real(sl(3)))**2 . +(aimag(w( 1 ) )-aimag (sl(2))>**2. 
err3=(aimag(w(2) )-aimag(sl(2) ) )**2. 
err4=(aimag(w(3) ) -aimag ( si { 3 ) ) )**2. 

er r*errl+err2+er r3+er r4  +  err8 


RETURN 

END 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  SUBROUT IKE  MAG  PH_NORM 

C 

C  This  subroutine  is  used  to  convert  a  complex  number 

C  into  Phasor  form  (magnitude  and  phase).  Once  in 

C  phasor  form  the  magnitude  and  phases  are  then  normalized 

C  with  respect  to  the  current  flowing  on  the  first  antenna. 

C  That  is  the  first  complex  number  in  cmplxl  which  is  cmplxl(l). 

C 

C 

C  Input  s  V, CMPLXI 

C 

C  Output  :  MAGTN, PHSEN 

C 

C 

SUBROUTINE  MAG_PH_NORM( V , CMPLXI , MAGTN , PHSEN) 

INTEGER  A, V 

REAL  MAGTN ( V ) , PHSEN ( V ) 

COMPLEX  CMPLXI (V) 

C 

C  CALCULATIONS  OF  THE  MAGNITUDE  AND  PHASE  OF  A  INPUT  COMPLEX  NUMBER. 

C 

DO  A«1,V 

MAGTN ( A ) -CABS ( CMPLX 1(A)) 

PHSEN ( A) «ATAN2(AIMAG< CMPLXI (A) ) , REAL (CMPLXI (A) ) ) 

END  DO 
C 

C  NORMALIZATION  ACCOMPLISHED  BY  DIVISION  OF  THE  MAGNITUDES 

C  AND  SUBTRACTION  OF  THE  PHASES. 

C 


DO  A-2 , V 


MAGTN ( A ) -MAGTN ( A ) /MAGTN ( 1 ) 
PHSEN( A ) -PHSEN ( A ) -PHSEN ( 1 ) 


END  DO 

MAGTN ( 1 ) - 1 . 0 
PHSEN ( 1 ) -0 . 0 


RETURN 

END 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTINE  CMPN 


THIS  SUBROUTINE  WAS  DESIGNED  TO  CALCULATED  THE  COMPLEX 
VALUE  FROM  INPUT  MAGNITUDE  AND  PHASE  VECTORS. 

GIVEN  THE  MAGNITUDE  AND  PHASE  OF  A  VECTOR  THE  COMPLEX 
NUMBER  REPRESENTATION  IS  FOUND  USING  THE  FORMULA 

M*COS ( 0 )  +  j  *M*SIN( 0 ) 

INPUTS  :  V,MAGN,PHSE 
OUTPUT  :  CMPLXN 


SUBROUTINE  CMPN(V,MAGN, PHSE , CMPLXN ) 

INTEGER  A, V 

REAL  MAGN ( V ) , PHSE ( V ) 

COMPLEX  CMPLXN {V ) 


DO  A-1,V 

CMPLXN(A) *MAGN( A) *COS( PHSE( A) )  +  ( 0 . 0 , 1 . 0 ) *MAGN( A) *SIN( PHSE ( A) ) 

END  DO 

RETURN 

END 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  SUBROUTINE  CURR 

c 

C  THIS  SUBROUTINE  PERFORMS  THE  EVALUATION 

C  OF  THE  15  LINEAR  COMPLEX  EQUATIONS  THAT  DESRCIBE  THE 

C  EQUILATERAL  THREE  ELEMENT  ARRAY.  THE  PI  NETWORK  PARAMETERS 

C  (3  CAPACITORS  AND  6  INDUCTORS)  ARE  INPUT  INTO  THIS  ROUTINE 

C  AND  ARE  CONSTANTS  PROVIDED  BY  THE  USER.  THE  PARTICULAR 

C  MAGNITUDE  OF  THESE  PARAMETERS  AND  THE  LENGTH  OF  THE 

C  TRANSMISSION  LINES  ARE  WHAT  DETERMINE  THE  COEFFICENTS 

C  OF  THE  15  UNKNOWN  VARIABLES.  THIS  ROUTINE  USES  A  LINEAR 

C  EQUATION  SOLVER  ACCESSED  FROM  THE  IMSL  LIBRARY  OF  FORTRAN 

C  SUBROUTINES. 


i 

I 


i 

I 


to, 


subroutine  CURR(pi .Tight , cktl , CKT2 , CKT3 , I t , B) 


declare  al-1  variables- 


Sj! 

v 

k 

$ 

'f!i 


INTEGER  N, IA,M, IB, I JOB, IER 

COMPLEX  A( 15, 15) , B< 15) , WA( 255 ) , CKT1 ( 3 ) , CKT2 ( 3 ) ,CKT3 ( 3 ) , XL13 , XC13 
COMPLEX  XC31,Z11,Z12,Z13,IT,C46,C02,PII,X1,X2,X3,X4,X5,X6 
REAL  wk ( 15) .Tight , pi , L46, L02, LI  I 


-matrix  parameters- 


IA«15 
IB-15 
N«15 
M=  1 


-mutual  impedances- 


Zll=( 36 . 5,21.0) 
Z13= (20.4,-14.18) 
Z12=Z13 


-calculation  of  constants  based  on  transmission  line- 


L4 6=50. 0*(sin(3.0*pi ‘Tight) ) 
L02=(sin(3.0*pi ‘Tight) )/50.0 
LII=cos(3.0*pi ‘Tight) 

C46=cmplx( 0 . 0 , L46 ) 

C02=cmplx( 0 . 0 , L02 ) 

PI  I=cmplx(LI 1,0.0) 


m 


c 

c 


- set  constant  matrix  b  back  to  zero - 


A(7,7)=(l.,0.) 
A( 7 , 10 ) =X4 
A( 7 , 11 ) =X6 


» 


SECOND  SET  OF  EQUATIONS 


XL13=CKT2 ( 1 ) 
XC13=CKT2(2) 
XC31=CKT2 ( 3 ) 


CALL  COEFF(C46 , C02 , PI  I , XL1 3 , XC13 , XC3 1 , XI , X2 , X3 , X4 , X5 , X6 ) 


- LOAD  COEFFS  INTO  MATRIX  "A" 

A(8 , 1 ) =Z11 
A(B,2)=Z12 
A( 8 , 3 ) =Z13 
A( 8 , 12 ) =X1 
A(  8 , 13 ) =X2 


A(9,4)»(l.,0.) 
A( 9 , 12 ) =X3 
A(9,13)=X4 


A( 10 , 1 ) «Z12 
A(10,2)«Z11 
A( 10, 3 )=Z13 
A( 10 , 12 ) “X2 
A(10,13)«X5 


A(ll,5)-(1.,0.) 
A( 11 , 12 ) =X4 
A( 11 , 13 ) =X6 


THIRD  SET  OF  EQUATIONS 


XL13=CKT3 ( 1 ) 
XC13=CKT3 ( 2 ) 
XC31=CKT3(3) 


CALL  COEFF { C46 , C02 , PI  I , XL13 , XC1 3 , XC31 , XI ,X2 ,X3 , X4 , X5 ,X6 ) 


A(  12, 1)=Z12 
A(  12 , 2)  «=Z11 
A( 12 , 3 ) =Z13 
A( 12 , 14 ) =X1 


a 


A(  12 , 15 ) =X2 


A(  13 ,8 )  =  ( 1. , 0. ) 
A(  13 , 14 ) =X3 
A(13,15)=X4 


!»v* 

% 

lift 


A(  14 , 1 )=Z12 
A(  14 , 2 ) =Z13 
A(  14 , 3  >=Z11 
A( 14 , 14 ) =X2 
A(  14 , 15 ) =X5 


A(15,9)=(l. ,0. ) 
A( 15 , 14 ) =X4 
A( 15 , 15 ) =X6 


- SET  FUNCTION  DEFINITION- 


I JOB=0 


-CALL  MATRIX  INVERSION  IMSL  ROUTINE- 


CALL  LEQ2C( A,N, I A , B , M , IB, IJOB,WA,WK, IER) 


RETURN 

END 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  SUBROUITNE  COEFF 

C 

C  THIS  ROUTINE  CALCULATED  THE  COEFFICENTS  OF  THE 

C  UNKNOWN  VARIABLES  BASED  ON  THE  LENGTH  OF  THE  TRANSMISSION 

C  LINE  BETWEEN  THE  ARRAY  ELEMENT  AND  THE  PI  NETWORK  AND  THE 

C  MAGNITUDE  OF  THE  NETWORK  PARAMETRS 


INPUT:  C46 ,C02 , PI  I ,XL13 ,XC13 , XC31 
OUTPUT:  X1,X2,X3,X4,X5,X6 


SUBROUTINE  COEFF(C46 , C02 , PI  I , XL1 3 , XC1 3 , XC31 , XI , X2 , X3 , X4 , X5 , X6 ) 


-DECLARE  COMPLEX  VARIABLES- 


non 


»n*  rj<  nj<  njt  -jm JfR  1.VUVIM 


COMPLEX  T1 ,T2 ,X1 , X2 , X3 ,X4 ,X5 ,X6,C46 ,C02 , PI  I , XL13 ,XC13 , XC31 


-COMPLEX  ALGEBRA  CALCULATIONS- 


T1=(PII+(C46/XL13) ) 
T2= (C02+(PI I/XL13 ) ) 
Xl=- (T1+ (C46*XC13 ) ) 
X2= (C46/XL13 ) 

X3=- (T2+ (PI  I *XC13 ) ) 
X4= (PI I/XL13 ) 

X5=- (T1+ (C46*XC31 ) ) 
X6=- (T2+ (PI I*XC3 1 ) ) 


RETURN 

END 
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CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
SUBROUTINE  PATT3 

This  subroutine  calculates  the  radiation  pattern 
for  the  three  element  equilateral  triangular  array  used 
to  demonstated  the  design  methods  of  this  thesis. 

The  amplitudes  and  phases  of  the  elements  of  the  array 
are  input  from  one  of  the  design  procedures.  This  routine 
is  called  by  one  of  the  design  procedures  to  display  the 
pattern  generated  by  the  solved  current  distributions. 

INPUT  :  SCALE, AMP, PH, PI , LENGTH, NUM1, MIT, PHIT 

OUTPUT:  XY 


subroutine  patt3 ( scale, amp, ph, pi , length, numl ,Mi t , phi t ) 

real  amp(3) ,p(3),ph(3),x(3},y(3),xy(360,2),Etot(360) 
real  scale.d, length, kl, pi , Emax,mi t ,phi t , er ( 3 ) , ei ( 3) 
integer  numl 


These  are  the  coordinates  of  the  element  on  the  XY  plane 
The  length  variable  is  the  reciprocal  of  the  separation 
distance  between  the  elements.  KL  is  used  to  convert 
the  geometric  phase  shift  to  radians. 


kl=2.0*pi/length 

x(l)=0.0 

x( 2 ) =0 . 5 

x(3)=1.0 

y(l)=0.0 

y(2)». 866025404 

y(3 )=0. 0 


c 

c  Use  to  plot  the  pattern  of  a  single  dipole, 

c 


if (numl.eq.l)then 
n«  1 

amp(l)=Mit 
p< 1 >«phi t 
x ( 1 ) *0 . 0 
y(l)=0.0 

else 

n=3 

end  if 


Emax=0 ,0 
do  1=1,360 

Theta= float (I)*pi/180.Q 


do  J=l,n 

if (x(J) .eq.0.0  .and.  y(J) .eq.0.0)goto  14 


P  is  the  geometric  phase  shift  associated  with  the  physical 
location  of  the  elements  with  respect  to  each  other. 


P( J ) =kl*sqrt (X(J)**2.+Y(J)**2.) *cos (atan2(Y(J) ,X(J) ) -Theta ) 
p(l)=0.0 


Er(J)  is  the  real  component  of  the  complex  current  on 
the  element  and  Ei(J)  is  the  imaginary. 


er(J)=amp(J)*cos(ph(J)+P(J) ) 
ei(J)=amp(J)*sin(ph(J)+P(j) ) 

end  do 


This  is  for  the  single  dipole  calculations. 


if(numl.gt.l)then 
do  J=2 , 3 

er(l)=er(l)+er(J) 

ei(l)=ei(l)+ei(J) 

end  do 

end  if 


Etot(I)  is  the  magnitude  of  the  radiating  current 
as  seen  at  all  I  angles  of  observation. 


Etot (I)=Sqrt(er(l)**2.+ei(l)**2.  ) 
if (Etot ( I ) .gt .Emax )Emax=Etot ( I ) 


end  do 

write (7,*)’  Emax=',Emax 


This  360  degrees  loop  calculates  the  XY  coordinates  of  the 
radiation  pattern  of  the  array. 


do  1=1,360 

theta= float (I) *p i/180.0 
ETOT ( I )  =  ( Etot ( I ) *4 1 . 0 ) 
xy ( I , 1 ) =ETOT( I ) *cos( theta ) 
xy(I,2)=ETOT( I )*s in (theta) 

end  do 


Calls  the  plotting  routine  which  display  the  pattern  to  the  screeen 
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This  subroutine  uses  the  PlotlO  graphics  library 
on  the  Vax  11/780  at  the  University  of  Louisville. 

It  is  used  to  plot  radiation  patterns.  A  matrix 
containing  the  X  an..  Y  coordinates  of  the  pattern 
is  passed  into  this  routine  and  move  and  draw  commands 
are  used  to  draw  this  pattern. 

Input  :  XY,NUM1 


subroutine  plot ( xy , numl ) 

dimension  xy(360,2) 
real  radd, ang , x, y,pi 
integer  numl 


call  grstrt (4014 , 1 ) 

i  f ( numl . eq . 1 ) then 

call  newpag 

end  if 

call  window( -750. 0,750. 0,-750. 0,750.0) 

call  move(-750. 0,0.0) 
call  draw(750. 0,0.0) 
call  move(0. 0,750.0) 
call  draw(0. 0,-750.0) 

call  move ( xy ( 1 , 1 ) , xy ( 1 , 2 ) ) 

do  1=2,360 

x=xy( 1,1) 
y=xy(  T,2) 
call  draw(x,y) 

end  do 

call  grstop 

return 

end 
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